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ABSTRACT 

of 

FILM  AND  TRANSPIRATION  COOLING 
of 

NOZZLE  THROATS 

Analytical  studies  vere  made  of  liquid  film  cooling,  gas  film 
cooling,  and  liquid  or  gas  transpiration  cooling  of  hypersonic  nozzles. 
Experimental  studies  of  gas  film  cooling  were  performed  on  a  small 
nozzle  using  air  as  both  the  mainstream  and  coolant  gases.  Based  on 
the  experimental  results  obtained  and  the  development  of  satisfactory 
calculations!  techniques  to  implement  the  analyses  given  in  the 
present  study,  the  following  conclusions  were  drawn: 

1.  Gas  film  cooling  can  be  used  to  measurably  lower  the  wall 
temperatures  and  the  wall  heat  fluxes  in  the  converging  section  and 
at  the  throat  of  a  high-pressure  high-temperature  nozzle.  Some  gross 
mixing  occurs  at  the  injection  point,  the  exact  amount  being  some  as 
yet  undetermined  function  of  the  injection  geometry,  the  relative 
velocities  of  the  main  gas  stream  and  the  coolant  stream  at  the  in¬ 
jection  point,  and  the  entering  velocity  profiles  and  turbulence 
conditions. 

2.  A  straightforward  boundary  layer  type  analysis  was 
developed  and  programmed  which  predicts  with  reasonable  accuracy  the 
nozzle  wall  temperatures  and  wall  heat  fluxes  in  the  converging 
section  and  at  the  throat  for  gas  film  cooled  nozzles. 

3.  Calculational  techniques  have  been  developed  and  programmed 
for  predicting  the  effectiveness  of  liquid  film  cooling  and  liquid  or 
gas  transpiration  cooling  in  nozzles.  These  techniques  should  be 
checked  against  experimental  data  before  they  are  used  for  design 
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purposes ,  however. 
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CHAPTER  I 
INTRODUCTION 

I.  GENERAL  DISCUSSION 

Recent  technological  developments  have  caused  an  Increased  Interest 
in  analytical  and  experimental  studies  of  film  and  transpiration  cooling. 
Military  rocket  and  space  research  has  brought  about  extensive  investi¬ 
gations  into  the  problems  of  rocket  nozzle  design.  In  particular , 
greater  specific  thrust  can  be  achieved  by  using  higher  combustion  chamber 
pressures  and  temperatures,  however  this  requires  better  nozzle  throat 
cooling  in  order  to  reduce  corrosion,  oxidation,  and  thermal  stresses. 
Another  problem  involves  the  surfaces  of  re-entry  vehicles*  which  must  be 
protected  from  the  extremely  high  heat  loads  encountered.  Experimental 
studies  utilizing  wind  tunnels  which  simulate  re-entry  conditions  provide 
the  most  direct  approach  to  this  problem.  Wind  tunnel  systems  designed 
to  provide  these  test  conditions  must  use  air  at  extremely  high  stagnation 
pressures  and  temperatures  resulting  in  extremely  high  heat  loads  at  the 
nozzle  throat  which  must  be  protected  by  proper  heat  transfer  design. 

Examination  of  ordinary  backside  cooling  of  the1  nozzle  throat 
shows  that  it  may  become  inadequate  for  two  principal  reasons: 

(1)  Exceedingly  high  static  pressure  of  the  coolant  may  be  required  to 
prevent  boiling  (for  conventional  liquid  coolants)  with  the  subsequent 
danger  of  burnout,  (2)  backside  cooling  cannot  reduce  the  heat  load  to 
the  nozzle  surface  so  it  cannot  alleviate  the  problem  of  thermal 
stresses. 
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It  has  bees  proposed,  therefore,  that  rocket  combustion  chambers 
and  nozzles,  re-entry  vehicles,  high-pressure  high-temperature  wind 
tunnel  nozzle*  and  other  device*  with  the  same  problems  be  protected 
by  film  or  transpiration  cooling  or  a  combination  of  one  of  these  with 
backside  cooling,  Film  cooling  refers  to  the  injection  of  a  liquid  or 
gas  into  a  system  in  such  a  manner  that  it  forms  a  relatively  thin  pro¬ 
tective  layer  over  the  surface  to  be  cooled.  The  liquid  or  gaa  is 
injected  at  one  or  more  discrete  locations  on  the  surface  through  holes 
or  slots.  Transpiration  cooling  refers  to  the  injection  of  a  liquid  or 
gas  through  a  porous  surfaea  to  cool  it.  Gas  film  cooling  by  Injection 
through  holts  approaches  transpiration  cooling  as  the  size  of  the  holes 
decreases  and  the  number  increases,  however,  in  actual  practice  the  two 
are  usually  wall  dlatlngulshad.  Figure  1  shows  the  two  systems 
dlairasaatlcally. 

.  Zt  la  instructive  to  eonalder  the  desirable  properties  for  film 
and  transpiration  coolants.  In  order  to  maintain  a  low  surface 
temperature  with  liquids  they  should  have  s  low  saturation  temperature 
et  the  static  pressure  in  the  chamber.  They  should  also  have  e  high 
heat  of  vaporization  end  e  high  specific  heat  eepeeity  in  the  vapor 
state.  This  reduce*  the  emount  of  coolant  required.  The  primary  re¬ 
quirement  for  e  gaseous  coolant,  from  a  thermal  standpoint,  is  that  it 
have  a  high  hast  capacity,  This  lndicstss  that  Tight  gases,  such  as  heli¬ 
um  and  hydrogen  should  make  good  coolant'  fr  „*  flits  end  transpiration 
cooling.  There  are  other  requirements  to  be  considered,  however.  The 
introduction  of  the  coolant  into  the  system  must  not  produce  effects 
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Coolant  Liquid 


Liquid  Film  Cooling 


Transpiration  Cooling 


Fig.  1  Diagram  of  Liquid  Film  Cooling  and  Transpiration  Cooling 
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which  destroy  Che  usefulness  of  Che  system,  This  would  not  severely 
limlc  Che  choice  of  coolanCs  for  a  rocket  nozzle  or  a  re-entry  vehicle, 
but  it  muse  be  carefully  considered  for  a  wind  tunnel.  If  possible 
wind  tunnels  should  CesC  vehicles  In  the  atmosphere  in  which  they  are 
to  operate.  The  main  gas  stream  thus  will  normally  be  air.  Nitrogen 
might  be  used  in  some  cases.  If  the  cooling  process  requires  in¬ 
jection  of  an  amount  which  is  a  significant  fraction  of  the  main  gas 
mass  flow  rate  of  the  wind  tunnel  the  required  coolant  would  thus  have 
to  be  either  air  or  nitrogen. 

In  both  film  and  transpiration  cooling  the  protective  effect  of 
the  heat  absorption  by  the  coolant  la  augmented  by  the  heat  transfer 
blocking  action  of  the  coolant  gas  as  it  moves  into  the  main  stream. 

This  "blowing"  action  at  the  wall  or  coolant  film  surface  reduces  the 
velocity  and  temperature  gradients  at  the  surface  and  thus  decreases  the 
heat  transfer  rate. 

Backside  cooling  may  be  used  with  film  cooling.  This  combi¬ 
nation  reduces  the  nozzle  liner  temperature  and  allows  the  cooling  of  a 
greater  nozzle  wall  area  by  the  same  amount  of  film  coolant.  It  is  more 
effective  for  liquid  films  under  the  following  conditions: 

1,  The  film  saturation  temperature  is  considerably  higher 
than  the  backside  coolant  temperature. 

2.  There  is  a  high  thermal  conductance  between  the  coolant 
film  and  the  backside  coolant. 
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II.  REVIEW  OF  PREVIOUS  STUDIES 

Analytical  and  experimental  studies  of  nozzle  heat  transfer  rates 
have  been  made  by  a  number  of  investigators;  however,  studies  of  heat 
transfer  with  extremely  high  heat  loads  using  nozzles  instrumented  to 
give  a  good  profile  of  the  heat  transfer  rate  are  not  available.  Zn  the 
present  study  such  a  nozzle  was  used. 

Analytical  and  experimental  studies  of  film  and  transpiration 
cooling  have  also  been  made  by  many  Investigators,  The  Bibliography 
contains  references  to  several  of  these  studies.  An  extensive  review  of 
the  literature  on  film  cooling  studies  up  to  1957  la  given  in  Graham  and 
gucrov's  "Film  Cooling,  Its  Theory  and  Application,"  (16).  Llbrlzzl  and 
Cresci’s  "Transpiration  Cooling  of  a  Turbulent  Boundary  Layer  In  an 
Axlaymetric  Nozzle"  (24)  lists  a  large  amount  of  work  which  has  been 
done  on  film  and  transpiration  cooling  since  1957.  A  careful  exami- 
nation  of  these  references,  however,  disclosed  that  none  of  these 
analyses  is  directly  applicable  to  the  present  problem. 

Since  the  present  solution  to  these  problems  Is  based  on  adap¬ 
tation  of  previous  works  (except  for  gas  film  cooling),  a  somewhat  more 
extensive  review  of  these  works  is  In  order. 

A  widely  used  equation  for  rapid  calculation  of  heat  transfer 
coefficients  in  nozzle  throats  without  film  or  transpiration  cooling 
was  developed  by  Bartz  (2),  Re  started  from  the  basic  premise,  based 
on  more  elaborate  boundary  layer  calculations  performed  earlier  (3) 
that  for  a  turbulent  boundary  layer  the  dominant  variable  factor  is 
the  mass  flow  rate,  l.e. 
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h  -V  (pu)n  (1) 

He  Chen  rewrote  Eq.  (1)  In  the  form 

Nu  -  C(Re)m  (Pr)n  (2) 

which  is  the  same  form  as  the  Widely  used  Dittus-Boelter  equation  for 
correlation  of  fully  developed  turbulent  heat  transfer  In  pipe  flow. 

By  regrouping  the  variables  in  Eq.  (2),  assuming  y*\.  T6* ,  p  \  Cp  30(1 

Pr  constant,  and  Introducing  a  factor  (D^/r^)  to  account  for  curvature 
effects  at  the  throat,  Bartz  was  able  to  write  Eq.  (2)  In  the  form 


The  constant  C  In  the  equation  was  determined  from  experimental  data 
for  heat  transfer  at  nozzle  throats  to  be  0,0265  In  a  nozzle  with  con¬ 
traction  and  expansion  half  angles  of  30  and  15  degrees,  respectively. 
This  is  only  slightly  higher  than  that  usually  given  for  the  constant  C 
when  Eq.  (3)  is  applied  to  fully  developed  turbulent  pipe  flow, 
namely  0.023.  An  examination  of  the  velocity  profiles  for  turbulent 
flow  with  negative  axial  pressure  gradients  (39)  shows  that  the  pressure 
gradient  has  the  effect  of  increasing  the  wall  shear  stress,  and  thus 
an  Increase  in  the  heat  transfer  coefficient  is  to  be  expected. 
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Emmons  (13)  developed  a  simplified  analysis  for  liquid  film 
cooling  of  a  cylindrical  rocket  motor  combustion  chamber  with  backside 
cooling.  It  was  based  on  the  following  assumptions: 

1*  The  gas  stream  pressure  and  velocity  do  not  vary 
appreciably  in  the  axial  direction, 

2,  The  wall  temperature  does  not  vary  appreciably  in  the 
axial  direction. 

3,  The  temperature  of  the  liquid  film  surface  exposed  to 
the  hot  gas  is  at  all  points  equal  to  the  saturation 
temperature  of  the  liquid, 

4,  Temperature  gradients  in  the  radial  direction  through 
the  liquid  film  are  small. 

5,  The  liquid  film  temperature  does  not  vary  rapidly  in 
the  axial  direction. 

6,  The  heat  transfer  coefficient  between  the  gas  stream 
and  the  liquid  film  is  constant  and  is.  based  on  a 
steady  state  analysis. 

7,  The  amount  of  liquid  in  the  film  varies  linearly  with 
the  distance  in  the  axial  direction. 

In  the  case  studied  by  Emmons  these  were  all  reasonable  as¬ 
sumptions  provided  the  coolant  film  was  Injected  at  the  local  saturation 
temperature,,  It  may  be  noted  immediately,  however,  that  several  of 
these  assumptions  are  not  reasonable  for  the  case  of  injection  of  a 
coolant  liquid  well  below  saturation  temperature  into  a  high  pressure, 
high  temperature  nozzle  gas  stream.  In  a  supersonic  nozzle  there  are 
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large  axial  variations  in  static  pressure,  film  temperature,  heat 
transfer  coefficient,  film  thickness,'  and  film  shear  s tress «  The 
rate  of  evaporation  is  not  constant  and  thus  the  film  mass  flow  rate 
does  not  vary  linearly  with  axial  distance.  To  accurately  predict 
the  required  film  cooling  for  a  nozzle  using  Emmons'  approach  it 
was  thus  necessary  to  remove  the  restrictions  of  the  above  assump¬ 
tions. 

III.  PLAN  OF  THE  PRESENT  INVESTIGATIONS 

Objectives.  The  objectives  of  the  present  Investigation  were 
to  perform  analytical  and  experimental  studies  of  heat  transfer  rates 
and  required  coolant  injection  rates  for  a  hypersonic  nozzle  with 
gas  film  cooling  and  to  make  analytical  studies  of  heat  transfer  r*.«s 
and  required  coolant  flow  rates  for  a  hypersonic  nozzle  with  liquid 
film  and  transpiration  cooling. 

Theoretical  Investigations.  For  liquid  film  cooling  and  tran¬ 
spiration  cooling  the  results  of  previous  investigations,  particularly 
those  outlined  above,  were  examined,  suitably  modified  and  combined  in 
a  finite  difference  stepwise  calculation  through  the  nozzle.  This  was 
programmed  for  solution  using  an  IBM  7040  digital  computer. 

For  gas  film  cooling  a  method  of  solution  was  programmed  which 
consisted  of  determining  the  diffusivltles  of  momentum  and  heat  based 
on  a  two-region  universal  pipe  velocity  profile  for  turbulent  flow; 
then  performing  layerwlsa  momentum  and  energy  balances  to  determine 
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th«  temperature  profile  and  heat  transfer  rate  at  the  wall.  This 
program  was  also  run  on  the  IBM  7040  digital  computer  and  the  reaulta 
correlated  with  the  experimental  results. 

Experimental  Investigations.  A  temperature  instrumented 
nozzle  was  constructed  for  the  experimental  part  of  this  study.  It 
was  instrumented  to  measure  the  outside  nozzle  wall  temperatures  at 
one-eighth  inch  Intervals  along  the  nozzle.  These  temperature 
readings,  along  with  known  flow  conditions  in  the  backside  cooling 
annulus,  gave  sufficient  information  to  calculate  local  heat  transfer 
rates.  The  nozzle  was  also  constructed  with  an  air  film  injector 
which  Injected  a  controlled  film  of  coolant  air  along  the  convergent 
section  of  the  noz2le.  Using  the  experimental  information,  a 
comparison  could  be  made  to  the  theoretical  heat  fluxes  and 
temperature  profiles  for  various  nozzle  flow  conditions. 
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CHAPTER  II 
GAS  FILM  COOLING 

I.  GENERAL- DISCUSSION 

A  large  number  of  Investigators  have  workea  on  solutions 
of  the  equations  of  continuity  of  mass,  momentum,  and  energy.  Bird, 
Stewart  and  Llghtfoot  (7),  Bennett  and  Myers  (5),  and  other  authors 
Include  In  their  books  many  of  the  available  solutions.  Other  _ 
solutions  are  available  In  the  journals  concerning  fluid  flow  and 
heat  transfer.  A  careful  study  shows  that  closed  type  analytical 
solutions  have  been  accomplished  for  only  a  relatively  few  rather 
special  cases.  With  the  development  of  high  speed  computing  equip¬ 
ment  numerical  methods  of  solution  for  the  more  difficult  cases 
ere  now  feasible.  The  present  investigator  decided  to  make  use  of 
such  capability  In  the  study  of  film  and  transpiration  coding. 

For  gas  film  cooling  a  simplified  model  was  adopted  which 
disregards  diffusion  or  mixing  of  the  coolant  gas  with  the  main¬ 
stream  gas.  Since  Pr  \  1  for  most  gases  it  was  necessary  to  solve 
for  the  developing  velocity  and  temperature  profiles  simultaneously. 
This  was  done  by  making  certain  assumptions  about  the  viscosity 
and  thermal  conductivity  for  Individual  gas  layers  and  then  making 
layerwlse  momentum  and  energy  balances  over  each  axial  increment. 
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Libby  and  Pallone  (23)  worked  out  and  published  in  1954 
a  rather  extensive  Integral  method  for  handling  the  calculation 
of  velocity  and  enthalpy  profiles  for  a  gas  film  introduced  into 
a  laminar  compressible  boundary  layer.  Sixth  degree  poly¬ 
nomials  were  used  to  describe  the  velocity  and  enthalpy  profiles. 
Since  the  species  concentration  equation  was  not  used  the  injected 
fluid  was  assumed  to  be  the  same  as  the  mainstream  fluid.  The 
polynomial  coefficients  were  chosen  without  particular  regard  for 
the  asymptotic  (far  downstream)  results  and  thus  the  solution  is 
not  asymptotically  correct.  The  numerical  results  do  indicate 
a  Btrong  effect  of  the  Injected  fluid  for  some  distance  down¬ 
stream  from  the  point  where  injection  ceases. 

Baron  (1)  worked  out  and  published  in  1956  an  extension 
of  the  Libby  and  Pallone  method.  He  included  the  species 
concentration  equation  for  a  binary  fluid  boundary  layer  and 
adjusted  the  form  of  the  polynomials  bo  as  to  assure  correct 
asymptotic  solutions. 

Neither  of  these  solutions  Is  applicable  to  the  present 
case  of  nozzle  throat  cooling  with  Injection  into  a  turbulent 
boundary  layer. 
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Emmons  (13)  and  Hatch  and  Papell  (17)  used  a  very  crude 
flaw  model  for  tangential  coolant  gas  Injection  Into  turbulent 
pipe  flow.  The  coolant  gas  was  assumed  to  form  a  layer  at 
Che  wall  with  no  diffusion  or  mixing  with  the  mainstream.  A 
heat  transfer  coefficient  h  was  used  to  calculate  the  heat  trans¬ 
fer  rate  to  the  coolant  layer  from  the  mainstream  and  from  the 
coolant  layer  to  the  nozzle  wall.  The  heat  transfer  coefficient 
was  calculated  from  the  mainstream  conditions  as  though  the 
coolant  layer  were  not  present.  It  may  be  noted  immediately 
that  this  model,  besides  its  other  deficiencies,  is  not  asymptoti¬ 
cally  correct.  It  predicts  a  downstream  heat  transfer  rate 
which  la  only  one  half  the  correct  rate  for  the  far  downstream 
value. 


II.  GAS  FILM  THEORETICAL  ANALYSIS 

The  approach  taken  in  the  present  investigation  was  to  solve 
the  continuity,  momentum,  and  energy  balance  equations  for  the 
boundary  layer  written  in  finite  difference  form,  for  turbulent  main¬ 
stream  flow  with  tangential  Injection  of  a  turbulent  coolant  gas> 
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The  necessary  assumptions  are: 

1.  The  coolant  Injection  is  a  low  percentage  of  the  main¬ 
stream  flow  rate  and  has  little  effect  on  the  expansion 
process  of  the  mainstream, 

2.  The  injected  film  is  turbulent  and  a  laminar  sublayer  is 
immediately  established. 

3.  The  laminar  sublayer  thickness  and  the  eddy  dlffusivitles 
of  momentum  and  heat  may  be  calculated  from  the  two  region 

✓* 

universal  turbulent  velocity  profile  equations. 

4.  The  ratio  of  the  eddy  dlffusivitles  of  momentum  and  heat 
is  a  constant.  As  a  preliminary  assumption  this  constant 
is  taken  to  be  one. 

5.  The  vail  shear  stress  may  initially  ba  calculated  using 
the  Blaslus  turbulent  -  smooth  pipe  -  flow  formula. 

Based  on  the  above  assumptions  the  coolant  and  mainstream 
gases-  near  the  vail  may  be  divided  into  finite  layers  parallel  to  the 
wall  and  finite  Increments  axially  along  the  surface.  The  velocity 
and  temperature  profiles  are  then  calculated  in  a  stepwise  fashion 
through  the  nozzle.  The  coordinate  system  is  thus  one  which  moves 
with  the  fluid. 

As  was  stated  above,  the  analysis  depends  upon  the  solution 
of  the  continuity,  momentum,  and  energy  aquations.  The  continuity 
equation  is  satisfied  by  assigning  a  fixed  amount  of  fluid  to  each 
layer  except  for  the  laminar  sublayer  and  the  neighboring  turbulent 
layer.  The  exchange  between  these  layers  is  calculated  et  each 
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increment  end  thus  continuity  assured. 

A  steady  state  momentum  balance  for  the  coolant  layer  at 
the  coolant-mainstream  interface  aa  shown  in  Figure  2,  is 
written  as  * 


Pressure  forces 

Inertia  forces 

Shear  forces  on 

in  A  direction 
on  coolant  layer  1 

+ 

in  A  direction 
on  coolant  layer  1 

+ 

top  and  bottom  of 
coolant  layer  1 

Using  the  notation  shown  in  Figure  2*  noting  that  primed  quantities 
refer  to  the  coolant  stream  and  subscripts  refer  to  layer  position# » 
the  resulting  mathematical  momentum  equation  is 


+ 


J  rdf  (p!  -  P2>1  + 

wf  v'20 

l—  -J 

8c 

*df  A  (vlj  +  V21  -  y'z1  - 

*1  +  *f 

Piei  P ! e  f 

Jrdf  A  (v'l  +  y'22  ~  V'l,  -  v'2,) 
2  £  1 

«1  +  S 2 
Tief  P2e2 


0 


(6) 


The  differential  equation  and  development  is  shown  in  Appendix  6. 
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Solving  Eq.  (6)  for  the  coolant  velocity  at  the  outlet  of  the  first - 
layer,  gives 


wi  v'l  +  A{vi  +  V2.  -  v'l  )  +  B(v'i  +  v"2  -  v'i  )  +  P 

x  1  111  2^1 


^  +  A  +  B 
«c 


where 


<$1 

+  6i 

Plel 

PieJ 

irdfA 

■jL 

+  _ 

*>ut 

P$e£ 

(8) 


(9) 


P  =  Tr(d{)  (<5J)(pi  -  P2> 


(10) 


A  EDC-TR-66-86 


Nozzle  Center  Line 


Fig.  2 


Velocity  and  Temperature  Notation  for  Gas  Film  by  Layers  and  Increments 
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A  set  of  such  equations  is  written  for  the  fluid  layers  in  the 
coolant  from  the  coolant-mainstream  interface  to  the  nozzle  vail. 
Another  set  is  written  for  the  mainstream  fluid  layers  from  the 
coolant-mainstream  interface  to  a  point  in  the  mainstream  where 
there  is  negligible  variation  in  exit  velocity  between  layers.  The 
wall  boundary  condition  requires  zero  velocity  at  the  wall.  The 
mainstream  boundary  condition  may  be  handled  either  by  using  enough 
layers  of  fluid  so  they  contain  several  times  the  mass  flow  rate  of 
the  coolant  or  by  specifying  a  smaller  number  of  layers  and  adding 
layers  if  the  chosen  number  does  not  cover  the  velocity  boundary 
layer.  The  outer  mainstream  fluid  layer  is  then  fixed  at  the  main¬ 
stream  velocity  calculated  from  frozen  flow  conditions  established 
by  an  lsentroplc  expansion  with  a  constant  ratio  of  specific  heats. 
The  Gauss-Siedel  method  of  successive  substitution  is  used  to  solve 
the  equations  iteratively. 

In  order  to  make  the  ebove  calculations,  the  dlffusivlty  of 
momentum  must  be  evaluated  at  each  layer.  In  the  laminar  sublayer, 

r 

0<y  <10,  the  molecular  dlffusivlty  of  momentum  is  used.  In  the 
turbulent  core  assumption  3  is  used  and  the  eddy  diffuslvities 
calculated  from  the  universal  velocity  profile  equation  for  this 
region,  which  ib: 


u 

* 

u 


* 

2.5  In  +  5.5 
u 


Differentiating  Eq.  (11)  with  respect  to  y  gives : 


£-2.5*- 
dy  y 


(11) 


(12) 
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Now  for  the  turbulent  core 


TSC  .  du  / 

— £  -  e  —  *  (u  ) 
P  dy 


(13) 


thus 


£  - 


p 

~*3iT 

dy 


(14) 


Assuming  the  velocity  boundary  layer  terminates  near  the  wall 


T  “  T 


W 


(15) 


so 


e  ■ 


2.5 


(16) 


or 


e  -  0.4  u  (r^  -  r) 


(17) 


In  the  notation  of  Figure  2,  page  16,  Eq.  (17)  becomes 


*  4  l 

-  0.4  u  (r, - =) 

1  L  O 


(18) 


or 


*  ^1 
:t  -  0.4  u  (zi  -  f>) 


(19) 
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In  order  to  evaluate  the  diffusivities  using  either  Eq.  (18) 

or  Eq.  (19),  u  and  thus  is  required  at  each  increment.  Assumption  5 

is  used  to  evaluate  t  for  the  initial  increment,  i.e. 

w  ’ 


T 


W 


.053  (Re) 


2 

pu 


(20) 


Following  the  initial  increment,  may  be  evaluated  at  each  suc¬ 
cessive  increment  either  from  Eq.  (20)  or  from  the  wall  velocity 
gradient  at  the  exit  of  the  previous  increment. 

Energy  balance  equations  are  written  for  the  same  fluid 
layers  in  order  to  calculate  temperatures  at  the  exit  of  the 
increments  -  (Refer  to  Figure  2  for  notation) .  The  steady  state 
energy  equation  for  the  coolant  layer  at  the  coolant-mainstream  inter¬ 
face  is  written  as 


+ 


+ 


Enthalpy  change 
in  A  direction  of 
coolant  layer  1 


Energy  transfer  from 
gas  layer  1  to  coolant 
layer  1  by  diffusion 

Work  at  boundaries 
on  coolant  layer  1. 
Shear  stress  effects 


+ 


+ 


0 


Kinetic  energy  change 
in  A  direction  of 
coolant  layer  1 

Energy  transfer  from 
coolant  layer  2  to 
coolant  layer  1  by  diffusion 


(21) 
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The  mathematical  statement  of  Eq.  (21)  is 


1 —  v  ,*71 

■ 

2  2 

wlcl (T/l1  -  T  2X) 

+ 

wf(v'*l1  -'V^j) 

L  _ 

2SCJ 

-  T/l]_  -  T/21) 

si  «r 

—  ■¥  — 

PlClEHJ)  PlCfGflll 

irdiA(T'l2  +  T'2g  -  T/l1  -  T'2^ 

fit  S  2 

P  lcl  e  14*  +  P2C2*2* 


TTdfdSiPitiCvli  +  V2i  -  v'l  -  v'2  ) 


Jg„(26J  +  d|  +  6l) 


-  0 


(22) 


Solving  Eq.  (22)  for  T/2^,  the  exit  temperature  of  the  first  coolant 
layer  gives: 


T/l1(w'*c'‘  -  A  -  B)  +  +  T21)  +  Bte'lg  +  T%)  E  t  F 

v'c'  -  A  +  B 


(23) 


where 

A  - 


irdj'A 

"fil  '  +  "fit 

Piciei^  Ptcfet^ 


(24) 
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TTdi'A 


P  l’ci'e  i'll*  P2c£c$^ 


(25) 


WlCv'lj  -  v'2  ) 

2JT  J_ 


(26) 


irdiASiPie^vLi  +  v2j  -  v'lg  -  v'2g) 

■  - —  1  — 1  2 - 

Jg.(28f  +  fi2  +  dj) 


(27) 


Since  the  energy  equations  contain  the  entrance  and  exit 
velocities,  the  velocity  solution  is  performed  first  using  layer 
properties  based  an  the  exit  temperatures  from  the  preceding  incre¬ 
ment.  The  velocity  equations  are  dependent  on ■ temperature  only  , 
through  the  variation  of  the  fluid  properties  with  temperatures. 

If  desired  the  fluid  properties  could  be  re-evaluated  after  com¬ 
pletion  of  the  temperature  iteration. 

When  the  velocity  iteration  is  completed  the  temperature 
iteration  is  performed  in  like  manner  by  successive  substitution. 

The  analysis  described  above  has  been  programmed  in. 

Fortran  IV  language  for  execution  on  a  digital  computer.  The 

program  is  included  in  Appendix  C. 

Included  in  the  program  are  provisions  for- wall  shear  stress 
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calculation  either  from  the  Blasius  Equation  or  from  the  past 
Increment,  multiple  injection  calculations,  and  provision  for 
approaching  a  transpiration  calculation  by  injecting  a  very  small 
amount  of  coolant  at  each  Increment.  New  coolant  layers  are  added 
for  the  multiple  injection  calculation.  The  fluid  is  added  at 
a  specified  increment,  or  increments  and  the  wall  side  eoolant 
layer  is  recalculated  and  layers  added  as  needed.  Provision  is 
made  for  backside  cooling  in  the  case  of  single  and  multiple 
injection,  while  an  adiabatic  wall  may  be  specified  for  any  type 
injection. 

111.  EXPERIMENTAL  INVESTIGATION 

Design  Requirements.  An  examination  of  some  published 
data  for  film  cooled  nozzles  (6),(19),  (24),  (25),  indicated  that 
in  most  cases  either  insufficient  data  are  available  to  make  a 
detailed  comparison  with  the  present  analysis  or  the  fluid 
stagnation  conditions  are  not  in  the  range  of  present  interests. 
Therefore  an  instrumented  nozzle  was  designed  and  fabricated  for 
this  investigation.  The  general  geometrical  requirements  for  this 
nozzle  were  that  it  have  an  entrance  diameter  of  1.5  inches,  a 
throat  diameter  of  0.5  Inches,  and  an  expansion  angle  of  15°.  Ad¬ 
ditional  requirements  were  that  it  be  mechanically  designed  to  handle 
air  as  the  expanding  fluid  with  stagnation  pressure  up  to  ^  50  atmos¬ 
pheres  and  stagnation  temperatures  up  lo  u  Id, 000  °R.  In  order  to  make 
the  required  correlation  with  the  analytical  study,  instrumentation  for 
the  nozzle  was  required  to  produce  data  necessary  for  calculation  of 
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che  local  nozzle  wall  Inside  temperature  and  heat  loads  as  distributed 
axially  through  the  nozzle. 

Nozzle  Design.  Details  of  the  materials  study,  stress 
analysis  and  backside  cooling  problems  have  been  previously  published 
(27).  It  was  found  that  the  combination  of  high  thermal  stresses 
and  hoop  stresses  in  the  nozzle  could  best  be  handled  in  the  range  of 
interest  by  using  a  thin  vailed  nozzle  liner  turned  from  copper- 
zirconium  alloy.  In  higher  pressure  and  heat  load  regions  of  the 
nozzle  the  wall  thickness  was  0.040  Inches.  Downstream  from  the  throat 
the  liner  gradually  thickened  to  approximately  0.15  inches  near  the  exit. 
The  liner  was  surrounded  by  a  steel  split  sleeve  turned  to  provide  a 
backside  coolant  flow  channel  of  constant  flow  area  around  the  liner. 

The  liner  and  sleeve  were  surrounded  by  a  steel  cylindrical  outer  casing. 
The  outer  casing  was  provided  with  passages  for  thermocouple  lead  access 
to  the  nozzle  liner  through  the  split  sleeve. 

Twenty-four  duplex  copper-constantan  thermocouples  were  led 
in  on  each  side  of  the  liner,  passing  through  a  narrow  slot  milled  in 
the  split  sleeve,  through  small  holes  drilled  in  a  positioning  rib,  and 
attached  to  the  nozzle  liner  with  600  °F,  softening  temperature  solder. 

The  thermocouples  started  0.875  inches  upstream  from  the  point  of  smallest 
radius  in  the  nozzle  and  were  spaced  0.125  inches  apart  on  the  liner 
surface.  0-ring  seals  were  provided  on  either  side  of  the  thermocouple 
passage. 

A  full-scale  sectional  view  of  the  liner,  sleeve,  and  casing 
shoving  the  thermocouple  installation,  is  given  in  Figure  3 . 
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Film  Inlector  Design.  An  Injector  plate  wee  designed  to 
natch  the  nozzle  liner  and  casing.  It  consisted  of  a  single  piece 
turned  from  cop per- zirconium  and  drilled  to  provide  18  coolant  flow 
passages  evenly  spaced  about  ehe  circumference.  These  coolant  passages 
discharged  Into  a  snail  plenum  chamber  Just  before  the  film  entrance 
annulus.  Provision  was  made  for  varying  the  film  injection  slot 
width  by  placing  spacer  shims  between  the  injector  plate  and  the 
nozzle  casing.  Figure  4  shows  a  sectional  view  of  the  Injector  plate. 

Stilling  Chamber  Design.  Experience  with  the  previous 
instrumented  nozzle  (27)  showed  that  the  nozzle  wall  heat  load 
was  somewhat  lower  than  that  anticipated.  There  was  also  a 
swirl  component  of  velocity  of  the  mainstream  gas.  A  possible 
reason  for  the  low  heat  flux  In  the  first  experiments  was  a 
"natural**  film  cooling  effect  caused  by  the  peaking  of  the 
temperature  at  the  centerline  of  the  arc  heater  tube,  Thla,  together 
with  the  centrifugal  effect  of  the  swirl  flow,  would  tend  to 
produce  lower  temperatures  near  the  wall  and  thus  lower  heat  fluxes. 

In  order  to  remove  the  mainstream  swirl  and  provide  a  more  uniform 
temperature  profile  at  the  nozzle  entrance,  a  stilling  chamber  was 
designed,  fabricated,  and  installed  between  the  arc  air  heater  and  the 
film  Injector  station.  This  chamber  consisted  of  a  brass  liner 
0.125  Inches  thick  and  a  steel  outer  casing.  It  was  cylindrical  In 
shape  with  an  inside  diameter  of  ^5  Inches  and  an  Inside  length  of 
^,8  inches*  The  liner  consisted  of  a  cylindrical  sleeve  and  two 
curved  end  pieces.  Cooling  water  paasages  were  provided  between 
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Fig«  4  Injection  Plot® 
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Che  casing  and  Che  liner  and  a  pressure  cap  was  provided  at  the  downstream 
end.  Figure  5  shows  a  half-scale  sectional  view  of  the  stilling  chamber,, 

Arc  Air  Heater.  The  air  for  expansion  through  the  nozzle 
was  heated  by  passage  through  an  electric  arc  heater  supplied  hy  the 
Linde  Company,  Indianapolis,  Indiana-  The  heater  has  a  capacity  of 
1-3  pounds  per  second  of  air  flow  at  pressures  up  to  ~1400  psia  with 
a  3 /8  inch  diameter  constrictor  nozzle .  Energy  can  be  supplied  to  the 
air  up  to  a  maximum  rate  of  ~4  megawatts-  The  air  heater,  stilling 
chamber,  injection  plate,  and  nozzle  are  shown  assembled  in  Figure  6. 

Instrumentation,  The  experiment  was  set  up  to  be  operated 
remotely  due  to  the  high  pressures  and  temperatures  Involved.  Instrumentation 
for  the  operation  of  the  air  supply  and  air  heater  was  designed  by  the 
staff  at  Arnold  Engineering  Development  Center  and  will  not  be  further 
discussed  In  the  present  report.  Instrumentation  for  the  stilling 
chamberg  injection  station  and  nozzle  consisted  of  temperature  and 
flow  measuring  equipment  for  water  and  atr  as  supplied  to  the 
particular  pieces  of  equipment. 

Water  flow  rate  to  the  stilling  chamber  was  measured  by  a 
calibrated  flow  rate  meter0  The  inlet  to  exit  temperature  difference 
of  the  stilling  chamber  cooling  water  was  measured  by  thermistors „  The 
pressure  at  the  stilling  chamber  exit  was  measured  by  a  pressure 
transducer. 

Air  mass  flow  rate  to  the  coolant  injector  was  established  by 

choked  flow  through  a  venturi .  The  inlet  temperature  was  measured  by 
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Fig-  5  Stilling  Chamber  Assembly 
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Fig.  6  Photograph  of  Experimental  Assembly 
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a  thermocouple. 

Water  flow  rate  to  the  nozzle  backside  coolant  channel  was 
measured  by  a  calibrated  flow  rate  meter.  The  inlet-to-exit  temperature 
difference  was  measured  by  thermistors.  The  nozzle  wall  thermocouple 
leads  were  attached  to  an  automatic  scanning  Austin  system. 

Experimental  Procedure.  The  experimental  nozzle  was  fabri¬ 
cated  and  assembled  at  The  University  of  Tennessee,  Knoxville,  Tennessee, 
under  the  supervision  of  the  investigator.  It  was  then  taken  to  the 
Arnold  Engineering  Development  Center,  Tullahoma,  Tennessee,  where  the 
actual  operation  of  the  experiments  was  performed  by  a  group  now  called 
the  experimenter. 

The  Investigator  specified  gas  film  cooling  ratios  of  approxi¬ 
mately  20Z,  10X,  and  5Z  with  mass  flow  rates  ranging  from-0.4  lbm./sec. 
to  0.6  lbm./sec.  and  three  enthalpy  values  ranging  up  to  the  maximum 
value  available  from  the  heater  at  the  given  flow  rate.  In  addition 
to  the  cases  at  the  above  film  cooling  ratios  the  nozzle  was  operated 
for  the  first  run  with  a  gas  film  cooling  ratio  of  27.4Z. 

The  gas  film  Injection  slot  was  adjusted  for  each  film  cooling 
ratio.  In  the  range  of  pressures  and  temperatures  of  the  test  cases  the 
mainstream  Inlet  velocity  wad  from  100  to  250  ft. /sec.  To  minimize  mix¬ 
ing  of  the  film  with  the  mainstream  gas  at  injection  the  slot  was  ad¬ 
justed  to  give  film  Inlet  velocities  of  approximately  150  ft. /sec.*  For 
air  injected  at  approximately  room  temperature  this  required  a  slot  width 

This  film  velocity  approximately  matched  the  mainstream  velocity  at  the 
location  of  injection. 
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of  0.002  inches  for  each  52  of  coolant.  Spacing  shims  0.004  Inches  in  thick¬ 
ness  were  provided  by  the  investigator  and  one  shim  was  inserted  for  each  52 
coolant  ratio.  Due  to  the  30°  convergent  angle  of  the  nozzle  each  shim 
increased  the  slot  width  0.002  inches , 

When  a  test  case  was  run  an  ambient  scan  of  all  the  thermo¬ 
couples  was  first  recorded.  The  heater  was  then  fired  and  the  thermo¬ 
couples  were  scanned  approximately  20  times  during  the  run.  Since 
steady  state  for  this  system  was  reached  very  quickly  the  thermocouple 
readings  during  the  data  taking  period  were  constant  to  within  approxi¬ 
mately  5  °F,  except  for  a  few  runs  when  the  heater  was  slightly  un¬ 
stable  over  a  few  scanning  periods. 

The  stagnation  enthalpy  of  the  gas  was  calculated  by  the 
experimenter  using  an  energy  balance  for  the  system  including  the  air 
heater  and  the  stilling  chamber.  This  balance  Included  the  input 
electrical  energy  to  the  heater,  the  energy  removed  in  cooling  the  heater 
and  the  stilling  chamber,  and  the  change  of  energy  of  the  mainstream  gas. 

The  raw  Instrument  readings  recorded  on  tape  were  reduced  to 
pressure,  temperature,  and  flow  rate  values  by  the  experimenter  using 
standard  programs  and  then  furnished  to  the  present  investigator. 

IV.  DISCUSSION  OF  RESULTS 

The  described  analysis  was  programmed  in  Fortran  IV  language 
for  execution  on  a  digital  computer.  The  program  is  included  in  Appendix 

-J 

C.  Experiments  with  the  gas  film  cooled,  instrumented  nozzle  were  run 
at  film  cooling  ratios  of  approximately  302,  202,  102,  and  52  of  the 
mainstream  flow  rate.  Table  X  summarizes  the  experimental  conditions 
for  all  the  cases  run. 
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Experimental  Results.  From  the  twenty-eight  cases  run  five 
were  selected  for  analytical  correlation.  They  were  chosen  to  Include 
cases  from  each  of  the  film  cooling  ratios  given  above,  and  in  addition 
two  pairs  ware  chosen  because  their  stagnation  conditions  were  close 
to  the  same  while  they  had  different  film  cooling  ratios.  From  Table  I 
the  cases  chosen  ware  Hl-1,  H2-2,  H6-2,  H5-2,  and  H9-2,  Cases  Hl-1 
and  H5-2,  while  differing  some  in  stagnation  pressure  and  temperature 
represent  approximately  the  some  severity  of  heat  loading  on  the  nozzle. 
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AEDC* 

Case 

Number 

Mainstream 

Stagnation 

Pressure 

(psia 

Mainstream 

Stagnation 

Temperature 

CKo) 

Mainstream 
Mass  Flow 
Rate 

(lbm./sec.) 

Film  Coolant 
Mass  Flow 
Rate 

(lbm.  /sec.) 

Film  Coolant 
Inlet 

Temperature 

?F.) 

Backside 

Coolant 

Flow 

Rate 

(gpm.) 

Backside 

Coolant 

Inlet 

Temp. 

<°F.) 

Backside 

Coolant 

AT 

<°F.) 

Hl-1 

349 

3840 

0.394 

0.108 

69.0 

41.5 

65.0 

2.57 

H2-1 

341 

3160 

0.412 

0.080 

49.5 

32.0 

70.0 

4.64 

H2-2 

37T 

3660 

0.419 

0.082 

52.5 

32.0 

66.0 

7.34 

H2-3 

373 

3910 

0.419 

0.092 

55.0 

32.0 

70.5 

8.37 

H3-1 

409 

3290 

0.496 

0.094 

65.5 

30.6 

71.5 

5.31 

H3-2 

434 

3660 

0.504 

0.096 

65,5 

30.7 

71.5 

8,81 

H3-3 

447 

3800 

0.498 

0.095 

65.5 

30.7 

71,5 

10.13 

H4-1 

473 

3180 

0.606 

0.110 

69.0 

34.1 

72.0 

4.41 

H4-2 

492 

3320 

0.605 

0.108 

69.0 

35,7 

72.0 

6.74 

H4-3 

511 

3620 

0.605 

0.108 

69.0 

34.2 

72.0 

8.68 

H5-1 

299 

3180 

0.401 

0.039 

55.0 

30.1 

73.0 

5.75 

H5-2 

312 

3450 

0.409 

0.039 

57.5 

30.1 

73.0 

8.18 

H5-3 

324 

3930 

0.401 

0.039 

60.0 

38.4 

73.0 

7.29 

H6-1 

381 

3120 

0.517 

0.049 

60.0 

49,2 

73.0 

4.73 

H6-2 

407 

3460 

0.525 

0.049 

63.0 

49.2 

73.0 

6.14 

H6-3 

413 

3760 

0.525 

0,049 

67.0 

49.4 

73.0 

7.03 

H7-1 

432 

2920 

0.603 

0.059 

52.0 

55.7 

71.0 

3.45 

H7-2 

470 

3280 

0.603 

0.059 

52.0 

56,3 

71.0 

5.50 

H7-3 

490 

3600 

0.603 

0.059 

52,0 

55.8 

71.0 

6,40 

H8-1 

284 

3110 

0.401 

0.020 

61.0 

51.6 

78.0 

3.39 

H8-2 

301 

3440 

0.401 

0.020 

61.0 

51.6 

78.0 

4.56 

H8-3 

313 

3720 

0,403 

0.020 

61.0 

52.0 

78,0 

5,22 

00 

CO 


★Arnold  Engineering  Development  Center 


AEDC-TR-66-38 


w 


TABLE  I  (Concluded) 


AEDC* 

Case 

Humber 

Mainstream 

Stagnation 

Pressure 

(psia.) 

Mainstream 

Stagnation 

Temperature 

TK.) 

Mainstream 
Mass  Plow 
Rate 

(lbm./sec.) 

Film  Coolant 
Mass  Flow 
Rate 

(lba*/sec.) 

Film  Coolant 
Inlet 

Temperature 

(V) 

Backside 

Coolant 

Flow 

Rate 

(gam.) 

Backside 

Coolant 

Inlet 

Temp. 

(°F.) 

Backside 

Coolar.t 

AT 

(°F.) 

H9-1 

372 

3130 

0.526 

0.023 

61.0 

51.7 

74.0 

4*43 

H9-2 

402 

3450 

0*538 

0.023 

61.0 

51.9 

74.0 

6.26 

H9-3 

421 

3710 

0*545 

0*024 

61*0 

52.0 

78.0 

6*81 

H10-1 

425 

2960 

0*580 

0*028 

55.0 

51*4 

78.0 

5*74 

H10-2 

437 

3330 

0*581 

0*028 

55.0 

51.6 

78.0 

5.74 

H10-3 

459 

3650 

0C581 

0*028 

55.0 

51.6 

78.0 

6.78 

* Arno Id  Engineering  Development  Center 
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They  provide  a  good  comparison  of  film  cooling  effectiveness  since 
case  Hl-1  was  run  with  27,4  Z  film  cooling  while  case  H5-2  was  run 
with  only  9.73  Z  film  cooling.  Cases  H6-2  and  H9-2  have  almost 
identical  stagnation  conditions,  while  being  cooled  with  9.34  Z 
and  4.27  Z  film  cooling  respectively.  The  original  wall  temperature 
data  for  the  five  cases  correlated  are  given  in  Tables  II,  III,  IV, 

V,  and  VI. 

Values  for  the  backside  cooling  heat  transfer  coefficient 
and  local  values  of  the  heat  flux  were  calculated  at  the  various 
thermocouple  positions  using  the  measured  wall  temperatures,  the  local 
dimensions  of  the  coolant  flow  annulus,  the  coolant  mass  flow  rate,  and 
the  physical  properties  of  the  coolant  evaluated  at  the  local  film 
temperature.  The  Dlttus-Boelter  equation  was  used  for  calculation 
of  the  backside  heat  transfer  coefficients. 

Theoretical  calculations  of  the  backside  surface  temperature 
and  wall  heat  flux  distributions  were  made  for  the  five  cases  named 
above.  These  are  plotted,  along  with  the  measured  temperatures  and 
with  heat  flux  curves  calculated  by  smoothing  the  experimental 
temperature  distributions,.  In  Figures  7  through  16. 

Comparison  of  Theoretics!  and  Experimental  Results .  Attempts 
were  made  to  perform  the  theoretical  calculations  described  above 
using  the  eddy  diffusivltles  based  on  wall  shear  stress  calculated  from 
both  the  Blaslus  equation  and  from  the  past  increment  value  of  the 
wall  shear  atresa.  In  aeveral  cases  both  calculations  were  performed 
successfully.,  however  in  other  cases  the  calculation  using  the  well 
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TABLE  II 

ORIGINAL  TEMPERATURE  DATA  FOR  GAS  FILM  CASE  HI-1 


Thermocouple 

Averaged 

Thermocouple 

Averaged 

Number 

Temperature 

Number 

Temperature 

<v> 

(  ) 

1 

88.00 

13 

193.33 

2 

122.60 

14 

159.59 

3 

15 

149.02 

4 

126.76 

16 

171.49 

5 

122.45 

17 

161.07 

6 

128.97 

18 

7 

206.41 

19 

135.75 

8 

196.88 

20 

113.46 

9 

263.24 

21 

94.55 

10 

265.37 

22 

96.11 

11 

23 

129.68 

12 

222.02 

24 

*  A  dash  Indicates  a  thermocouple  which  was  not  reading. 
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TABLE  Ml 

ORIGINAL  TEMPERATURE  DATA  FOR  GAS  FILM  CASE  H2-2 


Thermocouple 

Averaged 

Thermocouple  Averaged 

Number 

Temperature 

CF.) 

Number  Temperature 

CF.) 

1 

147.63 

13 

301,16 

2 

197.99* 

14' 

229.16 

3 

15 

197,80 

4 

207.48 

16 

152,25 

5 

196.11 

17 

220.70 

6 

203.00 

18 

7 

382.93 

19 

175.19 

8 

341.27 

20 

117.73 

9 

424,68 

21 

104.17 

10 

406.02 

22 

105,47 

11 

23 

126.51 

12 

344.21 

24 

*  A  dash  Indicates  a  thermocouple  which  was  not  reading 
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TABLE  IV 

ORIGINAL  TEMPERATURE  DATA  FOR  GAS  FILM  CASE  H6-2 


Thermocouple 

Number 

Averaged 

Temperatures 

cr.) 

Thermocouple 

Number 

Averaged 

Temperatures 

1 

136.67 

13 

301.60 

2 

190.64* 

14 

128.71 

3 

15 

236,09 

4 

202.42 

16 

114.16 

5 

189.19 

17 

117,05 

6 

196.72 

18 

7 

350.63 

19 

98.07 

8 

328.76 

20 

98.71 

9 

424.00 

21 

94.22 

10 

385.58 

22 

91.16 

11 

276.36 

23 

90.70 

12 

306.11 

24 

*  A  daah  indicates  a  thermocouple  which  was  not  reading. 
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TABLE  V 

ORIGINAL  TEMPERATURE. DATA  FOR  GAS  FILM  CASE  H5-2 


Thermocouple 

Averaged 

Thermocouple 

Averaged 

Number 

Temperature 

(6f.) 

Number 

Temperature 

CF.) 

* 

1 

13 

269 o 84 

2 

152.15 

14 

136.44 

3 

101.18 

15 

289.31 

4 

166.60 

16 

113.87 

5 

157.74 

17 

136.55 

6 

163.68 

18 

7 

260.71 

19 

96.68 

8 

240.69 

20 

104.22 

9 

323.00 

21 

97.40 

10 

316.40 

22 

99.07 

11 

242.53 

23 

96,05 

12 

274.65 

24 

*  A  dash  indicates  a  thermocouple  which  was  not  reading. 
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TABLE  VI 

ORIGINAL  TEMPERATURE  DATA  FOR  GAS  FILM  CASE  N9-2 


Thermocouple 

Number 

Averaged 

Temperature 

CF.) 

Thermocouple 

Number 

Averaged 

Temperature 

CFJ 

1 

* 

13 

238.65 

2 

201.50 

14 

91,89 

3 

15 

116.64 

4 

209.25 

16 

93.32 

5 

202,23 

17 

6 

212.55 

18 

100.51 

7 

359,61 

19 

97,30 

8 

301.13 

20 

108.55 

9 

430.54 

21 

104.45 

id 

291.20 

22 

108.93 

U 

169.00 

23 

115.15 

12 

163.20 

24 

*  A  dash  indicates  a  thermocouple  which  was  not  reading * 


40 


0  1.0  2.0  3.0 

Distance  from  Injection  -  inches. 


Fig.  7  Gas  Film  Temperature  Correlation  With  Data  for  Case  HI -1 
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Fig.  8  Gat  Film  Heat  Flax  Correlation  With  Data  lor  Case  HI-1 
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Fig.  9  Gas  Film  Temperature  Correlation  With  Data  for  Case  H2-2 
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Fig.  10  Gas  Film  Heat  Flux 
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Fig.  11  Gat  Film  Temper  a  turn  Correlation  With  Data  for  Case  H6*2 
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Fig.  12  Gas  Film  Heat  Flux  Correlation  With  Data  for  Case  H6-2 


^  Fig^  13  Gfl*  Film  Temperature  Correlation  With  Data  for  Casa  H5-2 
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Fig.  14  Gas  Film  Heat  Flux  Correlation  With  Data  for  Case  H5-2 
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Fig.  15  Gas  Film  Temperature  Correlation  With  Data  (or  Cose  H9-2 
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Fig.  16  Gas  Film  Heat  Flux  Correlation  With  Data  for  Case  H9-2 
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shear  stress  from  the  post  Increment  became  unstable  and  either  the 
iterative  procedures  for  the  velocity  profile  or  the  temperature  profile 
would* not  converge#  Reduction  of  the  incremental  distance  achieved 
convergence  for  some  of  the  cases  which  were  initially  unstable,  however 
in  cases  wleh  very  low  ratios  of  film  coolant  even  the  reduction  to 
increment  lengths  of  the  order  of  1/10,000  of  the  nozzle  length  did  not 
bring  about  stability  of  the  solutions#  Calculations  using  the  eddy 
diffusivlty  based  on  wall  shear  Btress  calculated  at  each  increment  from 
the  Blasius  equation,  on  the  other  hand,  remained  stable  for  reasonable 
incremental*  lengths  for  all  cases.  These  results  are  the  ones  chosen 
for  comparison  with  the  experimental  results.  Mainstream  and  coolant 
fluid  properties  as  functions  off  temperature  were  evaluated  from  the 
drfte  presented  graphically  In  Appendix  E. 

The  analytical  program  written  gives  as  output  the  local  wall 
heat  flux  and  the  streamalde  nozzle  surface  temperature.  The  calcu¬ 
lated  heat  flux  was  used  with  the  wall  thickness  and  thermal  con¬ 
ductivity  to  calculate  backside  wall  temperatures  for  correlation  with 
the  measured  backside  temperatures. 

In  addition  to  the  data  taken  by  the  present  investigator  some 
recent  data  were  available  taken  by  Lieu  (25)  for  air  film  cooling  of 
a  supersonic  nozzle.  It  was  decided  to  make  correlations  of  the  Lieu 
data  for  a  case  with  stagnation  pressure  315  psla.  and  stagnation 
temperature  1032  °R.  The  analytical  and  experimental  results  for 
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three  different  coolant  Injection  rates  are  presented  in  Figures 
17,  18,  and  19. 

As  a  check  on  the  changes  which  would  occur  in  the  theoreti¬ 
cally  predicted  wall  temperatures  and  heat  fluxes,  if  it  were  assumed 
that  the  gas  film  coolant  were  totally  mixed  with  the  free  stream  gas 
in  a  very  short  distance  from  the  injection  point,  two  additional 
theoretical  calculations  were  made.  Total  mixing  was  assumed  for 
Case  Hl-1,  which  had  27.4  %  gas  film  coolant,  and  theoretical  calcu¬ 
lations  were  made  for  wall  temperatures  and  heat  fluxes.  Case  H2-2, 
which  had  9.73  %  gas  film  coolant,  was  also  recalculated  assuming  total 
mixing  at  the  injection  point.  These  theoretically  calculated  curves 
are  shown,  along  with  the  curves  calculated  for  no  film  cooling  and 
for  full  film  cooling,  in  Figures  7,  8,  13,  and  14. 

Without  discussing  each  of  the  Figures  7  through  16  in  detail 
several  general  observations  may  be  made.  First,  the  theoretically 
predicted  temperatures  for  no  film  cooling  fall  well  above  the 
measured  temperatures  for  the  larger  ratios  of  gas  film  coolant,  but 
as  the  amount  of  film  coolant  is  decreased  the  theoretical  curve 
approaches  the  measured  temperatures  until  the  correlation  is  very 
good  for  Case  H9-2,  with  approximately  5  %  film  coolant.  This 
indicates  that  the  predicted  temperatures  for  no  film  cooling  using 
the  present  analysis  should  be  reasonably  accurate.  As  a  further 
check  on  the  theoretical  model,  with  no  film  cooling,  comparison  was 
made  to  data  published  by  the- Jet  ' Propuls  ion  Laboratory  {42-j  and  the 
results  are  shown  in  Figure  20.  This  figure  further  demonstrates 
the  validity  of  the  theoretical  model  to  predict  temperature  or  heat 
flux  profiles,  at  least  through  the  throat.  With  the  accuracy  of 
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the  no  film  cooling  calculations  established,  examination  of  the  above 
named  figures  indicates  that  the  gas  film  coolant  does  measurably 
lower  the  wall  temperatures  and  heat  fluxes,  particularly  in  the  con¬ 
verging  section  of  the  nozzle.  Third,  from  Figures  7,  8,  13,  and  14, 
which  contain  curves  calculated  assuming  complete  mixing  of  the  coolant 
with  the  mainstream  at  injection,  it  may  be  seen  that  such  predicted 
values  fall  well  above  the  measured  values  of  wall  temperature  and  heat 
flux. 

The  effectiveness  of  the  gas  film  for  protecting  the  nozzle  wall 
is  further  shown  in  Figure  21,  which  shows  photographs  of  the  converging 
section  of  the  nozzle  and  the  inside  surface  of  the  injection  ring. 

These  photographs  were  made  without  cleaning  up  the  surfaces  following 
the  conclusion  of  the  experimental  program.  The  conical  converging 
section  of  the  nozzle  is  still  relatively  smooth  and  clean.  There  is  a 
noticeable  change  very  close  to  the  throat,  where  the  nozzle  is  somewhat 
roughened  and  discolored.  The  inside  of  the  injection  ring  is  badly 
eroded.  In  a  previous  nozzle  operated  without  gas  film  cooling  (27) , 
the  entire  convergent  portion  was  eroded  during  the  experiments. 

Figures  22  and  23  show  a  direct  comparison  of  the  wall  heat 
flux  for  the  two  pairs  of  cases  with  very  similar  stagnation  conditions. 
In  Figure  22  the  larger  heat  flux  for  the  case  with  approximately  10  % 
cooling  as  compared  with  the  case  for  approximately  30  %  cooling  is 
evident.  In  Figure  23,  for  cases  with  approximately  5  %  and  10  % 
cooling  there  is  very  little  difference,  but  an  examination  of  the 
tabulated  experimental  temperatures  for  the  two  cases  shows  that  the 
temperatures  near  the  nozzle  entrance  for  the  5  %  cooling  case  were 
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Fig.  21  Photograph  of  Noxile  Liner  and  Injection  Ring  after 
Termination  of  Experimentation 
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Fig,  22  Comparison  of  Hast  Flexes  for  Cates  HI-1  and  H5-2 


Fig.  23  Comparison  of  Heat  Fluxes  for  Cases  H6-2  and  H9-2 
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slightly  higher.  Comparison  of  the  backside  coolant  flow  rates  for 
the  two  cases  from  Table  I,  pages  33  and  34,  shows  that  the  5  %  case 
was  more  strongly  cooled.  Taken  together  this  is  an  indication 
that  the  gas  film  coolant  was  effective  in  the  converging  section. 

Study  of  the  measured  temperatures  in  Figures  7  through  16 
shows  a  quick  rise  in  temperature  near  the  nozzle  entrance,  then 
a  rather  steady  temperature  until  it  rises  rapidly  toward  the  throat. 
The  present  investigator  believes  this  is  due  to  a  certain  amount 
of  gross  mixing  at  or  near  the  injection  point.  This  "mixed  film" 
then  forms  the  actual  film  coolant. 

The  temperature  curves  also  reveal  that  as  the  nozzle  was 
operated  for  longer  periods  of  time  the  differential  thermal  expansion 
between  the  liner  and  the  instrument  rib  caused  the  far  downstream 
thermocouples  to  loosen  from  the  wall  and  read  low.  This  was  caused 
by  the  freeing  of  the  downstream  end  of  the-  liner  for  expansion. 

Notice  in  this  connection  thermocouples  14  through  24  in  Figures  11, 
13,  and  15.  This  also  explains  the  poor  correlation  between  the 
curves  in  Figure  22  for  the  section  downstream  from  the  throat. 

Comparison  of  the  calculated  temperatures  with  the  measured 
values  in  the  case  of  the  Lieu  nozzle,  Figures  17,  18,  and  19  on 
pages  52,  53,  and  54,  respectively,  show  reasonably  good  correlation 
for  the  4.35  %  injection  ratio  but  poorer  correlation  for  the 
higher  injection  rates.  The  present  investigator  noted  some  apparent 
discrepancy  between  Lieu's  data  and  the  present  data,  however.  For 
Lieu's  nozzle  the  maximum  temperature  showed  a  definite  shift  upstream 
as  the  film  cooling  was  increased.  Both  the  analytical  prediction  and 
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the  experimental  data  in  the  present  investigation  show  maximum 
temperatures  immediately  downstream  from  the  throat  with  very  little 
shift  in  position  as  the  injection  ratio  varies  from  ^5%  to  ^30%. 
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CHAPTER  III 
LIQUID  FILM  COOLING 

I.  THEORETICAL  ANALYSIS 


As  was  stated  in  Chapter  I,  Bartz's  equation  for  the  heat 
transfer  coefficient  in  a  nozzle  without  coolant  injection  is  given 
in  its  final  form  as 


h 


c 


in  which 

o 


.  8-to/i3 


1  + 


k-1  2 


oj/5 


(28) 


(29) 


The  main  point  of  interest  for  the  present  investigation  is  that  Bartz 
based  his  derivation  on  the  Dittus-Boelter  form  of  correlation  for 
fully  developed  turbulent  pipe  flow,  i.e. 

Nu  =  C  (Re)m  (Pr)n  (30) 

Emmons  (13)  developed  a  one  dimensional  analysis  for  cooling 
with  an  evaporating  liquid  film  in  fully  developed  turbulent  pipe  flow. 
It  was  based  on  the  assumptions  previously  listed  in  Chapter  I.  He 
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started  with  th«  following  basic  equations  for  a  turbulent  in- 
eosq>resslble  boundary  layer: 

Continuity 


(31) 


Momentum 


Energy 


Diffusion 


—  3u  ,  —  3u  3 
u  T-  +  ▼  TT  "  TT 
3a  9y  Jy 


ug  +  75-4- 

3*  3y  3y 


-W.  -Wa  3 

u~ +  v  — "  17 


<U  +  «H>  £ 


3WS 

<*a  +  «D>  17 


(32) 


(33) 


(34) 


The  barred  parameters  In  these  equations  represent  time  averaged 
quantities,  for  example: 

_  r  i  „ 

u(t)  dc 


=■/  h 

T-*~ 


(35) 


The  bars  are  now  dropped,  with  the  understanding  that  these  para¬ 
meters  are  time  averaged, 

Emmons  simplified  Eqs,  (31),  (32),  and  (33)  by  assuming  that 
terms  Including  gradients  in  the  x-dlrection  are  negligible  compared 
with  terms  involving  gradients  in  the  y-directioa.  From  the  continuity 
equation  this  gave  v  -  constant.  Humans  evaluated  this  constant  by 
assuming  that  the  gas  on  the  surface  of  the  liquid  was  pure  coolant 
vapor,  thus 


(36) 
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The  mo— ntm  and  energy  aquation*  were  then 
Mo— ntum 


Energy 


ff-tf-v  S] 


(37) 


(38) 


At  thl*  point  Emmons  Introduced  the  expression  for  the  eddy 
dif fuslvity  of  momentum  a*  developed  by  Rennie  (29)  and  modified  by 
Turcotte  for  the  case  of  wall  Injection  (40)* 

*M  *  »  ,lnh2(T^>  (39> 

Substituting  from  Eq.  (39)  Into  Eq.  (37)  and  integrating  assuming 
zero  velocity  at  y  »  0  gave 


V 

p  Q 


t“h !®) 

'  w  w 


(40) 


From  Turcotte' s  analysis  (40) 


t  -r  . 
v  ni 
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(41) 
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where 


0 


13.89 


* 

byu 


(43) 


and 


* 


u 


(44) 


Emmons  then  used  an  extended  Reynolds  analogy,  with  Pr  •  1 


“d  «H"«M 


to  get  the  temperature  distribution  as 


T 
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p  u  Q 

w  gx 


ja 

exp(0Qe  2 
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(45) 


Use  of  this  temperature  distribution  gave  the  following  equation  for 
the  Stanton  Number 


JS 

2  ni  u 


(46) 


Emmons  then  noted  that  for  the  case  u  ■  u  and  Q  -  0  this  reduces  to 

m  g 

the  result  obtained  by  the  simple  Reynolds  analogy.  He  then  proceeded 
to  use  the  heat  transfer  coefficient  obtained  from  Eq.  (45),  specifically 


h 


kHTni  -8Q 
0 

vu8 


(47) 


without  further  modification. 
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Ie  was  noted  by  the  present  Investigator  that  the  exponential 
tera  In  Eqs.  (46)  end  (47)  could  be  looked  upon  as  a  simple  correction 
term  applied  to  the  non-injection  heat  transfer  coefficient  to  take 
Into  account  the  fluid  injection,  either  by  e vapor htl on  or  by  trans¬ 
piration  at  the  surface. 

There  ere  better  correlating  equations  available  for  predicting 
the  convective  heat  transfer  coefficient  without  injection  than  that 
which  coma  from  the  simple  Reynolds  analogy.  One  of  thesh  is  the 
Dlttus-Eoeltsr  equation  for  fully  developed  turbulent  pipe  flow,  the 
form  of  which  was  used  by  Barts  In  his  analysis  (2) ,  It  was  thus 
decided  to  program  e  stepwise  finite  difference  analysis  for  digital 
computer  solution  using  the  Dlttus-Boalter  correlation  at  each  step 
but  aubstltuclng  the  coefficient  "C"  ee  determined  by  Berts.  The 
convective  heat  transfer  coefficient  thus  determined  would  then  be 
modified  using  Easons'  "exponential  correction"  factor.  The  program 
also  provides  for  the  direct  use  of  experimentally  determined  heat 
transfer  rates  without  injection  and  modifying  these  using  Emmons- * 
"exponential  correction"  factor. 

Emmons1  analysis  assumed  that  the  surface  of  the  coolant,  film 
wee  at  saturation  temperature  immediately  upon  injection.  If  the 
coolant  has  a  high  heat  capacity  and  high  saturation  temperature  this 
is  not  e  goad  assumption,  even  for  very  high  heat  transfer  rates.  It 
wee  thus  decided  to  forego  this  assumption  end  calculate  the  temperature 
profile  of  the  liquid  film  starting  at  the  Injection  point  and  taking 
into  account  the  evaporation  of  the  film  while  it  Is  increasing  in 
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temperature  up  Co  the  saturation  temperature. 

In  order  to  calculate  the  developing  temperature  profile  it  is 
necessary  to  assume,  or  calculate,  the  velocity  profile.  Since  the 
ratio  of  the  dlffusivlty  of  momentum  to  the  diffuslvity  of  heat,  l.e*, 
Pr  Is  >  3  at  room  temperature  for  water,  the  main  coolant  liquid  under 
consideration,  It  was  decided  to  assume  instantaneous  development  of 
the  velocity  profile  upon  injection.  A  condition  of  constant  shear 
stress  across  the  liquid  film  was  also  assumed.  This  latter  assumption 
appaars  to  be  reasonable  for  such  a  thin  film  aa  that  used  for  film 
cooling. 

Consider  a  fluid  layer  as  shown  in  Figure  24,  with  mass  flow  rate 
w.  Now  for  constant  properties  and  constant  shear  stress  across  this 
layer 


n+1 


(48) 


and 


U_.  .  +  U« 

*  -  2lV“u  (--2— S) 
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Combining  Eqa.  (A 8)  and  (49)  and  rearranging 
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The  positive  value  on  the  radical  Is  used  since  d&  oust  be  positive. 
Starting  with  a  known  value  of  uft  and  estimating  initial  values  of 
P*  M»  andr  Eq.  (51)  may  be  solved  for  d  ,  then  u. . «  calculated  from 

W  u  Qtl 

Eq.  (19) .  As  better  average  values  of  temperatures  become  known  from 
the  iterative  solution  of  the  heat  balance  equations  p  and  p  are 
corrected  and  the  velocity  profile  recalculated.  It  is  not  necessary 
to  ra-evaluate  the  velocity  profile  for  each  temperature  iteration 
as  it  produces  only  vary  slight  effects.  However,  when  the  con¬ 
vergence  criterion  for  the  temperature  Iteration  is  satisfied,  the 
velocity  profile  should  be  corrected  and  convergence  checked  again. 

This  assures  the  dse  of  values  of  p  and  u  which  are  temperature 
corrected. 

The  aquations  necessary  for  the  temperature  profile  calculation 
are  derived  by  taking  e  heat  balance  on  the  nth  layer.  The  heat 
balnea  for  steady  state  requires  that: 

Hat  Heat  Rate  Convected  In  +  Net  Heat  Rate  Conducted  In  ■  0  (52) 

Referring  to  Figure  24,  page  69. 

Net  Heat  Rata  Convected  In  ■  w.c  (T1  -  T2  )  (53) 

nun  n  ' 

Hat  Hast  Rate  Conducted  In  -  2trr& 
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Let 


<T1n±I  +  TW 


a+1 


(55) 


and  similarly  for  and 
Also  let 


2irr 


rr+1 


n+1  d  ,  d 

(-S±L-  +  - 2-) 
'‘2k. .  2k  ; 
n+1  n 


(56) 


and  similarly  for  and  Kn_^. 


Now  substituting  from  Eqa.  (53),  (54),  (55),  and  (56)  into  (52) 


V.  <T1n  *  TV  -  4WVl  -  V  +  “.-lff,-rV-°  «7) 


Resubstituting  from  Eq.  (55)  for  T  and  solving  for  T2 

n  n 


T2  - 
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Tl  Tw  c  -  +  K  ,)1  +  4(K  .T  . .  +  K  -T  *) 

n  L  n  n  2  n+1  n-l  J  n+1  n+1  n-1  n»l# 


wnCn  +  f  (Kn+l  «n-i> 


(58) 
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Now  Eq»  (58)  applies  for  the  Interior  layers  of  the  film  but 
ocher  equations  must  be  developed  for  the  heat  balance  on  the  stream" 
aide  and  wall-aide  layers „  For  the  stream-aide  layer 

V  c  CT1  -  T2  )  +  AK,  ,  (T  ,  -  T  )  +  2irr  Ah(T  -  ?  )-2rr  Oh.  -0  <59> 

8  8  8  8  8—1  8—1  8  a  8  8  fg 
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w  c  -  “(K  ,+2trr  h)l  +  <K  ,T  ,+2rr  hT  )-2»r  Qh, 
,8  8  2  3-1  8  J  8-1  8-1  8  °*  8^  fg 


(60) 


Vs  +  2<K8-1+2  r8h> 


For  the  well-side  layer 

V.(n,  -  TV  +  -  V  +  4VTb  -  V° 
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Thus 
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(63) 


Iterative  calculation  of  Q  la  described  starting  in  paragraph  2,  page 
77 4  The  coefficient  h  is  determined  from' the  modified  Dlttus-Boelter 
equation  which  will  subsequently  be'  discussed  in  detail.  The  co¬ 
efficient  h.  is  determined  from  the  backside  coolant  flow  conditions. 
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Since  the  backside  coolant  flow  will  be  turbulent  for  practically 
all  cases,  a  regular  Dittus-Boelter  calculation  will  normally  be  used. 

Application  of  Eqs,  (58),  (60),  and  (63)  across  the  film  results 
in  a  set  of  nonlinear  simultaneous  algebraic  equations.  The  non¬ 
linearity  comes  from  the  dependence  of  c,  d,  and  k  on  T,  Since  this 
dependence  is  relatively  weak,  however,  the  equations  may  be  solved  by 
the  Gauss-Seldel  iteration  method  with  periodic  corrections  of  c,  d,  and 
k  for  temperature  change.  This  is  done  in  the  same  manner  as  the 
periodic  correction  of  P  and  P  in  the  velocity  profile  calculation. 

The  temperature  profile  is  first  estimated.  Substitution  is  then  made 

into  the  explicit  equations  for  T2  , - ,  T2  and  the  amount  of  change 

w  s 

in  T2  checked  against  some  convergence  criterion  until  it  la  satisfied 
n 

for  all  the  layers.  Periodically  c,  d,  and  k  are  re-evaluated,  with 

a  final  evaluation  made  as  a  check  when  T2  has  converged. 

n 

Other  methods  could  be  used  for  the  solution  of  Eqs.  (58),  (60), 
and  (63).  For  example,  by  rearrangement  they  could  be  put  into  the 
matrix  form: 

A  T  -  B  (64) 

By  estimating  T,  A  could  be  calculated.  Then 

T  -  A~l  3 

A  could  then  be  recalculated  using  the  new  values  for  T  and  the  method 
repeated  until  T  converges.  The  Gauss-Seidel  iteration  method  was' used 
because  the  coefficient  matrix  la  sparsely  populated  and  this  method 
does  not  require  manipulation  of  the  zero  elements. 
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The  incremental  film  temperature  calculation  la  continued 
a*  described  above  until  T2  >  T  .  Normally  after  some  particular 

8  ®  SAC 

incremental  calculation  T2  >  T  „ .  In  order  to  determine  the  position 

s  sat 

where  T2fl  ■  an  interpolation  Is  made*  It  was  found  that  if  the 
initial  Incremental  steps  are  small  a  straight  line  interpolation 
here  is  sufficiently  accurate  without  further  iterative  calculation 
of  T2  . 

The  evaporation  rate  Q  must  be  calculated  both  in  the  region 
before  T2  •  T  and  following  this  region,.  The  only  difference  be* 
tween  the  two  calculations  is  that  in  the  region  before  saturation 
temperature  is  reached  the  rate  of  evaporation  is  controlled  by  the 
available  partial  pressure  of  the  film  vapor  and  the  mass  transfer 
coefficient  at  the  surface,  whereas  in  the  region  following  the 
rate  of  evaporation  is  controlled  by  the  excess  thermal  energy 
deposited  in  the  film  by  heat  transfer  from  the  mainstream.In  both 
cases  the  heat  transfer  coefficient  h  must  be  corrected  for  the 
evaporation  effect  and  the  loss  of  liquid  film  due  to  evaporation  must 
be  calculated.  Examination  of  these  effects  shows  that  they  are 
coupled  and  thus  require  on  iterative  solution*  The  evaporation  rate 
for  a  given  evaporation  enthalpy  and  temperature  difference  varies 
with  h.  On  the  other  hand,  the  following  equation  for  the  heat  trans¬ 
fer  coefficient  shows  the  dependence  of  h  on  the  evaporation  rate* 

h  »  ,0265  ^  (Re) 08  (Pr)  ‘333pxpr13ffi-Q)j  T65) 

In  Eq,  (65)  Q  is  the  mass  evaporation  rate  of  the  film  coolant. 

In  the  region  before  saturation  temperature  is  reached  Q  may  be  cal¬ 
culated  directly  using  Eq.  (69),  page  75  ,  and  the  vapor  partial  pressure. 
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For  tbe  rag ion  following  the  attainment  of  saturation  temperature,  here¬ 
after  called  the  boiling  region,  there  la  a  simple  relationship  between 


Q  and  h. 


qhfgA  =.  hA(TAdw 


T1  +  T2 
_ a _ a 


')  -  Q. 


In  Emmons'  analysis  it  was  assumed  that  the  fluid  in  the 
boundary  layer  over  the  liquid  surface  Is  pure  coolant  vapor  and  that 
the  thermal  conductivity  in  Eq .  (68)  should  be  that  of  the  pure 
coolant  vapor  evaluated  at  T  ,  It  has  been  demonstrated  (12)  that 

Sow 

constant  property  solutions  for  skin  friction,  Nusselt  number  and 

recovery  factor  nay  be  used  even  when  large  variations  In  physical 

properties  occur  provided  the  properties  are  evaluated  at  a  reference 
* 

temperature  T  which  is  expressed  explicitly  as 

T*  •  Tg  +  0.5(Tw  -  T  )  +  0.22(Tr  -  T  )  (69) 

* 

It  can  be  readily  seen  that  T  la  the  widely  used  film  temperature 


(or  average  of  f  and  T  )  plus  s  term  to  account  for  temperature 

8 

recovery  In  high  speed  flow;  This  term  seems  to  effect  s  good  property 


correction  even  for  dissociation  effects «  In  the  present  solution 

* 

fluid  properties  are  evaluated  at  T  , 
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If  ehs  thermal  conductivity  of  the  binary  mixture  differs  gnat* 
ly  from  that  of  tha  coolant  vapor  an  effort  abould  ba  mads  to  evaluate 
^  for  tba  binary  mixture.  If  tha  eencantratlon  of  coolant  vqpor  at 
tha  liquid  film  surf ape  can  ba  calculated  tha  thermal  conductivity  for 
a  gas  mixture  at  half  this  concentration  should  ba  a  batter  value  to 
usa  than  that  of  tha  pun  coolant  vapor.  This  conductivity  should  ba 
evaluated  at  T*  as  Indicated  above. 

Examination  of  Eqs. (32) ,(33) «  and  (34)  shows  that  the  same 
analogy  may  ba  made  between  mass  transfer  and  momentum  transfer  as  that 
which  was  mads  between  heat  transfer  and  momentum  transfer.  Experimental 
nsults  indicate  that  by  use  of  tha  proper  form  for  the  dlffuslv£ty  of 
mass  and  use  of  the  Schmitt  number.  Sc,  in  piece  of  Pr,  a  coefficient 
of  mass  transfer  for  turbulent  flow  may  be  calculated  ae 


(Sc) 


.333 


(70) 


It  is  therefore  proposed  that  in  accordance  with  tha  above  mentioned 
analogy  Eq.  (-70)  be  modified  and  written  in  the  fora  of  Eq.  (68)  ** 


(Sc) 


■13.890^ 

byu*  J 


(71) 


If  the  total  coolant  Injaetlon  rate  is  only  e  smell  pareqatage  tha 
main  gaa  flow  rata  any  concentration  of  coolant  in  the  main  stream 
may  ba  naglectad  so  for  this  '’•ase 


How 


T  ♦  T 
'c  'a 


(72) 

(73) 
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Substituting  Eqs.  (72)  and  (73)  Into  Eq.  (71)  gives 


.02651^ 

D 


(Re)48  (Sc)*33txp 


-13.890 

bu*(h^'HV 


(74) 


The  Iterative  procedure  for  determining  h  is  then  as  follows: 

* 

1.  Calculate  and  Sc  at  T  , 

2.  Calculate  Q  based  on  value  of  h  for  previous  Increment, 

3.  Starting  with  Q  from  step  2  iteratively  calculate  h  using 

H 

Eqs,  (71),  (72),  (73),  and  (74). 

4.  Calculate  h  from  Eq,  (65) , 

5.  Calculate  Q  from  Eq.  (67), 

6.  Calculate  from  Eq,  (72). 

7.  Re-evaluate  kg  based  on  y  from  step  6. 

c 

8.  Re-evaluate  h  using  Eq.  (68)  and  compare  with  h  from  step  4, 

9.  Repeat  steps  3  through  7  until  h  converges. 

10.  Re-evaluate  Q  and  compare  with  Q  from  step  2. 

11*  Repeat  iteps  3  through  10  until  Q  converges. 


Two  additional  problems  are  encountered  in  the  f.  on  -dure  outlined 
above.  One  of  these  is  the  determination  of  t.ie  thermal  conductivity  of 
a  binary  gas  mixture  at  a  given  temperature  and  the  other  is  the 
determination  of  the  dlffusivlty  of  mass  for  the  coolant  vapor  through 
the  mainstream  gas  for  a  known  pressure  and  temperature. 


Thermal  conductivities  of  gas  mixtures- may  be  estimated  from 

Q 


mix 


Vi 


(75) 
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whs  re  ths  x  and  are  oole  fractions  end  thermal  conductivities, 
respectively,  of  the  pure  gas  species* 

Given  the  dif fusivity  of  mass  for  a  binary  gas  mixture  at 
oile  condition  of  pressure  and  temperature  the  dif fusivity  at  other 
pressures  and  temperatures  may  be  estimated  from  the  following 
equation  (33) 


c  t3//2 

P 


(76) 


thus 


where  C  is  a  constant; 


(dab)t2,p2  “  (dab5t,p(t  j  1  (p^ 


(77) 


In  the  present  study  Eq.  (77)  was  used  to  predict  mass  dlffusjLy- 
itiea  at  the  local  reference  temperature  T*  and  the  locpl  static 
pressure*  dependence  of  the  dif fusivity  on  concentration  was  neglect¬ 
ed* 

Now  for  the  diffusion  of  mass 


dPA 

lk  "  "DAB  dy 


(78) 


where  1^  Is  the  mass  current  in  lbm«/ft.Sir.  and  dp^/dy  is  the 
density  gradient  in  lbm./ft,  ,  In  forced  convection  mass  transfsr 
the  diffuslvity  0^  plays  the  same  rola  as  that  played  by  the  thermal 
conductivity  In  conduction  boat  transfer.  Eq.  (71)  may  thus  be 
rewritten  as 
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*m 


.0265  °AB 
D 


(**r8(sc) 


rl3,«|.V 

bu 


(79) 


Is  order  to  use  the  above  equations  at  each  axial  position  the 
mainstream  velocity,  temperature,  and  density,  must  be  calculated. 
Froxen  flow  through  the  nozzle  is  assumed  with  a  specific  heat  ratio 
determined  by  the  bulk  assn  temperature  at  the  nozxlst  entry.  The 
aquation  which  relates  the  local  diameter  to  throat  diameter  ratio, 
specific  heat  ratio,  and  Mach  number  is 


Jfrl . 

2  (k— 1) 


(80) 


Now  (D/Dt>  and  k  are  known  and  M  is  to  be  calculated.  Since  Eq.  (goy 
cannot  be  solved  explicitly  for  H  soma  iterative  procedure  must  be 
used  to  solve  for  M„ 

To  solve  Eq.  (go)  for  M  the  following  rapidly  converging 
Iterative  procedures  were  devised.  Designate  the  Mach  numbers  in 
Eq0  (80)  differently  by  subscripts  as  follows: 


Now  for  Iteration  in  the  subsonic  portion  of  (tie  nozzle,  1.  e., 
upstream  from  the  throat,  Eq.  <5L)  Is  solved  explicitly  for  Kft, 
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H 


a 


it+i 


2(k-l) 


<h2 


(82) 


Tha  Iteration  may  be  started  by  substituting  any  cumber  from  0 
to  1  Into  Eq,  (82)  for  M^.  The  value  of  calculated  Is  then 
substituted  back  in  for  and  the  procedure  repeated  until  the 
difference  between  and  meets  a  specified  criterion*  For  the 
supersonic  portion  of  the  nozsle,  1.  a. ,  downstream  from  the  throat, 
Eq«  (81)  Is  solved  explicitly  for 


(83) 


This  Iteration  la  started  by  substituting  1  or  a  number  greater  than 
1  into  Eq.  (  81)  for  After  Is  calculated  it  is  substituted  back 
In  for  and  the  procedure  repeated  as  described  above. 

Vlth  the  Mach  number  calculated  the  local  pressure,  temperature 
and  density  for  the  mein  stream  gas  may  be  calculated. 

The  general  caleulational  procedure  for  liquid  film  cooling 
proceeds  as  follows: 

1.  For  the  known  axial  position  calculate  the  ratio  0/Dg. 

2.  Doing  Eqe,  (82)  end  (83)  end  the  iterative  procedure 
described  calculate  the  Mach  number. 
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3.  Using  Che  Mach  number  calculate  the  local  (free  stream) 
pressure!  temperature  and  density, 

4,  Evaluating  properties  at  the  reference  temperature 
calculate  Re,  Prt  Sc,  and  D^, 

S'.  Calculate  T  .  and  T  .  for  the  film, 
ni  sat 

6.  Check  to  see  If  the  film  surface  is  up  to  T  .  If  not 

8  Qb 

calculate  the  film  velocity  and  temperature  profiles 
correcting  for  evaporation  effects  and  go  to  step  (9), 

If  the  film  surface  is  at  or  above  T  go  to  step  (7), 

7.  Set  film  surface  temperature  equal  to  T  and  calculate 
film  velocity  and  temperature  profiles, 

8.  Calculate  h  iteratively  correcting  for  the  change  in 
thermal  conductivity  due  to  mass  transfer. 

9.  Subtract  the  amount  of  film  coolant  evaporated  over  the 
axial  increment  from  the  maas  flow  rate  In  the  film  surface 
layer,  (Soma  fraction  of  the  layer  is  set  as  a  maximum 
amount  to  be  evaporated  per  axial  Increment  and  the 
lncreaant  else  reduced  if  thla  la  exceeded) . 

10.  Increment  the  axial  position  and  return  to  step  (1), 

The  above  procedure  has  been  programmed  in  Fortran  IT  language  for 
digital  computer  execution  and  the  program  is  included  as  Appendix  B, 

The  above  procedure  may  also  be  applied  for  calculation  of 
coolant  evaporation  rater  if  experimental  heat  transfer  values  without 
film  cooling  are  available.  If  the  flow  condition  is  turbulent,  so 
that  Turcot  to’ s  formula  for  c^(  Eq,  (?9),  la  a  good  approximation 
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f— rm 1  exponential  correction  f Actor  should  ba  direct ly  applicable. 

The  procedure  la  aa  follow*} 

lf2l3,4l5*  Tbeae  atepa  are  tba  aaaa  aa  for  calculation  of  the 
beat  load  fro*  • trees  variable*. 

6.  Check  to  aae  if  the  fils  a or face  la  op  to  K  not  use 

the  experimental  heat  load,  ^  and  calculate  film  velocity  and 
tasperatnra  profiles  correcting  for  evaporation  effects  and 
go  to  step  (9).  If  the  film  surface  la  -at  ox  above  Tft^t  go 
to  atap  <<77. 

7*  Sat  the  file  surface  tespavatura  equal  eo  and  calculate 
the  file  velocity  end  tenparature  pro tiles. 


8,  Using 


Q  - 


”  A 

~v 


(Si) 


and 


q 


^ex 


(85) 


calculate  the  base  load  to  the  file,  q,  Iteratively 
correcting  fob  the  change  in  the  real  conductivity 
due  to  ness  transfer^ 

9,  10*  Sene  as  for  direct  calculation  fron  stress 
variables* 
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II.  DISCUSSION  OF  ANALYTICAL  RESULTS 

The  above  analysis  of  liquid  film  cooling  was  programmed  in 
Fortran  IV  language  for  execution  on  a  digital  computer.  The  program 
la  Included  as  Appendix  B.  Since  experimental  studies  of  this  type 
cooling  were  not  to  be  performed  a  test  case  was  postulated  in  order  to 
check  out  the  analysis.  In  this  teat  case  the  nozzle  configuration, 
mainstream  stagnation  conditions,  and  the  backside  cooling  heat  transfer 
coefficient  were  similar  to  those  for  the  experimental  studies  of  gas 
film  cooling.  It  was  found  that  a  finite  distance  13  required  for  the 
liquid  film  to  rise  to  saturation  temperature,  and  that  a  significant 
fraction  of  the  liquid  evaporates  during  this  interval.  In  the  case 
chosen  this  amounted  to  approximately  30  per  cent  of  the  injected  film 
mass  flow  rate.  Figure  25  shows  the  increase  of  the  film  surface  tempera¬ 
ture  and  the  backside  wall  temperature  as  the  film  progresses  through 
the  nozzle.  The  increased  heat  flux  at  the  throat  and  the  rapidly  de¬ 
creasing  pressure  cause  the  liquid  at  the  surface  to  reach  the  satura¬ 
tion  temperature  just  downstream  from  the  throat  for  this  particular 
case.  The  liquid  surface  then  remains  at  the  saturation  temperature 
in  the  diverging  portion  of  the  nozzle  as  the  pressure  continues  to  drop, 
and  the  rate  of  film  evaporation  is  controlled  largely  by  the  heat  trans¬ 
fer  rate  available  from  the  mainstream  gas.  The  sample  computer  calcu¬ 
lation  for  this  case  is  given  in  Appendix  F. 
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CHAPTER  IV 

TRANSPIRATION  COOLING 
I.  THEORETICAL  ANALYSIS 

Transpiration  Based  on  Emmons 1  Liquid  Film  Analysis .  A  careful  compar¬ 
ison  of  liquid  film  cooling  with  transpiration  cooling  shows  that  they 
differ,  in  only  two  ways,  namely,  in  transpiration  cooling: 

1.  There  is  nc  liquid  film  to  be  transported  or  to  absorb  heat. 

2.  The  coolant  flow  rate  may  be  controlled  locally  independent¬ 
ly  of  the  heat  transfer  rate. 

The  basic  analysis  for  liquid  film  cooling  in  the  locations  where  boil¬ 
ing  occurs  may  thus  be  utilized  for  transpiration  cooling. 

In  the  caae  of  transpiration  cooling  with  a  liquid  it  is  neces¬ 
sary  for  the  liquid  to  be  vaporized  at  the  surface  or  before  reaching 
the  surface.  The  beat  procedure  is  to  specify  a  local  surface  temp¬ 
erature  which  is  equal  to  or  greater  than  the  local  saturation  tempera¬ 
ture.  The  parameter  to  be  calculated  is  the  coolant  mass  flow  rate 
per  unit  area  required  to  maintain  the  surface  at  the  specified  temp¬ 
erature.  This  may  be  done  iteratively  using  the  following  procedure: 
1,2, 3, 4.  Same  as  for  liquid  film  cooling. 

5.  Calculate  T  and  set  the  surface  temperature  either 
equal  to  T  or  to  some  specified  temperature  higher 

Saw 

than  T 

sat 

6.  Starting  with  the  heat  transfer  rate  without  transpira¬ 
tion  use  this  to  calculate  the  required  rate  to  hold 
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the  surface  at  the  set  temperature. 

7.  Correct  the  heat  transfer^ rate  for  transpiration  effects 
using  the  procedure  outlined  for  liquid  film  cooling. 
Recalculate  the  required  transpiration  rate.  Repeat  this 
step  until  either  the  change -in  heat  transfer  rate  or  the 
change  in  transpiration  rate  meets  the  set  convergence 
criterion. 

8.  Increment  the  axial  position  and  return  to  step  (1). 

The  procedure  is  identical  for  gas  transpiration  cooling  except  the 
calculation  of  T  is  not  required.  It  has  been  programmed  in  Fortran 

Sflb 

IV  language  for  execution  on  a  digital  computer.  The  program  Is  In¬ 
cluded  in  Appendix  0. 

The  same  type  of  calculation  may  be  made  using  heat  transfer 
rates  determined  experimentally  without  transpiration  cooling.  The 
changes  are  the  same  as  those  for  liquid  film  cooling  except  the  calcu- 
1 at ion  of  film  evaporation  is-  not  required. 

Analysis  Based  on  Gas  Film  Cooling.  The  analysis  given  for  gaa  film 
cooling  may  be  used  for  gas*  transpiration  cooling  provided  It  is  assumed 
that  the  -same  model  may  be  used  in  the  calculation  of  eddy  diffuslvitles 
for  both  types  of  cooling.  This  is  done  by  specifying  a  very  small  axial 
increment  for  the  calculation  and  aj  small  amount  of  coolant  fluid  to  be 
added  at  each  increment.  A  new  coolant  stream  layer  is  added  whenever 
second  coolant  layer  from  the  nozzle  wall  contains  twice  as  much 
fluid  as  it  originally  contained.  The  calculation  is  performed  assuming 
an  adiabatic  wall. 
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It.  TRANSPIRATION  CALCULATION  RESULTS 

Since  experimental  studies  of  transpiration  cooling  were  not 
to  be  performed  test  cases  were  postulated  in  order  to  check  out  the 
analyses . 

The  results  of  the  calculation  using  the  Liquid  or  Gas  Transpira¬ 
tion  Cooling  Program  are  shown  graphically  in  Figure  26.  Note  that  as 
the  calculation  proceeds  into  the  nozzle  an  increasing  coolant  flow  rate 
is  required  to  hold  the  nozzle  wall  at  500  °F.  The  calculation  was  con¬ 
tinued  only  far  enough  to  insure  proper  operation  of  the  program  and 
convergence  of  the  iterative  procedures.  The  sample  computer  calculation 
for  this  case  is  given  in  Appendix  F. 

The  results  of  a  transpiration  calculation  using  the  Gas  Film 
Cooling  Program  are  shown  graphically  in  Figure  27.  Note  the  shape  of 
the  wall  temperature  profile  curve  for  a  fixed  injection  rate  at  each 
increment.  The  maximum  temperature  occurs,  as  expected,  very  near  the 
throat.  The  sample  computer  calculation  for  this  case  is  given  in 
Appendix  F. 
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Fig-  26  Graph  of  Sample  Transpiration 
Calculation  Results  Using 
Liquid  or  Gas  Transpiration 
Cooling  Program 
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Fig-  27  Graph  of  Sample  Transpiration  Calculation  Results  Using  Gas  Film  Cooling  Program 
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CHAPTER  V 

CONCLUSIONS  AND  RECOMMENDATIONS 

Based  on  the  experimental  results  presented  and  the  develop¬ 
ment  of  satisfactory  calculatlonal  techniques  to  Implement  the 
analyses  given  in  the  present  study,  the  following  conclusions  can 
be  drawn: 

1.  Gas  film  cooling  can  be  used  to  measurably  lover  the  vail 
temperatures  and  the  wall  heat  fluxes  in  the  converging  section  and 
at  the  throat  of  a  high-pressure  high-temperature  nozzle.  For  thlB 
nozzle  a  change  of  gas  film  coolant  injection  ratio  from  5%  to  10|f 

produced  a  measurable  decrease  In  the  wall  temperatures  and  the  wall 

\ 

heat  flux  up  to  a  point  near  the  throat.  Higher  gas  film  cooling 
ratios,  such  as  30% ,  lower  the  temperature  and  heat  flux  in  the  con¬ 
verging  section  and  at  the  throat  appreciably  and  have  a  measurable 
effect  extending  into  the  supersonic  expansion  section.  Some  gross 
mixing  occurs  at  the  injection  point,  the  exact  amount  being  as  yet 
some  undetermined  function  of  the  Injection  geometry,  the  relative 
velocities  of  the  main  gas  stream  and  the  coolant  stream  at  the 
injection  point,  and  the  entering  velocity  profiles  and  turbulence 
conditions. 

2,  A  straightforward  boundary  layer  type  analysis  has  been 
developed  and  programmed  which  predicts  with  reasonable  accuracy 
the  nozzle  wall  temperatures  and  wall  heat  fluxes  in  the  converging 
section  and  at  the  throat  for  gas  film  cooled  nozzles.  The  calculated 
wall  temperatures  and  heat  fluxes  using  the  full  amount  of  film  coolant 
Injected  In  the  experiment  and  the  values  calculated  assuming  no  film 
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coolant  bracketed  the  measured  values  of  temperature  and  heat  flux 

* 

for  all  cases  correlated.  The  analysis  must  be  modified  in  the  super¬ 
sonic  portion  of  the  nozzle  if  accurate  prediction  of  wall  temperatures 
and  heat  fluxes  is  required  in  this  region. 

3.  Calculational  techniques  have  been  developed  and  programmed 
for  predicting  the  effectiveness  of  liquid  film  cooling  and  liquid  or 
gas  transpiration  cooling  in  nozzles.  These  techniques  should  be 
checked  against  experimental  data  before  they  are  used  for  design 
purposes,  however. 

It  is  recommended  that  gas  film  cooling  experiments  In  super¬ 
sonic  nozzles  be  performed  with  various  Injection  geometries  and  main¬ 
stream  to  coolant  flow  velocity  ratios  in  order  to  evaluate  the 
amount  of  gross  mixing  which  occurs  at  injection  as  some  function  of 
these  parameters. 

It  is  also  recommended  that  liquid  film  cooling  experiments  be 
performed  on  a  supersonic  nozzle  operating  at  high  stagnation- pressures 
and  temperatures  in  order  to  evaluate  the  accuracy  of  predictions  using 
the  liquid  film  cooling  analysis  In  the  present  investigation. 

It  is  further  recommended  that  nozzles  fully  cooled  up  to  the 
throat  by  transpiration  be  constructed  and  tested  experimentally  in 
order  to  evaluate  the  accuracy  of  predictions  using  the  analyses 
presented  in  the  present  investigation. 

* 

A  possible  modification  uses  information  as  shown  in  Fig.  7-58  of 
Turbulence  by  J.  0.  Hinze,  a  text  published  by  McGraw-Hill. 
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APPENDIX  A 
NOTATION 


A 

A, 

a 

A 

B 

b 

C 


\ 

j 


nozzle  cross  sectional  area  or  incremental  heat  transfer  area 
nozzle  throat  cross  sectional  area 
two-dimensional  coefficient  matrix 
column  matrix 

Turcotte's  constant  to  correct  diffusivity  for  mass  injection 
constant 

specific  heat  at  constant  pressure 

sonic  velocity  at  nozzle "throat 

heat  capacity  of  fluid  in  layer  at  wall 

nozzle  Inside  diameter 

binary  mass  diffusivity 

nozzle  throat  inside  diameter 

nozzle  throat  inside  diameter 

liquid  coolant  layer  thickness 

local  nozzle  liner  thickness 

friction  factor  with  no  injection  or  evaporation 
gravitational  constant 
convective  heat  transfer  coefficient 
enthalpy  of  evaporation 

coefficient  of  mass  transfer 
Joules  constant  =  778  ft  lh/Btu 

ratio  of  specific  heats,  1  /c^ 
nozzle  liner  thermal  conductivity 
thermal  conductivity 
"conductivity”  of  mass 
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M 


M 


A 

Nu 


?2 

Pr 

Q 

Qs 

Re 


I 

Sc 

St 

T 

T_ 


aat 


Adw 


u„ 


Mach  number 

molecular  weight  of  species  A 

molecular  weight  of  species  B 

Nusselt’s  number 

stagnation  pressure 

pressure  at  increment  entrance 

pressure  at  increment  exit 

Prandtl's  number 

liquid  film  evaporation  rate 

energy  stored  in  increment  of  liquid  film 

Reynold's  number 

radius 

radius  of  curvature  of  nozzle  throat 

inside  radius  of  nozzle 

Schmidt  number 

Stanton's  number 

temperature 

temperature  of  gas  stream 
stagnation  temperature 
recovery  temperature 

temperature  at  liquid  film  surface 
saturation  temperature,  coolant  vapor 

wall  temperature 

column  temperature  matrix 
adiabatic  wall  temperature 

time 

time  averaged  velocity  in  x-direction 
friction  velocity  A/p 

* 

dimensionless  velocity  u/u 
gas  stream  velocity  at  throat 
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u 

s 

u 

m 

V 

W 

a 

w 

w 

w 

x 

y 

+ 

y 

GREEK 


a 

A 

5 

a 

€ 

£D 

eh 

eM 

Y 


Y 


C 


Y 


3 


V 


U 


P 


P 


w 


T 


ni 


T 


W 


gag  stream  velocity 
stream  mean  velocity 

time  averaged  velocity  in  the  y-direction 
time  averaged  density  of  molecular  species  a 
mass  flow  rate 

mass  flow  rate  of  coolant  layer  at  the  wall 

axis  parallel  to  nozzle  wall 

axis  perpendicular  to  nozzle  wall 

: k 

dimensionless  y  variable  yu  /u 


molecular  diffugivity  of  heat 

increment  length  of  liquid  film  cooling  analysis 

molecular  diffusivlty  of  mass 

general  diffusivlty 

eddy  diffusivlty  of  mass 

eddy  diffusivlty  of  heat 

eddy  diffusivity  of  momentum 

gas  density  at  liquid-gas  interface 

coolant  gas  partial  density 

mainstream  gas  partial  density 

coefficient  of  viscosity 

molecular  diffusivlty  of  momentum 

density 

density  at  wall  or  surface 

wall  shear  stress  without  injection  or  evaporation 
wall  3hear  stress 


4)  ratio  of  diffusivity  of  heat  to  diffusivity  of  momentum 

<u  empirical  exponent  relating  viscosity  to  temperature 
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APPENDIX  B 

LIQUID  FILM  COOLING  PROGRAM 


C  FIRST  INPUT  CARD  40  BLANKS  FOR  CASE  TITLE. 

C  SECOND  INPUT  CARD 

C  FI — -F5  COEFFS  FOR  FILM  VI  SC. .LBM/FT.ScC. 

C  THIRD  INPUT  CARD 

C  El - E5  COEFFS  FOR  FILM  HT.  CAP. .B/LBM.F. 

C  FOURTH  INPUT  CARD 

C  CN1 - CN5  COEFFS  FOR  FILM  COND . , B/HR • FT «F * 

C  FIFTH  INPUT  CARD 

C  Gl - G5  COEFFS  FOR  STREAM  VlSC..  LBM/FT.SEC. 

C  SIXTH  INPUT  CARD 

C  R 1 - R5  COEFFS  FOR  FILM  DENS*  .  LBM/CU.FT. 

C  SEVENTH  INPUT  CARD 

C  VK1 - VK5  COEFFS  FOR  VAP.COND.,  B/HR.FT.F. 

C  EIGHTH  INPUT  CARD 

C  AK1 - AKS  COEFFS  FOR  STREAM  COND*.  B/HR.FT.F. 

C  NINTH  INPUT  CARD 

C  PR  1 - PR 5  COEFFS  FOR  MAINSTREAM  PRANDTL  NO. 

C  TENTH  INPUT  CARD 

C  P2  =  INLET  GAS  PRESSURE  P— ZERO « LB ./IN.  SO* 

C  TZ  55  INLET  GAS  TEMP.  T-ZERO,  DEGREES  F. 

C  DT  «  THROAT  DIAMETER  D*.  INCHES 

c  qi=initial  coolant  flow  rate.lbm/sec. 

C  TCI  =  INITIAL  COOLANT  TEMPERATURE.  DEGREES  F. 

C  ELEVENTH  INPUT  CARO 

C  GK  =  MAINSTREAM  RATIO  OF  SPECIFIC  HEATS* 

C  RHZ  =  INITIAL  MAINSTREAM  DENSITY.  LBM./CU.FT. 

C  B  =  TURCOTTES  CONSTANT  FOR  WALL  SHEAR  CALCULATION. 

C  RG  =  GAS  CONSTANT  R  AT  THROAT  CONDITIONS 

C  PPZ  =  INITIAL  PARTIAL  PRESSURE  OF  COOLANT.  PS  I  * 

C  TWELFTH  INPUT  CARD 

C  WCON  =  NOZZLE  LINER  THERMAL  CONDUCT I V I TY . B/HR.FT.F. 

C  WTMOL=  MOLECULAR  WEIGHT  OF  THE  FILM  LIQUID 

C  TCOM  =  CONVERGENCE  CRITERION  ON  TF  ITERATION 

C  T9N  =  CRITERION  FOR  THERMAL  BOUNDARY  LAYER.  DEG.  F. 

C  DIF  =  VAPOR-AIR  DIFFUSION  COEFF . ,  FT.SQ./SEC* 

C  THIRTEENTH  INPUT  CARD 

C  OX=  DISTANCE  INCREMENT.  INCHES. 

C  SMI =  DISTANCE  TO  NOZZLE  THROAT.  INCHES 

C  QOIV=MAX.  FRACTION  OF  LAYER  EVAP.  PER  SECTION 

C  RGC  =  COOLANT  VAPOR  GAS  CONST*  R 

C  fourteenth  INPUT  CARD 

C  M2  =  MAXIMUM  NUMBER  OF  INCREMENTS. 

C  KS  *  NO.  OF  TABLE  ENTRIES  FOR  D  AND  S. 

C  N  =  NUMBER  OF  FILM  LAYERS 

C  KPG  =  NUMBER  OF  TABLE  ENTRIES  FOR  TSAT.HF,  AND  HG. 

C  NEXT  K5  CARDS 
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C  S — IN**  0 — IN.*  WTHCK — IN.*  HW — 6/HR *SQ * FT . F . 

C  NEXT  KPG  CARDS 

C  ON  EACH  CARO  IN  ORDER  AS  FOLLOWS — 

C  PC — PSIA,  TSAT — DEG.  F« ,  HF — 8TU/LB*  HG — BTU/LB* 

C  OUTPUT  DESCRIPTION 

C  S  =  DISTANCE  AlONG  WALL  FROM  INJECTION*  INCHES. 

C  Q/A  =  WALL  HEAT  LOAD*  B/HR. SQ.FT* 

C  TB  “  LOCAL  OUTSIDE  WALL  TEMP . ♦ DEG .F . 

C  NBOIL  =  NUMBER  TO  SHOW  BOILING  CONDITION 

C  0  =  NO  BOILING 

C  1  BOILING 

C  H  a  MAINSTREAM  CONVECTION  COEFF.*  B/HR. SQ.FT *DEG.F. 

C  EXP9  =  EXPONENTIAL  ATTENUATION  FACTOR 

C  TF2(NI>  =  TEMP.  OF  INSIDE  COOLANT  LAVER.  DEG.F. 

C  N I  —  NUMBER  OF  INSIDE  COOLANT  LAYER 

100  FORMAT  (1H  *  5 (E 1 2 .3 )  ) 

101  FORMAT  < 1H  ,5(15)) 

102  FORMAT  '  1  X  *  34HPRESSURE  EXCEEDS  CRITICAL  PRESSURE) 

103  FORMAT  (5X,4HQRAT* 14X.2HQ1 * 15X.5HGFL0W) 

104  FORMAT  ( 1X*F12®8*5X*F12.8*5X,F12.B) 

105  FORMAT  (5(E12.3) ) 

106  FORMAT  < 1H1 ,6HOUTPUT) 

107  FORMAT  (1H1,34HH.  C.  ROLAND  *  LIQUID  FILM  COOLING) 

I  03  FORMAT  < 1  HO ,  1  OH  I NPUT  DATA) 

109  FORMAT  (5(15)) 

110  FORMAT  (3(E12.3 ) ,3( 14 ) ) 

111  FORMAT  (IX*  40H 

112  FORMAT  (  IX, 3E1 2*3*315 ) 

113  FORMAT  (15HFTLM  EVAPORATED) 

114  FORMAT  ( 3 (FI  0 .5 ) ) 

115  FORMAT  (1X*3(F10.5)> 

116  FORMAT  (7X, 1HS* 1 2X « 3HQ/A, 9X . 2HTB * 5X . 5HNB0 I L > 

117  FORMAT  MX,25HMAX.  DISTANCE  IS  EXCEEDED) 

118  FORMAT  (7X* 1HH* 12X,4HEXP9. 5X.7HTF2 (N1 ) .4X.2HN I ) 
DIMENSION  TF2 ( 21  ) *  V (21  ) * FMU ( 2 1  ) *RHF(2l  ) . TF ( 21  > . CF < 21  ) * 

1 TL (21  ) » CPF (21  ) *RT(21  ) *STAB(40 ) *DTAB(40) «WTA6(40 ) * 

2RL (21  ) *QX (21  ) «  PGT ( 40 ) .TSATT (40 ) *HFTC40) *HGT (40  >  * 

3HWTAB (40 ) »CTl (21  ) *CT2 (21  ) *  CT3 (2 1  ) *CT4  <21  ) *CT5(21  ) 

I  S=0. 000000000 1 

READ  111 

READ  105*  FI,  F2,  F3*  F4 ,  F5 
READ  105*  El*  E2»  E3*  E4 *  E5 
READ  105*  CM *CN2,CN3.CN4*CN5 
READ  105*  G 1 ,  G2,  G3,  G4 ,  G5 
READ  OL-  M*  R2*  R3*  R4  ♦  R5 
REAO  105*  VK1 *VK2* VK3*VK4»VK5 
READ  105*  AK1 *AK2»AK3,AK4*AK5 
READ  105*  PR1 *PR2»PR3*PR4*PR5 
READ  105*  PZ*  TZ*  DT*  Q1  »  TCI 
READ  105*  GK,  RHZ  *  B*  RG  *  PPZ 
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READ  1  05  <>  WCON ,  WTMOL «  TCOM  *  TBN,  DIF 
READ  1 05 «  DX,  SMI,  QDIV,  RGC 
READ  109,  M2,  KSt  N,  NTYPE ,  KPG 
DO  150  J=1,KS 

150  READ  1  05  ,  STAS(J),  DT-AB  [  J  )  .  WTAQ(J),  HWTAB  (  J  5 
DO  200  J=1,KPG 

200  READ  1 05,  PGT < J > ,  TSATT(J>,  HFT(U),  HGT(J) 
PUNCH  107 
PUNCH  1 1 1 
PUNCH  108 

PUNCH  100.F1 ,F2,F3,F4.F5 

PUNCH  100*  El»  E2 ,  E3,  E4 .  E5 

PUNCH  1 00  ,  CN 1  ♦ CN2 , CN3  4  CN4 , CN5 

PUNCH  100,  G1 *  G2.  G3 ,  G4 ,  G5 

PUNCH  1 00, R1 ,R2 ,R3 ,R4 ,R5 

PUNCH  ICO,  VK1 , VK2 , VK3,VK4.VK5 

PUNCH  100.  AK1 , AK2, AK3.AK4, AK5 

PUNCH  lOOoPRl ,PR2,PR3,PR4,PR5 

PUNCH  100»  PZ,  TZ,  DT,  QI  ,  TCI 

PUNCH  100.  GK,  RHZ,  B ,  RG  ,  PPZ 

PUNCH  lOO.WCON,  WTMOL,  TCOM,  TBN »  DIF 

PUNCH  ICO,  DX,  SMI,  QDIV.  RGC 

PUNCH  101,  M2,  KS,  N,  NTYPE.  KPG 

DO  230  J=1 , KS 

230  PUNCH  100.  STAB(J),  DTAB(J>,  WTAB<J>,  HWTAB < J > 
DO  195  J=1»KPG 

195  PUNCH  100,  PGT(J),  TSATT(J>»  HFT  <  J } ♦  HGT ( J > 

169  EXP2=GK/(GK-I • ) 

EXP3=1 ./(GK-1 • ) 

EXP4=2./(GK+1 . ) 

EXP5= (GK— l • (GK+1 • > 

EXP6=.5/EXP5 
EXP7=2.*EXP5 
EXP9=1  • 

DXF=DX/12. 

FN  =  N 

M=N+  1 

QL=OI/FN 

QL2=QL*QDI V 

TSURF=TCI 

NBOIL=0 

NI  =1 

167  DO  240  J=1,M 
TF< J)=TCI 
TF2  (  J  )  =TC  1 
QX(J)=QL 
RL  < J ) -0 • 

CT 1 ( J ) =0 • 

CT2 ( J)=0. 

240  V(J>=0* 
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S=S+DX/2 • 

QTOT=0 • 

DO  500  I  si »M2 
IF  ( S— S7 AB { KS  >  )  43  ♦ 9 ,  550 
9  D=DTABCKS> 

WTHCK=WT  AB  { KS >/12« 

HU/  =  HWT AH  (  KS  ) 

60  TO  11 

43  DO  6  J=1 »KS 

IF  CS-5TAB<jn  7 «  B  ,  6 

7  D»DTA6< J-l  )+< S-STAB CJ-1  ) > * C OTAB ( J ) -OTAB ( J- I  )  )/ 

1 (STAS (J  5-STAB (J-l  )  ) 

WaWTAB< J-l  )  +  <  S-STAB  t  J-l  > ) * C WTAB ( J ) -WTAB ( J- 1  ) )/ 

1 ( STAB  <  J ) -ST AB ( J- 1  )  1 
WTHCK»W/I2« 

HW  =  HWT AS ( J— I  )  +  ( S-STAB (J- 1  ) ) * ( HWTAB ( J > -HWTAB t J-l  )  )/ 
1 CSTAB< J1-STA0< J-l ) 1 
GO  TO  11 

8  D=DT  AB ( J ) 

WTHCK=WTAB< J )/l  2* 

HW  =  HWTAB( J  > 

GO  TO  11 
6  CONTINUE 

1 1  DRATaO/OT 
D=D*. 083333 
VRAT1*1 • 

IF  CS-SM1 )  12, 12* 13 

12  VM=( <EXP4+EXP5*< VRAT1**2) ) **EXP6 >/ ( DRAT**2 > 

GO  TO  14 

13  VM»SQRT( ( <VRAT1*<DRAT**2 J ) **EXP7-EXP4 ) /EXP5 ) 

14  VCRIT«.002#VM 
VCOMoABS (VM-VRATl } 

IF  (VCOM-VCRITJ  16*16,15 

15  VRATl=VM 

IF  (S-SM1 >  12*12*13 

16  PAR=l*+< (GK-1 • )/2. )*{VM**2) 

TG=  <  TZ+460  « ) /PAR 

TGS=TG— 460  « 

RHGoRHZ/'  tPAR»*EXP3  ) 

VG=VM*<  <GK*32*2*RG*TG  ) ■*♦•5 > 

IF  (S— SMI >  30*30*31 

30  TGM= (TG+TSURF+460* ) /2 .+ ,22* .445*TG* ( GK- 1 • 1*(VM**2 ) 
TGMaTGM— 460 • 

TGaTG* ( 1 •  +  *445*  <  GK— 1 • )*<VM**2  )  ) 

TG»TG— 460 • 

GO  TO  32 

31  TGM= ( TG+TSURF+460 • )/2 .+#22**430*TG* < GK-l « )*(VM**2) 
TGM*TGM— 460 • 

TG=TG*C 1 «  +  *430*  <GK-l • )*<VM**2 ) ) 

TG*TG— 460  * 
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32  PG=PZ/ ( PAR**EXP2 ) 

26  IF  (PG-PGT (KPG ) )  84,499,499 
84  CONTINUE 

DO  19  J= 1 , KPG 
IF  (PG-PGTfJM  18,17,19 

17  TSAT=TSATT(J) 

HFS=HFT { J) 

HG=HGT ( J) 

GO  TO  20 

18  PRAT= I PG-PGT { j-i  )  )/<PGT( J)-PGT( J-l  )  ) 

TSAT^TSATT  f  J- 1  1+PRAT* ( TSATT ( J l-TSATT < J- 1  ) ) 
HFS-HFT ( J— 1 J+PRAT+ <HFT ( J ) — HFT ( J-l ) ) 

HG=HGT ( J— l )+PRAT*(HGT( J)-HGT( J-l ) ) 

GO  TO  20 

19  CONTINUE 

20  CONTINUE 

DO  680  J=1,KPG 

IF  (TF2INI > -TSATT < J) >  681,662,680 

682  PPRES  =  PGT(J) 

HFGP=HGT  <  J ) —HFT { J 5 
GO  TO  683 

681  TRAT= ( TF2 (N 1  )-TSATT( J-l  } }/( TSATT (J  >-TSATT( J-l  )  ) 
HFSP  =  HFT( J-l  ] +TRAT  * ( HFT  <  J )— HFT ( J— 1  )  ) 

HGP=HGT< J-l >+TRAT*(HGT( J)-HGT( J-l ) ) 

HFGP=HGP-HFSP 
GO  TO  683 
680  CONTINUE 

683  CONTINUE 

AK= AK1 +AK2*TGM+AK3*TGM**2+AK4*TGM**3+AK5*TGM**4 
GMu  =  Gl+G2*TGM+G3*TGM**2+G4*TGM**3+G5*TGM**4 
PRG=PR  1  -i-PR2*TGM+PR3*TGM**2+PR4*TGM**3+PR5*TGM**4 
VK=VK1  +VK2*TGM-rVK3*TGM*+2+VK4*TGM**3+VK5*TGM*+4 
CON=AK 

RHV=144<,*PG#WTM0L/C  1  545.*  (  TSAT+460  •  >  > 

REG= (RHG*VG*D l/GMU 
DA=3. 142*D*DXF 

DMF  =  DIF*<  <  ITGM+460. )/530. >**1 .5>*< 14.7/PG1 
SCN=GMU/(RHG*DMF ) 

IF  U~1  )  21,21,22 

21  GFLOW=  RHG*VG*.785* (D*+2 ) 

QRAT=OI XGFLOW 

PUNCH  106 
PUNCH  1 03 

PUNCH  104,  QRAT,  01,  GFLOW 

22  TW1  =.0265*<RHG/32.21* (VG**2 )/ (REG**. 2  > 

H- .0265*CON*(REG**o8  >*(PRG**.333  >*EXP9/0 
RHST=QTOT/ ( » 785*D*#2*VG > 

PPST*RHST*RGC*(TG+460. )/l 44. 

IF  (NBOIL)  690,690,25 
690  QEV5  =  H*<TG-TF2(NI  )  ) / ( HFGP*3600 •  ) 
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QEV=QEV5 

HMASS=H*OMF/CON 

RH0S=QEV5/HMASS 

PPRES=PPST+RH0S*RGC*TF2 ( N 1  3/144* 

DO  940  Jel.KPG 

IF  CPPRES  -PGT(J))  941,941*940 

941  PRAT* (PPRES— PGT  < J-l ) 3 / < PGT ( J ) -PGT ( J- 1  ) ) 

TF2CNI  )*TSATTtJ-l 3+PRAT* ( TSATT ( J  )-TSATT < J- 1  ) 3 
GO  TO  942 

940  CONTINUE 

942  QTOT*QTOT+QEV*DA 

PRINT  105,  PPST,  RHST ,  QTOT,  PPRES 
QX  <  N I  )  *QX  C  N I  ) — OEV*DA 
IF  ( QX  C N I  ) )  990  *  990  » 25 
990  MI*NI+1 

OX (MI  ) *QX (MI  )  +QX INI  ) 

NI *N[+1 
25  NCT* 1 
NCT3*0 
141  NCT2«0 
TW*TWl 
LX*M— N I 
DO  62  J*sl  , LX 
L=M-J 

TMsCTFCl  3+TF2  (L  )  )/2* 

C2  =  TM 

C3*TM**2 

C4=C2*C3 

C5=C3*C3 

C6=C3*C4 

FMU  <  L  )  * F 1  /C3+F2/C2+F3+F4*C2+F5*C3 

FMU ( L 3 =  FMU(L  3/32.2 

RHF C  L ) =R 1 +R2*C2+R3*C3+R4*C4+R5*C5 

CF  (L  )  *CN  1  +CN2*C2+CN3*C3+CN4*C4+CN5*C5 

CPF  <L ) =El+E2*C2+E3#C3+E4*C4+E5*C5 

AVa 1 .57l*D*RHF(L)*TW/FMU(L) 

BV*3.142*D*V<L+1  )*RHF  <L ) 

CV*-OX  CL ) 

TL  <  L  3  * (-BV+SQRTI (BV**2 3-4 ,*AV*CV > )/<2«*AV 3 
V  t  L  3  = V  <  L+ 1  3 +TL ( L ) *TW/FMU ( L 3 
IF  (L-N )  63,64,64 

64  RL  <  L  3 *TL ( N 3 / ( 2 • *CF  <N  3  *DA  3  +WTHCK/ < WCON*OA >+!*/<  HW*DA ) 

GO  TO  62 

63  RLCL3«TL<L)/C2**CF<L3#DA3+TL<L+1 3/<2**CP<L+l )*0A) 

62  CONTINUE 
82  NCT* 1 

DO  80  J*N I , N 
J1 *J-1 
J2=U+1 

IF  (NBOILl  405,405,407 
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407  TF2  (NI  ) =  TSAT 
GO  TO  00 

405  IF  tJ-N] )  96*96.70 
96  TF2(NI)=TSET 
GO  TO  80 

70  IF  (NCT2-1  )  700*701*701 

700  CT1  < J)  =  l ./(QX< J >*7200.*CPF ( J)*RL(J1  >  ) 

CT2 ( J )  =  1 •/ < QX  <  J ) *7200 .*CPF ( J  >  *RL <  J }  ) 

CT3  (  J  )  =  1 .+CT1  (  J ) +CT2 ( J ) 

CT 4 ( J ) — CT 1  ( J  )/CT3  <  J 1 
CT5( J)=CT2( J )/CT3< J  > 

701  TFC0M=TF2(J) 

TF2< J1=TF{ J)/CT3(J)+CT4< J)*(TF2( J1  )+TF(Jl  1-TFlJ)  ) 
1+CT5( J )*(TF2 ( J2 )+TF( J2  > -TF { J)  ) 

74  TCT=A8S(TF2< Jl-TFCOM) 

IF  (TCT-TCOM)  80*80*76 
76  NCT -NCT  + 1 

80  CONTINUE 
NCT2=NCT2+1 

IF  <NCT2— 1 0  >  89*141,141 
89  IF  (NCT-1 )  81.81 *82 

81  IF  (NCT3 )  87,87*92 
87  NCT3=1 

GO  TO  141 

92  IF  (NBOIL)  408*408,409 

409  QST  =  o  •' 

DO  410  J  =  NI  *  N 

410  QST=QST+QX( J)*CPF(J)*<TF2(J)-TF(J) ) 
QWT=(TF(N>+TF2<N)-2**TCI ) / < 7200 . *RL ( N > ) 

TW=TW1*EXP9 

UST= (TW*32.2/RHV)**.5 

HM I = • 0265*DMF* ( REG** • 8  >  * ( SCN**.  333 ) *EXP9/D 
QE VI  = ( ( K* ( TG-TSAT ) /3600 • ) -OST/DA-QWT/DA )/ < HG-HFS  1 
422  GMAS=(PG*l44.-QEVl*RGC*(TSAT+460. )/HMl )/(RG* 

1 ( TSAT+460 • ) > 

EXP9=EXP(-13.89*QEV1/<B*UST*<QEV1/HM1+GMAS) > ) 

HM=  «0265*DMF* { REG** *8 )* (SCN***333 )*EXP9/D 
IF  (  ABS (HM-HM1 )-.02*HM)  420,420*421 
421  HM1=HM 

TW=TW1 *EXP9 
UST *  < TW*32 • 2/RHV ) *♦  .5 
GO  TO  422 
420  CONTINUE 

436  HI  =  .0265*CON* ( REG** *8 ) *  < PRG** .333  >*EXP9/D 

QEV=(  (Hl*(TG-TSAT)/3600  •  ) -QST/DA-QWT/DA  >/< HG-HFS  )  - 
GMA=GEV/HM 

PV=GMA* 1 545 .* (TSAT+460. )/ ( 1 44.*WTM0L ) 

IF  (PV-PG)  430*431,431 
431  PV=PG 

GV=PV* 1 44  •/ t  RGC* (TSAT+460 .  5  ) 
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GVM=GV/2. 

PVM=GVM*RGC* (TGM+460. )/l44. 

PSM=PG-PVM 

GSM=PSM# 1 44 ./ ( RG* < TGM+460  * ) ) 

GVM=GVM*RGC/1 545. 

GSM=GSM+RG/ 1 545  * 

XC=GVM/ (GVM+GSM  > 

XS«GSM/ ( GVM+GSM ) 

CON=XC*VK+XS*AK 

H= .0265*C0N* (REG**. 8 ) * < PRG** • 333 ) *EXP9/D 
QEV»  <  (H  *  <  TG-TS AT ) /3600 • }-QST/DA-QtfT/DA  >/(HG-HFS  > 
GO  TO  432 
430  GVM  =  GMA/2  » 

PVM=GVM*RGC* (TGM+460 . )/144. 

PSM=PG-PVM 

GSMoPSM*l 44./ (RG*  f TGM  +  460* ) ) 

GVM=GVM*RGC/1545. 

GSM=GSM#RG/I 545. 

XCsGVM/1  <  GVM+GSM  > 

XS  =  GSM/ C  GVM+GSM ) 

C0N2=XC*VK+XS*AK 
H  =  H1*CON2/-CON 

IF  < A8S(H-HL >-.02*H)  432,432.434 
434  CON=CON2 
GO  TO  436 

432  IF  <ABS(QEV-QEV1  )-.02*QEV>  495*495.494 

494  QEVlsQEV 
GO  TO  422 

495  QT=QEV*C3.14*D*DXF) 

IF  (QT-GL2)  440,440*441 
441  S3S— ( L .-QL2/QT  3  *DX 
DX=QL2*DX/QT 
DXFsDX/12. 

DO  443  J=NI ,N 

443  TF2(J)=TF(J>+(QL2/'QT)*(TF2f  JI-TFC  J)  ) 

QT®QL2 

440  IF  <QT-C5X(Nm  446,447.448 

446  0X(NJ}=QX<N1 )-OT 
GO  TO  300 

447  QX(N!)«C« 

NI -NI+1 

IF  (NI-N)  300,300,575 

448  QEX*»QT-QX(NI  > 

QX(NI )»0. 

N I  "N I  + 1 

QX ( NI  >oQX  <NI  J-QEX 
IF  (Nl-N)  300,300,575 
408  IF  (TF2(NI)-TSAT )  300,302,302 
302  RAT= (TSAT-TF (NI  ) ) / ( TF  2  ( N I  >-TF(NI  )  ) 

N0O I La 1 
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00  305  J«NJ«N 

305  TF2 ( J) =TF( J)+«AT*(TF2 ( J»-TF  C J) > 
S«S-DX+RAT*0X 
300  TMs(TF2(N)+TF(N> >/2* 

QA«(TM-TCI  J/MRU  CN)*DA  ) 

TB-TCI+OA/HW 

PUNCH  116 

PUNCH  112*  S'  QA<  T0 ♦  NBOIL 
PUNCH  118 

PUNCH  l 12«  H*  EXP9«TF2<NI )«NI 
S»S+DX 

450  DO  455  J«NI .N 

TD-IF=TF2CJ)-TF(J) 

TF(J)*TF2<J) 

455  TF2 ( J  >  *TF ( J ) +T0 1 F 
500  CONTINUE 
GO  TO  1 

499  PUNCH  1 02 
GO  TO  1 

575  PUNCH  1  1 3 
GO  TO  1 

550  PUNCH  117 
GO  TO  1 
ENO 
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APPENDIX  C 

GAS  FILM  COOLING  COMPUTER  PROGRAM 


C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


FIRST  INPUT  CARD  4C  BLANKS  FOR  CASE  TITLE. 

SECOND  DATA  INPUT  CARD 

PZ  =  INLET  GAS  PRESSURE.  LB./SQ.iN. 

TZ  *  INLET  GAS  TEMPERATURE.  DEGREES  F. 

DT  *  THROAT  INSIDE  DIAMETER.  INCHES* 

N  *  NUMBER  OF  COOLANT  LAYERS 

M  a  INITIAL  NUMBER  OF  STREAM  LAYERS 

NSH  a  NO.  TO  SPECIFY  TYPE  WALL  SHEAR  CALCULATION. 

0  s  WALL  SHEAR  FROM  BLASIUS  EQUATION. 

1  o  WALL  SHEAR  FROM  PREVIOUS  INCREMENT  PROFILE. 
THIRD  OATA  INPUT  CARD 

GK  =  GAS  CONSTANT.  CP/CV 

RHZ  »  INITIAL  GAS  DENSITY.  LBM./CU.FT. 

QRAT1  =  INJECTION  RATIO  AT  ENTRANCE 
QRAT2  a  SECOND  INJECTION  RATIO 
QRAT3  =  THIRD  INJECTION  RATIO 
FOURTH  DATA  INPUT  CARD 
RG  »  GAS  CONSTANT  R 

WCON  a  LINER  THERMAL  CONDUCT  I  V  I TY .B/HR.FT .DEG .F . 

PR  =  PRANDTL  NO.  OF  COOLANT  GAS. 

M2  =  MAX.  NO.  OF  INCREMENTS 

KS  =  NO.  OF  TABLE  ENTRIES  FOR  D  AND  WTHCK  FROM  S. 
NT  *  NO.  TO  DESIGNATE  TYPE  WALL  COOLING 
0  =  BACKSIDE  COOLED 
1  =  ADIABATIC  WALL 
FIFTH  DATA  INPUT  CARD 

DX  *  DISTANCE  INCREMENT.  INCHES 

EHM  =  RATIO  OF  THERMAL  TO  MOMENTUM  D I FFUS I  V I T I ES 
SMI  a  DISTANCE  TO  NOZZLE  THROAT.  INCHES 
NTYPE  =  CASE  TYPE 

0  a  SINGLE  INJECTION 

1  =  MULTIPLE  INJECTION 

2  a  INJECTION  AT  EACH  INCREMENT 

NU  =  INCREMENT  NO.  FOR  SECOND  INJECTION 
N I 2  a  INCREMENT  NUMBER  FOR  THIRD  INJECTION 
SIXTH  DATA  INPUT  CARD 

RGC  a  COOLANT  GAS  CONSTANT  R. 

GAP  =  COOLANT  INPUT  SLOT  WIDTH,  INCHES 
VCR  IT  =  CONVERGENCE  CRITERION  ON  ITERATIONS 
QINJ  a  INJECTION  RATIO  AT  ENTRANCE 
VI NJ  a  INJECTION  VELOCITY 
SEVENTH  DATA  INPUT  CARO 

TCI  a  INITIAL  COOLANT  TEMPERATURE,  DEGREES  Ff 
DXPT  =  DISTANCE  BETWEEN  PRINTOUTS,  INCHES 
EIGHTH  OATA  INPUT  CARD 
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C  AKI - AK5  COEFFS  FOR  A<  (MOLECULAR ) ,B/HR. FT. DEG. F. 

C  NINTH  DATA  INPUT  CARD 

C  Gl - G5  COOLANT  GMU  COEFFS < MOLECULAR ), LSM/HR .FT • 

C  TENTH  DATA  INPUT  CARD 

C  VI - V5  STREAM  GMU  COEFFS (MOLECULAR > .LBM/HR.FT. 

C  ELEVENTH  DATA  INPUT  CARD 

C  CSl - CS5  COEFFS  FOR  CPS*  B/LBM.DEG.F. 

C  TWELFTH  DATA  INPUT  CARD 

C  CC1 - CC5  COEFFS  FOR  CPC*  3/LBM.OEG.F. 

C  THIRTEENTH  DATA  INPUT  CARD 

C  NPT  «  CONTROL  NUMBER  FOR  PRINTOUT 

C  0  =  VELOCITY  AND  TEMPERATURE  PROFILES  PRINTED 

C  1  =  NO  PROFILE  PRINTOUT 

C  NEXT  KS  INPUT  CARDS- 

C  S — IN..  D — IN.,  WTHCK — IN.*  HW — B./HR.SQ.FT .F. 

C  NEXT  N  INPUT  CARDS 

C  VC1 - FT. /SEC.  TCI - DEG.-  F. 

C  NEXT  M  INPUT  CARDS 

C  VS1 - FT. /SEC.  TS1 - DEG.  F. 

C  OUTPUT  DESCRIPTION 

C  QINJ  a  COOLANT  1NUECTION  RATE.LBM/SEC. 

C  GFLOW  *  MAINSTREAM  FLOW  RATE*  LBM/SEC. 

C  S  =  DISTANCE  ALONG  LINER.  INCHES* 

C  TSURF  =  LINER  SURFACE  TEMPERATURE*  DEG.F. 

C  O/A  =  HEAT  TRANSFER  RATE  TO  LINER.  B/HR.SQIN. 

C  DISP.THCK.  ■  DISPLACEMENT  THICKNESS  AT  S*  INCHES. 

DIMENSION  STAB (60 )  , DTAB (60  > ♦ WTAB (60  >  * TC 1 ( 60 ) » TS 1  ( 20 ) * 
1TC2<60)*TS2(20) ,DC(60  > , DS (20 ) *QS (20 > *QC(60) *VS1  (20>* 
2VS2  <20  }  * VCl  (60 ) * VC2 ( 60 ) »RHS ( 20  >  *RHC(60 ) .EDS (20  )  * 

3EDC (60 ) *DIC<  60 ) *OIS (20 5  .CPC (60  >  * CPS (20 ) . AC  <  60 ) * BC (60 )  * 
4 AS (20) *BS<20  >  *EC (60 ) *ES (20 ) *FC (60 )  ,FS (20) *PC(60 ) * 
3PSI20 ) ,DYC (60 ) ,DYS(20 ) * YC ( 60 ) . YS  < 20  >  *HWTAB(60  > 

100  FORMAT  (5X.4HQINJ* 16X.5HGFL0W) 

101  FORMAT  (3X* 1HS»9X*5HTSURF,9X*3HQ/A,6X, 10HDISP.THCK. ) 

102  FORMAT  (S(E12.3)) 

103  FORMAT  ( 1  HO  *  1  OH  I NPUT  DATA) 

104  FORMAT  (3E12. 3*315) 

105  FORMAT  (F 12  * 5 * BX *F 1 2 .5 > 

106  FORMAT  ( 1  Hi  * 6H0UTPUT ) 

107  FORMAT  ( 1 H 1  *  3 1 HH .  C.  ROLANO  *  GAS  FILM  COOLING) 

108  FORMAT  ( 1 X  » 1 9HAD I  ABA  T I C  WALL  CASE  5 

109  FORMAT  (IX*  20HBACKS I DE  COOLED  CASE) 

I  10  FORMAT ( 1X.40H 

111  FORMAT  (1H  .5 (El 2.3)) 

112  FORMAT  (IX. 3E 12.3,315) 

113  FORMAT  (1X.30HWALL  SHEAR  FROM  PAST  INCREMENT) 

114  FORMAT  (1X.32HWALL  SHEAR  FROM  0LASIUS  EQUATION) 

115  FORMAT  (3F10.5) 

116  FORMAT  (IH  .3F10.5) 
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117  FORMAT  (IX, 2SHM AX •  DISTANCE  IS  EXCEEDED) 

118  FORMAT  UX.16HS INGLE  INJECTION) 

119  FORMAT  < 1 X  *  1 8HMULT I PLE  INJECTION) 

120  FORMAT  ( 1 X , 27H I NJECT I  ON  AT  EACH  INCREMENT) 

621  FORMAT  ( 1X,6(E12.3> ) 

140  FORMAT  (IX*  35HTRANSP 1  RAT  I  ON  WITH  BACKSIDE  COOLING) 
640  FORMAT  (IS) 

647  FORMAT  (IX, 15) 

990  FORMAT  (IX,22HQJ  IS  LESS  THAN  2QC(N)) 

998  FORMAT  (1H  ,IS,5X,E12.3) 

1  S*. 0000000001 

SPT  =0 • 0 
READ  110 

READ  104,  PZ,  TZ,  DT.  N,  M,  NSH 

READ  102,  GK,  RHZ,  QRAT1*  QRAT2 ,  QRAT3 

READ  104,  RG,  WCON,  PR,  M2,  KS,  NT 

READ  104,  DX,  EHM,  SMI,  NTYPE »  Nil,  NI2 

READ  102,  RGC ,  GAP,  VCRIT*  QINJ,  VlNJ 

READ  102,  TCI,  DXPT 

READ  102,  AK1,  AK2,  AK3 ,  AK4 ,  AK5 

READ  102»  Gl,  G2,  G3 ,  G4«  G5 

READ  102,  VI,  V2,  V3«  V4 ,  V5 

READ  102,  CS1 ,  CS2,  CS3,  CS4,  CS5 

READ  102,  CC1 ,  CC2,  CC3 ,  CC4 ,  CCS 

READ  640,  NPT 

PRINT  107 

PUNCH  107 

PRINT  110 

PUNCH  1  1  0 

IF  (NSH)  572,572,573 

572  PRINT  114 
PUNCH  114 
GO  TO  574 

573  PRINT  113 
PUNCH  1 1 3 

574  IF  (NT)  570,570,575 
570  PRINT  109 

PUNCH  109 
GO  TO  580 

575  PRINT  108 
PUNCH  108 

580  IF  (NTYPE- 1 )  820,825,830 
820  PRINT  118 
PUNCH  118 
GO  TO  850 
825  PRINT  119 
PUNCH  1 1 9 
GO  TO  850 

830  IF  (NT)  150,150,151 
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151  PRINT  120 
PUNCH  120 
850  PRINT  103 
PUNCH  1 03 

PRINT  I  1 2 i  PZ,  TZ,  OTi  N,  M,  NSH 
PRINT  llli  GK «  RHZ,  QRAT1 ,  QRAT2  *  QRAT3 
PRINT  1 1 2 «  RG.  WCON, -PR,  M2»  KS.  NT 
PRINT  112#  DX4  EHM,  SMI,  NTYPE 4  Nil,  NI2 
PRINT  Ul»  RGC,  GAP  4  VCRIT,  QINJ,  VINJ 
PRINT  1 11 .  TCI ,  OXPT 
PRINT  111  4  AK1  4  AK2 4  AK3,  AK4  ,  A<5 
PRINT  111  4  Gl,  G2  4  G3  *  G4,  G5 

PRINT  111  4  VI 4  V2  4  V3  4  V4  4  V5 

PRINT  111  4  CSI4  CS24  CS3 4  CS4  4  CS5 

PRINT  111  4  CC1  ,  CC2 4  C C3 4  CCA.  CC5 

PRINT  647  4  NPT 

PUNCH  1124  PZ  4  TZ ♦  DT4  N4  M,  NSH 
PUNCH  111,  GK  4  RHZ,  QRAT 1 4  QRAT2 ,  QRAT3 
PUNCH  112,  RG,  WCON,  PR,  M2,  KS,  NT 
PUNCH  112,  OX,  EHM,  SMI,  NTYPE,  NM<  NI2 
PUNCH  111,  RGC,  GAP,  VCRIT,  Q I N J ,  VINJ 
PUNCH  I  I  I  ,  TCI ,  OXPT 

PUNCH  HI,  AK 1  »  AK2,  AK3  4  AK4,  AK5 
PUNCH  111  4  Gl ,  G2  4  G3  4  G4,  GS 

PUNCH  111  4  VI 4  V2  4  V3  4  V4  4  V5 

PUNCH  111  4  CS14  CS2 1  CS3 4  CS4 4  CS5 

PUNCH  111,  CC14  CC2  4  CCS  4  CC4,  CCS 

PUNCH  6474  NPT 

00  4  J=1 .KS 

2  READ  1024  ST  A8 ( J  1  4  DTABIJ).  WTABIJ),  HWTAB < J ) 
PUNCH  lilt  STAB(J>,  DTAB{J)4  WTABfJ).  HWTAB(J) 

4  PRINT  1114  STAB(J)4  DTABCJ),  WTABIJ),  HWTABCJ) 

DO  904  J=1 ,N 

902  READ  102,  VC1(J>4  TC1<J) 

PUNCH  HI  4  VCI(J>,  TCI(J) 

904  PRINT  til,  VC1  ( J ) »  TCl  <  J) 

DO  905  J= l , M 

903  READ  102,  VSltJJ,  TSI(J) 

PUNCH  111,  VS1(J),  TS 1  <  J ) 

905  PRINT  111,  VS1(J>,  TS 1  I J ) 

GK 1 =GK- I • 

GK2=GK+l • 

EXP2=GK/GK1 
EXP3= 1 • /GK 1 
EXP4«2  4 /GK2 
EXP5=GK1/GK2 
EXP6*,5/EXP5 
EXP7=2.*EXP5 

Cl  NJaCCl  +CC2*TC  I  +CC3*TC  I  **2+CC4*TC  I  **3-fCC5*TC  I  **4 
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TCR I T- VCR  I T 
ML=M-1 

DO  500  I =1 «M2 
N I  =  I 

65  D I =DX 

IF  <  S—ST AS  C <S )  )  5 «  5 ♦ 550 
5  00  6  JMt2 

IF  <S-STAB(KS))  43  ♦ 9 ♦ 550 
9  D=DTAB(<S> 

WTHCK=WTAB  CKS  >*.083333 
HW=HWTAB (KS  > 

GO  TO  1 1 

43  DO  26  K« 1 .KS 
K 1 =K-1 

IF  (S-STAB(K))  17,18.26 

17  0=DT A8  <  K 1  )+(S-STAB<Ki  > ) * ( DTAB tK > -DTA0 (K l  >  >/(STA8(K) 
1  STAB (K 1  J  > 

WTHCK=  (1UTAB(Kl  )+CS-STAB(Kl  >  )  *  <  WTAB  (K  )  — WTAB  CKl  )  >/ 

1 (STAB (K> -STAB <K1 ) ) >*.083333 
HW=HWTAB(K1 ) + < S-STAB < K1 > )*(HWTAB(K )-HWTAB<Kl ) )/ 

1  ( ST  AB ( K  >  —ST AQ ( K 1  )  ) 

GO  TO  11 

18  D=DTAB(K) 

WTHCK=WTAB (K  >*.083333 
HW=HWTA0 (K  > 

GO  TO  11 
26  CONTINUE 
I  1  DRAT=D/DT 

D=D*. 083333 
VRA  T 1  =  1 • 

IF  (S-SM1 )  12,12,13 

12  VM  = l (EXP4+EXP5*(VRAT1**2>  >**EXP6  >/ (DRAT**2 ) 

GO  TO  14 

13  VM=SQRT( ( (VRATl*<DRAT**2>  5 **EXP7-EXP4 )/EXP5 > 

14  VCOMP=.002*VM 
VCOMs ABS  t VM— VRAT1  > 

IF  (VCOM-VCOMP)  16,16,15 

15  VRA  T 1 - VM 

IF  (S-SM1)  12,12,13 

16  PAR=1 .+.5*GK1 *( VM**2 > 

TG= (TZ+460. > /PAR 
TGS=TG 

RHG«RHZ/ (PAR**EXP3 ) 

VG=VM*< (GK*32.2+RG*TG )**.5> 

IF  (S-SM1)  30,30,31 

30  TGM=  (TG+TC1  (NJ+460.  >*.5+.22*.445*TG*G'"l*  Vv**2) 
TG*TG* ( 1 .+.445*GK1 *(VM**2 ) ) 

GO  TO  32 

31  TGM=  <  TG+TC1  ( N ) +460 • >*.5+.22*.430*TG*GK I  * ( VM*+2 > 
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TG=TG*< l #+.430*GKl*(VM**2) ) 
32  PG=P2/ (PAR**£XP2 ) 

IF  <J-1 >  7.7,8 

7  PGl =PG 

TGI =TG«460* 

TGSI=TGS-460. 

TGMI aTGM-460, 

VGl =VG 
RHGleRHG 
01  *D 

WTHCKis»WTHC< 

IF  (NTYPE-1 )  940,943.635 

943  IF  (NI-NI1)  940,941.942 

941  S»S+DI/5* 

SI=S-DI/lO. 

GO  TO  6 

942  IF  (NI»Nll-5)  941.941,944 

944  IF  (NI-NI2)  .635,941.946 

946  IF  (NI-N12-5)  941,941,635 
940  IF  (NI-5)  630.630,635 
630  S=S+Dl/5. 

SlaS-01/10. 

GO  TO  6 
635  S=S+DI 

Sl=S-DI/2. 

GO  TO  6 

8  PG2=PG 
TG2=TG— 460, 

TG52=TGS-460. 

TGM2=TGM— 460  * 

VG2*=VG 

RHG2-RHG 

D2  =  D 

WTHCK2=WTHCX 
6  CONTINUE 

IF  (NI-1J  71,71,72 

71  GFLOW=RHG2*VG2*. 785* < 02**2 > 
QI  =QRAT1*GFI_0W 

PRINT  106 
PRINT  1 00 

PRINT  105,  QI.  GFLCW 
PUNCH  1 06 
PUNCH  1 00 

PUNCH  105,  QI.  GFLOW 
PUNCH  1 0 1 

IF  CNPT)  72,72,672 
672  PRINT  101 

72  RHG*  <RHG 1 +RHG2  >  /2  • 

PG= (PG1+PG25/2. 
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TG* (TGI+TG2 )/2. 

VG= (VG1 +VG2 )/2. 

TGM= ( TGMt  +TGM2 ) /2  * 

PH  I  NIjsPG*!  44  •/■  {  RGC*  (TCI+460*  )  ) 

D*  r0l+02)/2* 

WTHCKo  C  WTHCK 1 +WTHCK2 ) /2 • 

RHGC=PG* 1 44 ./ ( RGC*  C TGM+460 • >  > 

GNUS  =  < V1/TGM+V2+V3*TGM )  / C 3600  * *RHGC ) 

QI «QRAT1*GFL0W 
REG=VG*0/GNUS 
IF  (Nt-1  >  610«610«  605 
605  IF  CNSH)  610«610*615 
6 1 0  TSS=.053*RHG* CVG**2l/(64«4*REG**«2 ) 

GO  TO  620 

615  TSS  =  2«*EDC(N)-*VC2(N)*RHC(N)/(32.2*OC(N)  ) 

620  USTS=SQRT{TSS*32.2/RHGC) 

DO  10  J=1,M 

IF  <NI  -1  ).  35«35«36 

35  RHS ( J ) -PG*1 44 •/ ( RG* tTSl  < J  )+460.  )  ) 

CPS ( J > =CS1 +CS2*TSt ( J ) +CS3*TS 1 t J >**2+CS4*TSl (J)**3+ 
1 CS5*TS 1 (J)**4 
TS2 ( J)=TS1  <  J ) 

VS2 ( J ) =VS 1 ( J ) 

GO  TO  10 

36  TDIF*TS2M»-TS1  (J) 

TS1  C  J ) «TS2 ( J 1 

TS2 ( J ) =TS 1 < Jl+TDIF 

CPS< J >=CS1+CS2*TS1 CJ)+CS3*TS1 ( J >**2+CS4*T51 CJ)**3+ 
1CS5*TS1 (J)**4 

RHS ( J ) =PG*1 44 •/ (RG* ( TS1  (JI  +460*  >  > 

VDIF=VS2CJ)-VS1 <J) 

VS1 < J)=VS2tJ ) 

VS2 t J ) =  VS1  ( J  >+VD IF 
10  CONTINUE 

TS1 (M ) =TGS1 

TS2  <M)=TGS2 

VS2<M)=VG2 

DO  33  J=1 * N 

TF  (NI-l)  81,8t«93 

01  RHC ( J ) =PG*144./ (RGC* ( TC 1  ( J ) +460  * >  ) 

CPC(  J)=CC1+CC2*TC1  <J)+CC3*TC1  t  J  >  **2+CC4*TC  1  C J)**3+ 
1 CC5*TC 1  ( J  >**4 
TC2 ( J)=TC1 ( J) 

VC2 ( J ) *VC1 r J) 

GO  TO  33 

93  TD1  ”=7  )  >— TC  1  ( J  > 

TCI ( J)=TC2 (J) 

TC2{ J)=TC1 ( Jl+TDIF 
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CPC( J )=CC1+CC2+TC1 < J)+CC3*TC1 ( J > **2+CC4*TC  1 <J)**3+ 
lCC5*TCl ( J)**4 

RHC ( J ) =  PG* 1 44 •/ <  RGC*  (TCI  (J)+460* ) ) 

VDIF=VC2(J)-VC1 < J) 

VC1  <  J ) a VC2  <  J ) 

VC2(J)=VC1 ( J )+VDIF 
33  CONTINUE 
NL  =  N— 1 
FNL=NL 

IF  ( N I  —  1  )  1 9, 19*760 
19  VC1 (N) =5**USTS 
DIC(N)=D 

DC ( N  >  = 1 0  **GNUS/USTS 
YC  CN)=DC (N )/2* 

QC  ( N ) =  RHC <  N ) *3* 1 42*D I  C  <  N  5 *DC ( N ) *VC  1  ( N  ) 

IF  (QI~2.*QCCN) )  910,911.911 
911  QL  = (QI -QC (N))/FNL 
DO  50  J= 1 , NL 
50  QC(U)*=QL 

DO  21  J=1 ,M 
21  QS< J)*2,*QC( 1 > 

GO  TO  33 

760  OC (N)=1C.*GNUS/USTS 
YC(N)=DC(N)/2. 

DIC (N)=D 

QCOM=RHC ( N ) *3 • 1 42*D I C  (N ) #DC ( N  > * VC 1  ( N ) 

IF  (NTYPE-l )  20,921,777 
921  IF  (NI-NI1)  20,777,926 

777  QI NJ=GFL0W*QRAT2 
GO  TO  924 

926  IF  CNI-NI2)  20,770,20 

778  QI NJ=GFLOW*QRAT3 

924  IF  (QINJ-QCOM)  983,927,927 
983  QEX=QI NJ+QC (N)— QCOM 
IF  (QEX)  986.987,987 

986  CPS T a  (QC  (N)*CPC  <N).fQINU*C  INJ— QEX*CPC  (NL.  )  >/QCOM 
RHC ( N ) = ( QC ( N  > *RHC  <  N ) +Q I N  J*RM 1 NU-QEX*RHC (NL )  ) /QCOM 
VC1  <N)=(QC(N)*VC1  (N  )+Q  INJ*V  I  NJ— QEX*VCI  (NL)  >/QCOM 
TCI <N)  = <QC(N)*TCl (N)*CPC(N)+QINU*CINJ*TCI- 

1 QEX*TC 1  (NL )*CPC (NL)  )/ (QCOM*CPST ) 

CPC (N) «CPST 
GO  TO  988 

987  QST=QC (NL )+QEX 
CPST»(QEX*CPCCN)+QC<NL)*CPC<NL) ) /QST  ■ 

RHC  <  NL  ) « ( GEX*RHC  ( N  )  +QC  ( NL  )  *RHC  ( NL  )  >  /QST 
VC1  (  NL  )  *  ( Q£X*VC  1  (N  )  +QC  ( NL  > *N/C1  (NL  >  )/QST 

TCI  (NL  )  = (QEX*CPC (N)*TC1  (N >+QC (NL ) *CPC (NL ) *TCl  (NL)  )/ 
1 (QST*CPC (NL ) ) 

CPC(NL)=CPST 
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CPC  CN)= <QtNJ*ClNJ+(QC (N  3-QEX )  *GPC (N>  )/QCOM 

RHC (  N  )  =  {  Q I N J*RH I N J+ (OC(N) -QEX  > *RHC ( N  >  ) /QCOM 

VCI (N ) = ( Q I N J*V 1 N J+ ( QC <N)-QEX)*VC1 (N) l/QCOM 

TCI  (N)3(QINJ*CI1MJ*TCI+{QC(N)-0EX)*CPC(N)*TC1  (N)  )/ 

1 (QCOM*CPC(N) ) 

988  QC<NL>=QC(N1_)+QEX 
QC<N)®QCOM 

IF  ( QC  (  NL  )  /QC  ( 1  )  —2.  •  )  53.861*861 
861  M5=N+1 

TCI <M5)=TCl { N ) 

TC2(M5)*TC2(N> 

VC1 (M5)=VCl (N) 

VC2 (M5)=VC2<N) 

QC CM5)=QC  (N) 

CPC  <M5 )=CPC (N) 

RHC(MS>=RHC(N) 

DIC (M5)=D 
OC(MS)=DCCN) 

QC  <NL)=QC(NLl/2. 

TCI (N ) =TC 1 (NL ) 

TC2 (N)=TC2( NL  > 

VC1 (N)=VCl (NL 1 
VC2(N)=VC2(NL) 

QC ( N ) =QC  f  NL ) 

CPC(N)«CPC(NL) 

RHC (N) «RHC (NL ) 

N=M5 
NL“N-1 
GO  TO  53 

928  QST=5QC(NL)+QC(N)+QEX 

CPST=QC (NL  >*CPC (NL )+QC <N)*CPC IN  >+QEX*ClNJ 
TCI (NL)=(QC(NL>*CPC(NL)*TCl (NL)+QC(N)*CPC(N>*TC1 (N>+ 
1QEX*CINJ*TCI l/CPST 
TC2 (NL  >=TCl  (NL ) 

VC l  < NL  > - (QC ( NL ) *VC 1  ( NL  ) +QC ( N } *VC 1  ( N ) +QEX*V I NJ ) /QST 
VC2 (NL 1 = VCI (NL) 

CPC ( NL ) c  CPST/QST 

RHC (NL )  =  < QC ( NL ) *RHC (NL  >+QC ( N ) *RHC ( N  > +Q£X*RH I N J >/QST 
QC(NL )*QST 

IF  (QC (NL  )/QC ( 1  >-2. )  795.851.851 
851  QC < NL) =QC (NL )/2. 

TCI (N)=TC1 (NL) 

TC2 ( N ) =TC2 ( NL ) 

VCI  ( N ) =  VC  1  (NL) 

VC2 (N) =*VC2 (NL > 

QC(N)=QC(NL) 

CPC (N ) =CPC (NL ) 

RHC (N) =RHC (NL ) 

M5=N+1 
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OIC (M5  )=0 
DC  (  M5 } =DC ( N ) 

N»M5 
NL  =N—  1 

795  QC(N)aQCOM 
TCI C'N)»TC1 
TC2  ( N ) =TC I 
RHC(N)=RH1NJ 
CPC (N) oCINJ 
VC1 (N)=VINJ 
VC2 (N)*VINJ 
GO  TO  53 

927  QEX=QINJ-QC0M 
NXL=QEX/QC { 1 ) 

FNXL»NXL 

IF  (NXL )  928.928.123 
123  QSTaQC CNL>+OC<N) 

CPS T» { QC { NL ) *C PC < NL ) +QC  <N ) *CPC ( N )  > /QST 
RHC ( NL  3  * ( QC ( NL ) *RHC ( NL ) +QC  C  N ) *RHC ( N )  )  /QST 
TCI <NL)=<QC<NL)*CPC(NL)*TC1 t NL ) +QC CN )*CPC ( N ) *TC1 CN))/ 
1 (QST#CPST) 

VCl  ( NL ) ~ (  QC ( NL ) * VC 1  ( NL ) +QC  <  N  >  *VC 1  ( N  >  > /QST 
OC (NL )»QST 
CPC (NL  JeCPST 
M5=N+NXL 
DO  85  J*N,M5 
QC< J)=QEX/FNXL 
RHC ( J)=RHINJ 
TCI (J)*TCI 
TC2  <  J ) aTC I 
VCl (J)=VINJ 
VC2 ( J ) = VI NJ 
CPC< J)«CINJ 
Q5  CONTINUE 

QC  CM5 ) =QCOM 
DC (M5 )*OC (N) 

DIC (M5 )=D 
N=M5 
NL»N-1 
GO  TO  53 

20  QC(NL)=QC(NL>+QC(N)-QCOM 
QC(N)*QCOM 

53  EDC  (N)  =  (G1+G2*TC1  (N)+G3*TC1  ( N  ) **2-M34*TC  1  (N)**3+G5* 

1 TCI  CN)*+4J/(3600**RHC(N>  > 

DO  24  Jsl.NL 

JL*N-J 

JL1 =JL+1 

DIC(JL)=DIC( JLl )—2**DC ( JL I ) 

DC(JL)=QC<JL)/(RHC{JL>*3.I42*DICCJL>*VC1 (JL) ) 
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YC( JL)-YC( JL1  ) +DC ( JL 1  ) /2 • +DC  <  JL ) /2 « 

24  EDCC JU)=.4+USTS*YC(JL) 

010  DO  75  J=1 .M 
J1 =J-1 

IF  (J1 )  22*22*23 

22  DIS(J)*DIC(  1  )-2.#0CU  ) 

DS ( J ) =QS (  J  )  /  ( RHS ( J ) *3 • 1 42*0  I S ( J ) *VS1  (J) ) 

YS ( J  )=YC ( 1  >  +  DC( 1  )/2  »  +  DS ( J >/2. 

GO  TO  75 

23  DISt J)=OIS( J1 )-2.*DS< J1 ) 

DS ( J ) =QS ( J  )  / (RHS ( J ) *3 • 1 42*D I S  <  J ) *VS1  (J)  > 

YS  ( J  )  =YS ( — J l  )+OS ( Jl  1/2.+DS ( J >/2. 

75  EDS ( J )  =  »4*USTS*YS ( J ) 

812  NI T=0 
63  NCT  =  0 

DO  34  J=1 *N 
J1 = J-l 
J2= J+l 
VC0M=VC2< J) 

IF  (NIT)  51 *51 .52 

5 1  PC<J>=3.142*DIC(J)*DC<J>*(PG1-PG2)*72.*<1 .+(D2/Dl  >**2  > 

52  IF  (Jl  )  37.37*38 

37  IF  I N I T )  61.61*62 

61  AC (J >=3.142* (DI/12. ) *D1 S ( J > / ( 32 .2* ( DS < J >/ (RHS ( J ) * 

1  EDS ( J )  ) +OC ( J  > / ( RHC ( J  >  *EDC ( J )  )  >  > 

BC ( J ) =3. 1 42* (Dl/12* )  *D 1 C ( J ) /  ( 32 .2* ( DC ( J 1/ ( RHC ( J  )  * 
lEDC(J)  )+DC(J2)/(RHC(J2)*EDC(J2>  )  )  ) 

62  VC2 ( Jl=(QC(J)*VCl  ( J1/32.2+AC ( J)*(VS1  ( J ) +VS2  <  J >— VC  1  CJ1  ) 
1 +BC ( J )* ( VCl  ( J2>+VC2( J23-VC1  (J) )+PC(J> )/ (QC  <  J  >/32»2 
2+AC ( J ) +BC ( J )  ) 

IF  ( VC2 ( J >  — VG2 )  300.300.301 

301  VC2 ( J ) = VG2 

300  IF  ( A8S ( VC2 ( J  )  — VCQM ) —VCR I T* VC2 ( J )  )  34*34.39 

38  IF  ( J-N>  41*42*42 

41  IF  (NI T ) . 171 ♦ 171 . 172 

171  AC(J)  =  .2610*DI *0  I C ( J 1  ) / ( 32 . 2* ( DC ( J 1 )/(RHC(Jl  )*EDS( J1  )  > 
1 +DC ( J ) / ( RHC ( J  >  *EDC ( J  >  > )  ) 

BC ( J ) =3 • 1 42* (DI/12. ) *D I C ( J ) / ( 32 .2* ( DC ( J  >/ (RHC ( J ) * 

1 EOC ( J >  >+DC ( J2 )/ (RHC ( J2 )*EDC ( J2  )  )  )  ) 

1 72  VC 2 ( J) ~ CQC ( J)*VC1  ( J )/32 .2  +  AC ( J  >* ( VCl  ( Jl  )+VC2 ( J1  >- 
1 VCl  ( J  5  ) +BC <  J ) * ( VCl  ( J 2 )  +  VC2 ( J2  >-VC 1  ( J  )  >+PC ( J  >  )/ 

2<QC( J)/32.2  +  AC( J)+BC( J)  ) 

IF  ( VC2 ( J ) — VG2 )  302*302*303 
303  VC2 ( J ) = VG2 

302  IF  (ABS(VC2( J ) -VCOM ) - VCR I T* VC2 ( J) )  34.34.39 

42  IF  (NIT)  181*181.182 

181  AC ( J >  =  .261 8*0 1*0 1 C ( Jl  )/ (32.2* (DC ( Jl  >/(RHC(Jl )*EDC( Jl  )  > 
1 +DC ( J ) / ( RHC ( J ) *EDC ( J ) 1 ) ) 

BC  (J)=3.142*(DI/12.  )  *0  I C ( J ) / ( 32.2* ( DC ( J ) / ( RHC ( J  >  * 
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1EDC(J) ) > ) 

182  VC2  ( J)=(QC( J)*VC1  ( J > /32 .2  +  AC C J >* < VC  1  ( J1  )+VC2< J1  >- 

1 VC1  ( J  )  ) -BC ( J )*VCl  { J )+PC ( J > 1 / (QC ( J )/32,2  +  AC  < J ) +BC ( J ) ) 

IF  (VC2(J)-VG2)  304*304,305 

305  VC2  C  J  > -  VG2 

304  IF  ( A8S ( VC2 ( J ) —VCOM ) — VCR I T*VC2 ( J ) )  34,34,39 
39  NCT =NCT  + 1 
34  CONTINUE 
816  DO  70  J=1,ML 
J1 = J-l 
J2  °  J+ 1 
VC0M=VS2 ( J  > 

IF  (NIT)  91,91,92 

91  PS( J)=3.142*DIS<J)*DS(J)*(PG2-PG2>*72.*( 1 .+(D2/D1 >**2  ) 

92  IF  (Jl  )  88,88,89 

88  IF  (NIT)  121,121.122 

121  AS( J)=*2618*DI+DIS( J2)/(32*2*(DSCJ2)/<RH5( J2)*EDS(J2) ) 
1 +DS <  J )  /  ( RHS ( J ) *EDS ( J ) ) )  ) 

BS ( J ) =3* 1 42* ( D I / 1 2  *  >  *0 1 S  C  J  > / ( 32 • 2*  <  DS  <  J  > / < RHS ( J  ) * 

1 EDS ( J ) ) +OC ( J ) / ( RHC ( J ) *EDC ( J )  )  )  ) 

122  VS2  CU>»  <QS  ( J  )*VS1  (J)/32»2+AS<  J)*(VS1  <J2  )+VS2(J2  >** 
1VS1  (J)  )+BS(J)*CVCi  <J)+VC2(J)-VS1  (J)  )+PS(J)  )/(QS(  J)/ 

232 • 2+AS  C  J ) +BS { J ) ) 

IF  ( VS2 ( J ) —VG2  >  306*306.307 

307  VS2(J)=VG 2 

306  IF  (ABS(VS2(J)-VC0M)-VCR1T*VS2( J)  )  70,70,69 

89  IF  (NIT)  131,131,132 

131  AS(J)=*2618*D1*DIS(J2)/C32.2*<DS(J2)/(RHS(J2)*EDS(J2) ) 

1  +0S< J)/(RHS( J)*EDSCJ> > > ) 

BS ( J ) *3 • 1 42* ( D J /1 2 • >*D1S  C J)/ (32 .2* (DSC J )/(RHS< J )* 
1EDSCJ) )+DS( J1 )/(RHS(JI >*EDS(J1 ) ) ) ) 

132  VS2 ( J)=(Q5( J)*VS1 l J5/32.2+ASC J)*(VSI (J2  1+V52CJ2  )- 
1VS1  ( J> ) +0S ( J ) * ( VS 1  (Jl  >+VS2< Jl  >-VSl  ( J)  )+PS  t  J  )  )/ 

2(QS C J)/32«2+AS< J)+B5< J) > 

IF  (VS2CJ )-VG2 )  308  #  30B , 309 
309  V52 ( J ) »VG2 

308  IF  (ABSCVS2(J)-VC0M)-VCRIT*VS2C J) )  70,70,69 

69  NCT=NCT+1 

70  CONTINUE 
818  NIT=NCT 

IF  (NCT)  80,80,63 
80  CONTINUE 
N I  T  =  0 
213  NCT  =  0 

00  205  J=1 , N 
Jl =J-1 
J2« J+I 
TC0M»TC2 ( J ) 

IF  (Jl  )  215,215,220 
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215  IF  (NIT)  216*216*217 

216  AC ( J )=«2618*DI *0IS ( J ) / (DS  C J >/ (RHS ( J  >*EHM*EDS ( J )* 

1 CPS { J )  ) +DC ( J  J / < RHC { J ) +EHM*EDC  <  J ) *CPC  <  J ) )  > 

BC ( J ) =3* 142*  <DI /12. )*DIC ( J  >/ (DC ( J )/ ( RHC ( J )+EHM*EDC ( J  > * 
1  CPC ( J  >  >+DC( J2)/(RHC(J21*EHM*EDC( J2)*CPC ( J2 )  )  ) 

EC ( J ) =QC ( J ) *  < VC  1  (J)**2“VC2( J>**2 ) / ( 2 .*778 .*32 .2 ) 

FC ( J  >= ( *2618*01 *DIC (J )*DC ( J ) *RHC ( J ) *EDC ( J )/( 4**778.* 
132  *2  >  )*(  ( VS1  ( J )  +VS2  ( J  )  —VC  1  ( J2)-VC2(J2)  )/(DC<  J>  + 

2DC ( J2)/2.+DS(J)/2. ) >**2 

217  TC2 ( J)=(TC1 (J)*(QC ( J ) *CPC ( J ) -AC ( J ) —BC ( J ) >+AC( J>* 

1  (TS1  ( J  3+TS2  < J ) )+BC ( J )*  < TCI  ( J2 )  +  TC2 ( J2  )  >+EC ( J  >+FC ( J ) ) 

2/ <  QC ( J ) *CPC ( J ) +AC ( J ) +BC ( J ) ) 

IF  ( ABS  <  TC2 ( J ) —TCOM )—TCRIT*TC2(J)  )  205*205*235 

220  IF  CJ-N)  221.222.222 

221  IF  (NIT)  223*223*224 

223  AC ( J )  =  *26 1 8*01 *0 IC ( J 1  )/(DC(Jl  )/(RHC(Jl  > *EHM*EDC ( J 1 ) 
1 *CPC ( J l  )  J+DC ( J )/ (RHC ( J )*EHM*EDC ( J ) *CPC ( J )  )  ) 

BC ( J ) =3 • 1 42* (01/12*  ) *D I C ( J ) / <  DC ( J ) / ( RHC ( J ) *EHM*EDC ( J ) * 
1CPC< J) )+0C ( J2 )/ (RHC ( J2>*EHM*EDC(J2)*CPC(J2> ) ) 

EC ( J )=QC ( J )* (VC1  ( J ) **2-VC2 ( J )**2 )/ (2  **778. *32 *2 ) 

FC ( J ) = ( . 261 8*0 1 *0  I C ( J  >  *DC ( J ) *RHC ( J ) *EDC ( J ) / ( 4 • *778*  * 

1 32  *2) )*( (VC1 (Jl  )+VC2( Jl  )-VCl  (J2)~VC2(J2)  ) / ( DC ( J  ) 

2+DC ( J2 )/2*+DC ( J1 )/2. ) )**2 

224  TC2 ( J ) = (TCI (J)*(QC( J)*CPC< J)-AC( J)-BC(J) )+AC(J)* 

1  (TCI  (Jl  )+TC2( Jl  ) ) +0C ( J ) * ( TC 1  ( J2  J+TC2 ( J2 ) ) +EC ( J ) +FC ( J )  ) 
2/ ( QC ( J  >*CPC ( J ) +AC ( J ) +BC ( J ) ) 

IF  ( ABS ( TC2 ( J )— TCOM ) -TCR I T*TC2 ( J  > )  205.205*235 

222  IF  (NIT)  225.225*226 

225  AC( J)=*26 18*01*01 C(Ji )/(DC( Jl )/(RHC(Jl )*EHM*EDC< Jl > 

1  *CPC ( J 1 ) ) +DC ( J ) *PR/ ( RHC { J ) *EDC ( J ) *CPC ( J ) ) ) 

IF  (NT)  450*450.455 

450  BC ( J )  =  .2618*0  I *0 1 C ( J )/ (DC ( J ) *PR/ (RHC ( J ) *E0C ( J )*CPC( J ) ) 
1 +3600 . *WTHCK/WC0N+3600 • /HW ) 

GO  TO  460 
455  BC ( J ) =0 • 

460  EC(J)=QC(J)*(VC1  ( J)**2~VC2 ( J  )**2 ) / ( 2  .*778 **32 *2  ) 

FC  C J )  =  <  *26 1 3*0 I*DIC(J)*DC(J) *RHC ( J ) *EDC ( J )/ (4**778.* 

1 32 . 2  )  ) * (  ( VC  1  ( J ) +VC2 ( J )  ) / ( DC ( J  3 /2 . J  ) **2 

226  TC2(J)-(TC1  ( J  ) * ( QC < J ) *CPC ( J ) -AC ( J ) -8C ( J ) >+AC(J>* 

1  (TCI  ( Jl  3+TC2 ( Jl  )  )+8C ( J  >* (2  **TCI  ) +EC ( J  ) +FC ( J )  )/ 

2 (QC ( J ) *CPC ( J ) +AC ( J ) +BC ( J ) ) 

IF  ( ABS ( TC2 ( J ) — TCOM ) — TCR I T*TC2 ( J ) )  205.205,235 
235  NCT=NCT+1 
205  CONTINUE 
821  DO  280  Js 1 ♦ ML 
^ i = J- j 
J2=J+1 
TC0M»TS2  < J ) 

IF  ( Jl  )  241 ,241 .242 
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241  IF  (NIT)  243  « 243 1 244 

243  AS( J)=.2618*DI*DIS< J2  )/(DS(J2  )/'(RHS(J2  )*EHM*EDS<  J2  ) 
l*CPS(J2) >+DSCJ)/(RHS(J)*EHM*EDS(J)*CPSCJ) ) ) 

BS ( J ) =3 • 1 42*  ( D I  / 1 2  • )  *  D  I  S  ( J  )  /  ( OS ( J  )  /  ( RHS  <  J ) #EHM*EDS  <  J ) * 
1  CPS ( J ) ) +DC( J)/(RHC ( J >*EHM*EDC ( J  >*CPC <  J  > ) ) 

ES  (  J  )  *CJS  (  J  )  *  ( VSl  ( J  )  **2-VS2  (  J  )  **2  )  /  ( 2  .*778  .*32.2  > 

FS ( J } =  < .2168*01 *DIS< J )*DSC J)*RHS< J)*EDStJ >✓ (4. *778.* 

132.2)  >*( (VS I < J2 )+VS2 ( J2 >-VCl (U)-VC2<J) )/ (OS < J >+0S ( J2 ) 
2/2.+DC ( J )/2. ) )**2 

244  TS2 ( J  )  =  ( TS 1 ( J )* ( QS ( J ) *CPS ( J ) -AS ( J ) -BS ( J ) )+AS( J>* 

1 (TSl (J2)+TS2( J2) )+QS( J)*(TC1 ( J>+TC2( J) )+ES ( J >+FS ( J ) )/ 

2 ( QS ( J ) *CPS ( J ) +  AS ( J ) +BS ( J ) ) 

IF  (A85(TS2( J)-TCOM)-TCRlT*TS2( J) )  280.280.279 

242  IF  (NIT)  247,247.248 

247  AS(J)=.26I8*DI*DrS(J2  )/(DS(J2  }/<RHS(J2  ) *EHM*EDS< J2 ) 
1*CPS(J2) )+DS ( J >/ (RHS( J)*EHM*EDS(J >*CPS( J) ) ) 

BS ( J  )  *3 . 1 42* (Dl/12.)*0IS(J)/<DS(J)/( RHS ( J ) *EHM*EDS ( J ) * 
1CPS(J) )+DS(Jl )/(RHS(Jl ) *EHM*EOS ( Dl )*CPS(J1 ) ) » 

ES ( J )*QS ( J )*  <VSI  ( J)**2-VS2( J)**2)/(2.*77Q.*32.2 ) 

FS ( J ) = ( • 2 1 68*0 1 *D I S ( J ) *DS ( J ) *RH5 ( J ) *EDS ( J )/ (4 .*778 • * 

132.2)  )*( CVS1 (J2)+VS2< J2 >-VSl (J1 )-VS2(Jl ) )/(DS(U>+ 

2DS ( J2 ) /2.+DS ( J 1 )/2» ) )**2 

248  TS2( J)= (TSl ( J)*(QS( J)*CPS( J)-AS(0)-BS( J) )+AS( J)* 

1  (TSl  (J2)+TS2(J2) ) +BS ( J ) * ( TS 1  (Ul  >+TS2(Jl  )  ) +ES ( J )+FS ( J  )  ) 
2/(QS(J ) *CPS  <  J ) +AS  <  J ) +BS ( J  > ) 

IF  (ABS(TS2( J)-TC0M)-TCRrT*TS2( J) )  280.280,279 

279  NCT=NCT+1 

280  CONTINUE 
826  NIT  =  NCT 

IF  (NCT)  290,290,213 
290  CONTINUE 

T A  V= . 5* ( TCI  ( N ) +TC2  <  N )  ) 

IF  (NT)  410.410,420 

41 0  CONI =PR*DC (N )/ (EDC < N ) *CPC < N ) *RHC ( N ) *7200 • ) 
C0N2=WTHCK/WC0N 
C0N3=t  */HW 

OA  = ( T A V— TC 1  >/ (  (CONI +C0N2+C0N3  ) * 1 44 .  ) 

TSURF»TCI+QA* (C0N2+C0N3  >*144. 

GO  TO  426 
420  TSURF«*TAV 
QA  =  0. 

426  NL-N-1 

DYC (N) «6.*DC (N ) 

DO  660  J- 1 , NL 

JL=N-J 

JL1 = JL+1 

660  DYC ( JL ) =DYC ( JUt >+6.*DC( JL l )+6**DC< JL ) 

DO  670  J=l ,M 
JI =  J-l 
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IF  <J1 )  675. 675*680 
675  OVS  < J ) =DYC ( 1  )+6#*DC ( I  )+6.*DS (  1  ) 


GO  TO  670 

680  DYS <  J )  =  DYS ( J 1  ) +6 • *DS ( J 1  ) +6 . *DS ( J ) 

670  CONTINUE 

DSPTK=DYC ( 1 )+6«*0C ( 1 1 
DO  720  J=  1  » M 

720  OSPTK=DSPTK+( I *-RHS< J )*VS2(J )/<RHG2*VG2 > ) *D I S (U  J *1 2«* 
IDS  ( J )/DlS< 1  ) 

IF  (NPT)  641.641.642 

641  IF  CS1-SPT)  500.501.501 

501  PRINT  621 .  TC2 
PRINT  621.  VC2 
PRINT  621 .  QC 
PRINT  101 

642  IF  (SI  — SPT )  500  *502  » 502 

502  PRINT  102.  SI.  TSURF.  QA «  DSPTK 
PUNCH  102.  SI.  TSURF.  QA .  DSPTK 
SPT=SPT+OXPT 

500  CONTINUE 

550  PRINT  117 
GO  TO  995 

910  PRINT  990 
GO  TO  995 

150  PRINT  140 

995  GO  TO  1 
END 
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APPENDIX  D 

TRANSPIRATION  COOLING  COMPUTER  PROGRAM 


C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


FIRST  INPUT  CARD  40  BLANKS  FOR  .CASE  TITLE. 

SECOND  INPUT  CARD 

G 1 - G5  COEFFS  FOR  STREAM  V  I  SC . , LBM/FT • SEC . 

THIRD  INPUT  CARD 

VK1 - VK5  COEFFS  FOR  VAP • COND • ♦ B/HR .FT • F • 

FOURTH  INPUT  CARD 

AKi - AK5  COEFFS  FOR  STREAM  COND • .B/HR .FT .F • 

FIFTH  INPUT  CARD 

PR1 - PR5  COEFFS  FOR  STREAM  PRANDTL  NO. 

SIXTH  INPUT  CARD 

PZ  *  INLET  GAS  PRESSURE  P-ZERD , LB/ I N. SQ • 

TZ  =  INLET  GAS  TEMP . T-ZERO ♦ DEG • F • 

DT  *  THROAT  DIAMETER  D* .  INCHES. 

TSURF  =  SPECIFIED  SURFACE  TEMP 
CPL  =  HEAT  CAPACITY  OF  TRANSP*  L I Q» , S/LBM.F* 
SEVENTH  INPUT  CARD 

GK  =  MAINSTREAM  RATIO  OF  SPECIFIC  HEATS 
RH2  =  INITIAL  MAINSTREAM  DENS I TY ♦ LBM/CU .FT . 

B  *  TURCOTTES  CONSTANT  FOR  WALL  SHEAR  CALC. 

RG  *  STREAM  G£  S  CONSTANT  R. 

CPV  =  HEAT  CAPACITY  OF  TRANSP. GAS  OR  VAP . . B/LBM . F . 
EIGHTH  INPUT  CARD 

WTMOL  =  COOLANT  MOLECULAR  WEIGHT.- 

DIF  =  COOLANT-AIR  DIFFUSION  COEFF • .FT. SO ./SEC . 

DX  =  DISTANCE  INCREMENT.  INCHES. 

SMI  a  DISTANCE  TO  NOZZLE  THROAt.  INCHES* 

NINTH  INPUT  CARD 

M2  =  MAXIMUM  NUMBER  OF  INCREMENTS* 

KS  «  NO.  OF  TABLE  ENTRIES  FOR  D  AND  S. 

KPG  =  NO.  OF  TABLE  ENTRIES  FOR  TSAT.HF.HG. 

NTYPE  *  CASE  TYPE  NUMBER 
0  =  GAS  TRANSPIRATION 
1  =  LIQUID  TRANSPIRATION 

NEXT  KS  CARDS 

ON  EACH  CARD  IN  ORDER  AS  FOLLOWS 
S — INCHES.  D — INCHES.  WTHCK — INCHES. 

NEXT  KPG  CARDS 

ON  EACH  CARD  IN  ORDER  AS  FOLLOWS 

PG — PSIA.  TSAT — DEG.F • ,  HF — BTU/LB,  HG — BTU/LB. 
OUTPUT  DESCRIPTION 

S  =  DISTANCE  ALONG  WALL  FROM  INJECTION.  INCHES. 

QEV  =  INJECTION  RATE  REQU I  RED ,  < LSM/SEC ) /SQ .FT * 

QTOT  =  TOTAL  INJECTION  TO  S.  LBM/SEC • 

QRAT  *  RATIO  OF  QTOT / f  GAS  FLOW  3  TO  S. 

TG  =  MAINSTREAM  STATIC  TEMPERATURE,  DEG.  F. 

PG  *  LOCAL  STATIC  PRESSURE,  PSIA. 
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C  EXP9  =  EXPONECTIAL  ATTENUATION  FACTOR* 

100  FORMAT  (IH  ,5(E12*3>> 

101  FORMAT  ( 1H  ,5(15)) 

102  FORMAT  (IX*  34HPRESSURE  EXCEEDS  CRITICAL  PRESSURE) 

103  FORMAT  ( 1 X  * 20HL I QU I D  COOLANT  TABLE) 

104  FORMAT  { IX, 12HN0Z2LE  TABLE) 

105  FORMAT  <5<E12.3)> 

106  FORMAT  ( 1 H 1 , 6H0UTPUT ) 

107  FORMAT  (  1  HI  ,  39HH*  C*  ROLANO*LIQ*  AND  GAS  TRANSPIRATION) 

108  FORMAT  < 1  HO , 1  OH  I NPUT  DATA) 

109  FORMAT  (5(15)) 

110  FORMAT  (3<E12.3).  3(15)) 

1 1 1  FORMAT  ( IX, 4CH  ) 

112  FORMAT  ( 1X,34HINC0MPAT1BLE  —  TSURF  LESS  I  HAN  TSAT > 

113  FORMAT  ( 7X , 1 HS , 1 1 X , 3HQEV , 9X , 4 HQ TOT , 8X , 4HQRAT  ) 

114  FORMAT  ( 1  HO , 5X , 5HGFL0W ) 

115  FORMAT  (7X.2HTG. 10X.2HPG,9X,4HEXP9) 

117  FORMAT  (1X»25HMAX.  DISTANCE  IS  EXCEEDED) 

DIMENSION  STAB (40) ,DTAB(40 ) , WTAB(40  )  ,PGT(40 ) ,TSATT(40 ) , 
1HFT (40 ) ,HGT (40 ) 

1  S=0*0000000001 

READ  1 t 1 

READ  105,  Gl*  G2,  G3,  G4 ,  G5 

READ  105,  VK1,  VK2.  VK3.  VK4,  VK5 

READ  105,  AK1  ,  AK2,  AK3 ,  AK4 ,  AK5 

READ  105,  PR  1 ,  PR2 ,  PR3 ,  PR4 ,  PR5 

READ  105,  P2, TZ,DT, TSURF ,CPL 
READ  105,  GK,  RHZ,  B.  RG ,  CPV 

READ  105,  WTMOL,  DIF,  DX ,  SMl 

READ  109,  M2,  KS,  KPG  , NTYPE 
DO  200  J=1 ,KS 

200  READ  105,  STA8(J),  DTAB(J),  WTA0 ( J ) 

DO  201  J=1,KPG 

201  READ  105,  PGT(J),  TSATT(J),  HFT(U),  HGT(J) 

PUNCH  107 

PUNCH  1 1 1 
PUNCH  108 

PUNCH  100.  Gl.  G2,  G3,  G4 .  G5 
PUNCH  100,  VK1,  VK2.  VK3 ,  VK4,  VK5 
PUNCH  100,  AK 1 ,  AK2,  AK3 ,  A<4,  AK5 
PUNCH  100*  PRl *  PR2,  PR3,  PR4,  PR5 
PUNCH  100,  PZ,  TZ,  DT.  TSURF, CPL 
PUNCH  100,  GK,  RHZ*  B,  RG  , CPV 
PUNCH  100,  WTMOL,  DIF,  DX .  SMl 
PUNCH  101.  M2r  KS.  KPG  , NTYPE 
PUNCH  104 
DO  202  J=1 , KS 

202  PUNCH  100*  STAB(J),  OTAS(J>»  WTABCJ) 

PUNCH  103 

DO  203  J=1«KPG 
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203  PUNCH  100*  PGT  <  J  )  *  TSATT ( J  )  *  HFT(J).  HGT ( J ) 

PUNCH  106 
EXP2=GK/<GK-1 • ) 

EXP3=1 ./{GK-1 • > 

EXP4=2*/(GK+1 • ) 

EXP5=  1GK-1 • )/( GK+1 • ) 

EXP6=.5/EXP5 
EXP7=2«*EXP5 
EXP9»1 • 

S=S+DX/2. 

DXF=DX/12* 

QT0T=0. 

DO  500  1=1 * M2 
IF  ( S— ST AS (KS 1 )  43*9*550 
9  D*DTAB(KS> 

WTHCK  =  U/TAB  <  KS  > / 1 2 . 

GO  TO  11 

43  DO  6  J= 1 «KS 
J1 = J-l 

IF  (S-STAB<U>)  7*8*6 

7  D«DTAB< J1  )  +  <  S— STAB  f  Ul  )  )  *  ( DTAB  <  J  1  -DTAB  C  J 1  1  )/ 

1 (STAB ( J l-STAB ( J1  1  ) 

W=WTAB< J1  )+ (S— 5TA8  <  Ul  >  >* (WTAB ( J )-WTAB C Jl  } )/ 
l  ( ST AB ( J )-ST AB  <  J 1  1  J 
WTHCK=W/12* 

GO  TO  11 

8  D»DTAB(J> 

WTHCK=WTAB ( J ) / l 2 ♦ 

GO  TO  1 1 

6  CONTINUE 

11  DRAT=D/DT 
D=0*. 08333 
VRA  T  1  =  l • 

IF  (S-SM1 )  12*12*13 

12  VMs( (EXP4+EXP5* ( VRAT1 **Z ) ) **EXP6 ) / f DRAT**2 ) 

GO  TO  14 

13  VM=SQRT( { (VRAT1*<DRAT**2> >**EXP7-ExP4 )/EXP5 1 

14  VCR  I T= • 002*VM 
VCOMaABS(VM-VRATl ) 

I F ( VCOM-VCR I T )  16*16*15 

15  VRAT 1 = VM 

IF  CS-SM1 >  12*12*13 

16  PAR=1  .+  (  (GK-1  •  )/2.  >*<VM-»-»2  > 

TG= (T2+460. ) /PAR 
TGS=TG-460* 

RHG=RH2/ (PAR**EXP3 1 
VG=VM*I  (GK*32«2*RG*TG  )**»S  1 
IF  (S-SM1 >  30*30.31 

30  TGM=  (TG+TSURF+460.  ) /2  .+  .22* .445*TG*  (GK- 1  .  )*  <VM**2) 
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TGM=TGM-460. 

TG=TG*( I .+.445*CGK-t • >*(VM**2> } 

TG=TG-460* 

GO  TO  32 

31  TGM= (TG+TSURF+460. > /2 . + .22* «430*TG* { GK- 1 • 

TGM=  TGM— 460 . 

TG  =  TG*< 1 • +4  30 •  *  (  GK— 1 • )*(VM**2>  ) 

TG=TG— 460« 

32  PG=PZ/(PAR**EXP2> 

26  IF  (PG— PGT (KPG  5  )  84*499.499 
84  CONTINUE 

00  19  J= 1  *  KPG 
IF  (PG-PGT(J)]  18.17*19 
17  TSAT  =  TSATT(J) 

HFS=HFT( J) 

HG*HGT< J) 

GO  TO  20 

10  PRAT= (PG-PGT ( J-l  ) )/lPGT(U)-PGT( J-i  )  ) 

TSAT  =  TSATT ( J-l  M-PRAT* (TSATT ( J  J-TSATT ( J-l  )  ) 
HFS=HFT( J-l )+PRAT*<HFT( JJ-HFTCJ-I ) ) 

HG=HGT< J-l  )+PRAT*(HGT(J )-HGT(J-l  )  ) 

GO  TO  20 

19  CONTINUE 

20  CONTINUE 

IF  (NTYPE)  700.700.701 
701  IF  (TSURF-TSAT 1  720.700*700 

700  AK=AK1 +AK2*TGM+AK3*TGM**2+AK4*TGM**3+AK5*TGM**4 
GMU=G  l  +G2*TGM+G3*TGM**2+G4*TGM**3+G5*TGM**4 
PRG  =  PR  1  +PR2*TGM+PR3*TGM**2+PR4*TGM**3+PR5*TGM-»*4 
VK=VK1 +VK2*TGM+VK3*TGM**2+VK4*TGM**3+VK5*TGM**4 
CONeAK 

0=  (DT*DRAT )/12. 

RHV=144.*PG*WTM0L/( 1 545 . * ( TSURF+460 •  )  ) 

REG“  CRHG*VG*0)//GMU 

DMF=DIF*( ( (TGM+460. >/530. >**1 .5>*( I4.7/PG) 
SCN=GMU/ <  RHG*OMF  > 

IF  ( 1-1  )  21.21.22 

21  GFL0W=RHG*VG*.785* (0**2 > 

PUNCH  114 

PUNCH  100.  GFLOW 

22  TW1 =*0265*(RHG/32.2 >* < VG**2 )/ I REG**. 2 ) 

H=  .0265*C0N*  (REG**. 8  )  *  ( PRG** .  333  )  +EXP9/D 
TW=TW1*EXP9 

UST= (TW*32.2/RHV>**.5 

HM1 =.0265*DMF*(REG**.8 )*(SCN**,333 )*EXP9/D 
IF  (NTYPE)  470.470.471 

470  OEV 1  =  ( H*  < TG-TSURF ) /3600 . ) / ( CPV* ( TSURF-TC I )  ) 

GO  TO  422 

471  QEV1=(H* (TG-TSURF >/3600. ) / ( CPL* ( TSAT-TC I 3+CPV* 
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1 (TSURF-TSAT ) ) 

422  GMAS«*  (  PG*1 44  ♦— QEVl^RGC*  < TSURF+46Q  •  1/HM1  >/(RG* 

1 (TSURF+460* ) ) 

EXP9=EXP (  —  1 3»89*QE Vl / CB*UST  4  C  QEV 1 /HM 1 +GMAS )  >  > 
HMs,0265*DMF*(REG***8 )*<5CN**.333 )*EXP9/D 
IF  (ABSIHM-HMl )-,02*HM)  420*420*421 

421 

TW=TW1*EXP9 
UST  a (TW*32 •2/RHV ) ** *5 
GO  TO  422 
420  CONTINUE 

436  HI  =  *0265*C0N*  <REG*+«8 )* (PRG*+*333 )*EXP9/D 
IF  < NT YPE  >  450*450*451 

450  QEVl=(HI*(TG-TSURF)/3600.>/CCPV*<TSURF-TCt ) ) 

GO  TO  452 

451  QEV1  •*(  HI  *(TG-TSURF  1/3600.  ) / ( CPL* < TSAT-TC I  )+CPV*(TSURF 
1-TSAT)  ) 

452  GMA*sQEV1/HM 

PV=GMA*1545.*(TSURF+460* ) / ( 1 44 .+WTMOL > 

IF  (PV-PGl  430*431*431 
431  PV=PG 

GV=PV*144*»tfTMOL/C1545**(TSAT+460. ) l 
GVM«GV/2. 

PVMaGVM*l'545.*  CTGM+460*  )/ (  1  44 .+WTMOL  ) 

PSM-PG-PVM 

GSM^PSM* 1 44  * / ( RG* ( TGM+460 • >  > 

GVM»GVM/WTMOL 
GSM=GS«*RG/'  1 545  ■ 

XC  »GVM/ ( GVM+GSM ) 

XSsGSM/'  ( GVM+GSM ) 

CON2eXC*VK+XS*AK 
H1=H]*C0N2/'C0N 
GO  TO  433 
430  GVM=*GMA/2. 

PVM®GVM*l 545»* ( TGM+460 • >/( 144**WTM0L ) 

PSMoPG-PVM 

GSM*PSM* 1 44 */<RG* (TGM+460* ) J 

GVM»GVM/WTMOL 

GSM»GSM*RG/ 1 545 • 

XC=GVM/ <  GVM+GSM ) 

XS=GSM/ (GVM+GSM  > 

CON2»XC*VK+XS+AK 
HI =Hl*CON2/CON 
433  IF  (NTYPE )  460*460*461 

460  0EV»(Hl*(TG-TSURF)/3600. >/( CP V* ( TSURF-TC 1 ) ) 

GO  TO  462 

461  QE  V°  <  H 1  * ( TG— TSURF ) /3600 • ) / l CPU* ( TS AT-TC 1  ) +CPV* ( TSURF 
1-TSAT) > 

462  IF  ( ABS  <  QEV— QEV 1  )  —  • 02+QEV  >  465*465*464 
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464  QEV 1 =QEV 
GO  TO  422 

465  CONTINUE 
QTOT=QTOT+QEV*3. 14*D*DXF 
QRAT=QTOT/GFLOW 

PUNCH  113 

PUNCH  lOOi  S<  QEV «  QTOT ♦  QRAT 
PUNCH  115 

PUNCH  100.  TG.  PG#  EXP9 
S=S+OX 

500  CONTINUE 
550  PUNCH  1  17 
GO  TO  1 

720  PUNCH  112 
GO  TO  1 

499  PUNCH  102 
GO  TO  1 
END 
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APPENDIX  E 
PHYSICAL  PROPERTIES 

The  eiscosity,  thermal  conductivity,  and  specific  heat  at 

constant  pressure  for  air  are  given  in  graphical  form  in  Figures  28, 

29,  and  30,  respectively.  The  data  for  temperatures  below  3000  °F. 

were  taken  from  tables  in  Kreith’s  Principles  of  Heat  Transfer.  The 
o 

data  above  3000  F.  were  taken  from  references  (15)  and  (28)  in  the 
Bibliography. 
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Fig-  29  Thermal  Conductivity  of  Air  as  a  Function  of  Temperature 
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APPENDIX  F 

SAMPLE  COMPUTER  CALCULATIONS 

This  Appendix  Includes  Tables  VII,  VIII,  which  are  the 
sample  computer  calculation  outputs  for  transpiration  cooling  using 
the  gas  film  cooling  program,  transpiration  cooling  using  the  liquid 
or  gas  transpiration  program,  and  the  liquid  film  cooling  calculation, 
respectively. 
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INPUT  DATA 
0.441E  03 
0*1 28E  01 
0.533E  02 
0. 100E-0 1 
0.533E  02 
0.100E  03 
0. 170E-01 
-C.262E  01 
-0.262E  01 
■  0.238E  00 
0.238E  00 
1 

0.' 

0.250E-01 
0.S0CE-01 
0.75CE-01 
0.100E  00 
0.125E  00 
0.150E  00 
0*1 75E  OC 
0.200E  00 
0 «  225E  00 
0.550E  00 
0.575E  00 
0  .600E  00 
0 .625E  00 
0.650E  00 
0.675E  00 
0 • 70 OE  00 
0.725E  CO 
C.75CE  00 
0.302E  01 
0 • 1 27E  03 
0.122E  03 
0*1 19E  03 
0. 1 16E  C3 
0.1 15E  03 
0.1 13E  03 
0.  A 1 1 1  03 
0.1 10E  03 
3.108E  03 
0.1C7E  03 


TABLE  VII 

SAMPLE  CALCULATION  FOR  TRANSPIRATION 
USING  GAS  FILM  PROGRAM 


0.300E  04 

0.500E  CO 

20  20 

0 

0.342E  00 

0.100E  00 

0.200E-03 

o. 

0.190E  03 

0.800E  00 

320  20 

1 

0.1SCE  01 

0.700E  OC 

2  320  320 

0.150E-01 
0.100E  00 

0.1  COE— 0 1 

0.100E  00 

C.750E  02 

0 .694E-05 

0.139E-07 

-0.717E-1 1 

0.125E-14 

0.700E-0 1 

0.233E-04 

n  • 

o. 

0 • 700E-0 1 

0.233E-04 

o. 

c. 

0.233E-04 

0. 

0  • 

o. 

0.233E-04 

o. 

o. 

o. 

0.150E  01 

u.ioOE  00 

0  «  50 OE  04 

0.150E  01 

0*1 80E  00 

0.500E  04 

0.1 49E  01 

0.160E  00 

0.500E  04 

0.1 48E  01 

0.140E  00 

0.5C0E  04 

0.1 47E  01 

0 . 1 1 OE  00 

0.500E  04 

0.14SE  01 

0.800E-01 

0.500E  04 

0.1 42E  01 

0  .500E-01 

0  * 500E  04 

0 • L39E  01 

0.450E-01 

0.500E  04 

0.134E  01 

0.400E-01 

0.500E  04 

0.130E  01 

0.350E-01 

0.500E  04 

0.660E  00 

0.300E-01 

0*5C0E  04 

0  s  6 1 OE  00 

0.300E-01 

0 • 500E  04 

O.570E  OC 

0  *300E— 0 1 

0  *50CE  04 

0.540E  00 

0.300E-01 

0  * 500E  04 

0.515E  00 

0  .300E-0 1 

0 • 5oOE  04 

0.505E  CO 

G  .300E-0 1 

0«E00E  04 

0.50CE  ^0 

0.300E-C1 

0 • 500E  04 

0.5C5E  00 

0  .300E-0 1 

0.500E  04 

O.S15E  OC 

0 .300E-C1 

n.500E  04 

0.200E  01 
0.135E  04 
0.1 27E  04 
0.1 18E  04 
C.i 1  IE  04 
O.103E  04 
0.96SE  03 
0.902E  03 
0.844E  03 
0.791E  03 
0.743E  03 

0.140E  CC 

0.500E  04 
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0*1 06E  03 
0.104E  03 
0*1 03E  03 
0.101E  03 
0 • 1 OOE  03 
0.980E  02 
0 .960E  02 
0.940E  02 
0 .930E  02 
0.510E  02 
0.134E  03 
0.148E  03 
0 • 1 62E  03 
0 • 1 71 E  03 
0*1 76E  03 
0*1 78E  03 
0*1 78E  03 
0*1 7SE  03 
0.178E  03 
C.l'7QE  03 
0.178E  03 
0.1 78E  03 
0.178E  03 
C.178E  03 
0.178E  03 
0.178E  03 
0.178E  03 
0.178E  03 
0.1 78E  03 
0.1 78E  03 
OUTPUT 

Q  INJ 

0.07537 

S 

C.500E-C2 
0.105E  00 
0.205E  00 
C.305E  OC 
0.405E  00 
0.505E  00 
0.6C5E  00 
0.705E  00 
0.805E  00 
0.905E  CC 
0.1 OOE  01 
0.1 10E  01 


TABLE  VII  (Continued) 


0.699E 

03 

0.659E 

03 

0.623E 

03 

0.590E 

03 

0.560E 

03 

0.534E 

03 

0.509E 

03 

0.484E 

03 

0.468E 

03 

0.237E 

03 

0.149E 

04 

0.1 68E 

04 

0.1 87E 

04 

0.204E 

04 

0.220E 

04 

0.235E 

04 

0.2488 

04 

0.259E 

04 

0.268E 

04 

0 . 275E 

04 

0 .281 E 

04 

0.285E 

04 

0.289E 

04 

0.291E 

04 

0.294E 

04 

0.296E 

04 

0.298E 

04 

0.299E 

04 

0.300E 

04 

0.300E 

04 

GFLOW 

0.75374 


TSURF 

Q/A 

oisp.thck. 

0.248E 

03 

0. 

0.157E-01 

0.312E 

03 

0. 

0.I44E-01 

C.353E 

03 

0- 

0.139E-01 

0.41  IE 

03 

0. 

0 . 134E-01 

0.484E 

03 

0. 

0. 122E-01 

0.582E 

03 

0. 

0.107E-01 

0.70  IE 

03 

0. 

0.865E-02 

0  *o04E 

03 

o. 

0.740E-02 

0 .800E 

03 

0. 

0.716E-02 

0.80  1  E 

03 

0* 

0.772E-02 

0 . 781 E 

03 

o. 

0.868E-02 

0.76BE 

03 

0. 

0.977E-02 

137 


A.EDC-TR-48-88 


TABLE  VII  (Concluded) 


U  •  1 20 E 

01 

0  a 756E 

03, 

0* 

;o. li2E-oi 

0 • 1 30 E 

01 

0.698E 

03 

-■  o 

oJ»422E-0  1 

0.140E 

01- 

0*6S6E 

03 

0* 

0>1 37E-0 1 

0.150E 

Cl 

0*673E 

03 

0* 

0 • 1 53E-0 1 

0*1 60 E 

01 

0  *603E 

03 

0* 

0*163E-01 

0*1 70E 

01 

0*594E 

03  ' 

o. 

. 0* 180E-01 

0*180E 

01 

0.544E 

03 

0* 

.0*191E-01 

0.190E 

01 

0*515E 

03 

0. 

vb*206E-01 

0*200E 

01 

0.510E 

03 

0*  . 

Q.*226E-0 1 

0.210E 

01 

0  «444E 

03 

0* 

tO .232E-0 1 

0.220E 

01 

0  *440E 

03 

0* 

' 0*25 IE— 01 

0  «230E 

01 

0  *43 1 E 

03 

0*. 

. 0.265E-01 

0.240E 

01 

G.371E 

03 

0* 

;0*275E-0l 

0 e25CE 

01 

0  •  J77E 

03 

0* 

0  *298E— 01 

CU260E 

01 

0*340E 

03 

0. 

o*bosE-oi 

0*270E 

,01 

0  *31 TE 

03 

0*- 

■'0>.320E-0 1 

0.280E 

01 

0  *322E 

03 

0* 

‘0.347E-01 

0.290E 

01 

Q*262E 

03 

0* 

';i:*344E-ol 

0  *3Q.0E 

ot 

.  0.263E 

03 

0* 

,0.*366E-0  I 
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TABLE  VMI 

SAMPLE  CALCULATION  FOR  TRANSPIRATION  USING 
LIQUID  OR  GAS  TRANSPIRATION  PROGRAM 


INPUT  DATA 

— 

0. 100E-04 

0. 100E-07 

O.OOOE-38 

0*000E-38 

0*000£— 38 

0.900E-02 

0.300E-04 

0  .000E-38 

0  *000E— 38 

O.OOOE-38 

0. 133E-01 

0. 170E-04 

O.OOOE-38 

O.OOOE-38 

O.OOOE-38 

0.700E 

00 

O.OOOE-38 

O.OOOE-38 

O.OOOE-38 

O.OOOE-38 

0.450E 

03 

0.2S0E 

04 

0.500E  00 

0.500E 

03 

0.1 OOE  01 

0.  I40E 

01 

0.400E 

00 

0. 1 10E  01 

0.533E 

02 

0.460E  00 

0« 1S0E 

02 

0.  100E 

01 

0.500E-01 

0.700E 

00 

110  10  23  1 

NOZZLE  TABLE 

0*000E-38 

0. 150E 

01 

0.400E-01 

0*1 OOE 

01 

0.590E 

00 

0.400E-01 

0*1  03E 

01 

0 . 550 E 

00 

0.400E-01 

0.1 05E 

01 

0.520E 

00 

0.400E*-01 

0. 10SE 

01 

0.495E 

00 

0.400E-01 

©•■1  1  OE 

01 

0  •  4QSE 

00 

0.400E-01 

0. 1  13E 

01 

0.480E 

00 

0.400E-01 

0 . 1  1  SE 

01 

0.48SE 

00 

0.400E-01 

0.  1 17E 

01 

0.495E 

00 

0.400E— 01 

0*31 3E 

01 

0.200E 

01 

0.160E  00 

LIQUID  COOLANT 

TABLE 

0.2O0E 

00 

0.531 E 

02 

0.212E  02 

0.1 08E 

04 

0.500E 

00 

0.796E 

02 

0.476E  02 

0. 1  1 OE 

04 

0.  I00E 

01 

0.102E 

03 

0.697E  02 

O.  1  1  IE 

04 

0  *200E 

01 

0. I26E 

03 

0.940E  02 

0. 1 12E 

04 

0.400E 

01 

0. 153E 

03 

0.121E  03 

0» 1 13E 

04 

0.700E 

01 

0*1 77E 

03 

0.145E  03 

0. 1 I4E 

04 

0.  I20E 

02 

0.202E 

03 

0.170E  03 

0*1 15E 

04 

0*200E 

02 

0.22BE 

03 

C.196E  03 

0 •  I  1 6E 

04 

0.400E 

02 

0.267E 

03 

0.236E  03 

0.1 17E 

04 

0.700E 

02 

0.303E 

03 

0.273E  03 

0 • 1 I  BE 

04 

0*1 OOE 

03 

0.328E 

03 

0.298E  03 

0  •  U  9E 

04 

0*1 50E 

03 

0.358E 

03 

0.331 E  03 

0.U9E 

04 

0.250E 

03 

0.40  I  E 

03 

0.376E  03 

0*1 20E 

04 

0.350E 

03 

0.432E 

03 

0.410E  03 

0* 120E 

04 

0.500E 

03 

0.467E 

03 

0.449E  03 

0. 120E 

04 

0.700E 

03 

0.503E 

03 

0.491E  03 

0.1 20E 

04 

0. 10QE 

04 

0.545E 

03 

0.542E  03 

0*1 I9E 

04 

0.1 50E 

04 

0.596E 

03 

0.612E  03 

0*1 I7E 

04 

0.200E 

04 

0.636E 

03 

0.672E  03 

0*1 14E 

04 

0.2S0E 

04 

0.66SE 

03 

0.731  E  03 

0« 1 09E 

04 

0 *290E 

04 

0.S90E 

03 

0.78SE  03 

0  *  1 04E 

04 

0.310E 

04 

0.700E 

03 

0.825E  03 

0.993E 

03 
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TABLE  VIII  (Concluded) 


0.321E 

04 

0.705E 

03 

0.903E  03 

0.903E  03 

OUTPUT 

GFLOU/ 

0.041  E 

00 

S 

QEV 

QTOT 

QRAT 

0.250E- 

01 

0  •  1 39E 

00 

0.223E-03 

0.265E-03 

TG 

PG 

EXP9 

0.250E 

04 

0.449E 

03 

0.589E  00 

S 

QEV 

QTOT 

QRAT 

0.750E- 

01 

0.1 89E 

00 

0.51  BE— 03 

0.61 6E-03 

TG 

PG 

EXP9 

0.250E 

04 

0 .448E 

03 

0.759E  00 

S 

QEV 

QTOT 

QRAT 

O.  125E 

00 

0.1 83E 

00 

0  .794E-03 

0.944E-03 

TG 

PG 

EXP9 

0.250E 

04 

O.440E 

03 

0.693E  00 

S 

QEV 

QTOT 

QRAT 

0.  t  75E 

00 

0.201E 

00 

0.109E-02 

0. 129E-02 

TG 

PG 

EXP9 

0.250E 

04 

0.448E 

03 

0.718E  00 

S 

QEV 

QTOT 

QRAT 

0.225E 

00 

0.212E 

00 

0.139E-02 

0.165E-02 

TG 

PG 

EXP9 

0.250E 

04 

0.448E 

03 

0.711E  00 

5 

QEV 

QTOT 

QRAT 

0.275E 

00 

0.227E 

00 

0.170E-02 

0.202E-02 

TG 

PG 

EXP9 

0.250E 

04 

0.447E 

03 

0.713E  00 

S 

QEV 

QTOT 

QRAT 

0.32SE 

00 

0.242E 

00 

0  .20 1 E-02 

0.239E-02 

TG 

PG 

EXP9 

0.250E 

04 

0.447E 

03 

0.713E  00 

S 

QEV 

QTOT 

QRAT 

0.375E 

00 

0.260E 

00 

0.234E-02 

0.279E-02 

TG 

PG- 

EXP9 

0.250E 

04 

0.446E 

03 

0.714E  00 

5 

QEV 

QTOT 

GRAT 

0.425E 

00 

0.280E 

00 

0.268E-02 

0.319E-02 

TG 

PG 

EXP9 

0.25CE 

04 

0.446E 

03 

0.715E  00 

S 

QEV 

QTOT 

QRAT 

0.475E 

00 

0.302E 

00 

0.303E-02 

0.361E-02 

TG 

PG 

EXP9 

0.230E 

04 

0.445E 

03 

0.716E  00 

S 

QEV 

QTOT 

QRAT 

0.525E 

00 

0.327E 

00 

0.340E-02 

0.404E-02 
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TABLE  IX 

SAMPLE  CALCULATION  FOR  LIQUID  FILM  COOLING 


INPUT  DATA 

0 

•OOOE-38 

0. OOOE-38 

0.200E-03 

0 

•  OOOE 

-38 

0. OOOE-38 

0 

•  100E 

01 

0. OOOE-38 

0 .C00E-38 

0 

.OOOE 

-38 

0. OOOE-38 

0 

•  375E 

00 

0 • OOOE— 38 

0. 000E-38 

0 

.OOOE 

-38 

0. OOOE-38 

0 

•  1  HE- 

■04 

0  o  850E 

-08 

0 .OOOE-38 

0 

•  OOOE 

-38 

0. OOOE- 38 

0 

•  600E 

02 

0 • OOOE— 38 

0. OOOE-38 

0 

.OOOE 

-38 

0. OOOE-38 

0 

•900E- 

■02 

0.300E 

-04 

0. 000E-38 

0 

•  OOOE 

-38 

0.000E-38 

0 

. 133E- 

■01 

0. 170E 

-04 

0. OOOE-38 

0 

•  OOOE 

-38 

0  .OOOE-38 

0 

•  70  OE • 

00 

0 • OOOE— 38 

0. OOOE-38 

0 

•OOOE-38 

0  .OOOE-38 

0 

•  300E 

03 

0 • 400E 

04 

0.480E  00 

0 

•  500E 

-02 

0.700E  02 

0 

.  140E 

01 

0 • 400E 

00 

0. 1 10E  01 

0 

•  533E 

02 

0.702E  01 

0 

•  1  90E 

03 

0. 180E 

02 

0.200E  01 

0 

•  1  OOE 

Ol 

0.1 OOE  01 

0 

• 30CE-0 1 

0.1 13E 

01 

0.2S0E  00 

0 

•  858E 

02 

120  10 

10  1 

20 

0 

.OOOE-38 

0. 1S0E 

01 

0.400E-01 

0 

.  500E 

04 

0 

•  1  OOE 

01 

0.590E 

00 

0.400E-01 

0 

.  SOOE 

04 

0 

•  1  03E 

01 

0.550E 

00 

0.400E-01 

0 

•  SOOE 

04 

0 

.  105E 

01 

0.520E 

00 

0.400E-01 

0 

.  5oOE 

04 

0 

.  10QE 

01 

0.495E 

00 

0.400E-01 

0 

•  500E 

04 

0 

•  1  10E 

01 

0.485E 

00 

0.400E-01 

0 

.SOOE 

04 

0 

.  1  13E 

01 

0.480E 

00 

0.400E-01 

0 

•  500E 

04 

0 

•  1  15E 

01 

0.485E 

00 

0.400E-01 

0 

•  SOOE 

04 

0 

•  II7E 

01 

0  »495E 

00 

0.400E-01 

0 

•  500E 

04 

0 

•  31  3E 

01 

0.200E 

01 

0.160E  00 

0 

•  500E 

04 

0 

•  200E 

00 

0 • 53 1 E 

02 

0*212E  02 

0 

.  1 08E 

04 

0 

•  500E 

00 

0.796E 

02 

0.476E  02 

0 

.  L  I  OE 

04 

0 

•  100E 

0  1 

0. 102E 

03 

0.697E  02 

0 

.  1 1  IE 

04 

0 

•  200E 

01 

0. 126E 

03 

0.940E  02 

0 

•  I  12E 

04 

0 

•  400E 

01 

0.I53E 

03 

0.121E  03 

0 

*1  13E 

04 

0 

•  70CE 

01 

0*1 77E 

03 

0.145E  03 

0 

.  1 14E 

04 

0 

•  I20E 

02 

0.202E 

03 

0.170E  03 

0 

.1  15E 

04 

0 

•  20  0E 

02 

0.22BE 

03 

0.1 96E  03 

0 

*  1  I6E 

04 

0 

•  40  OE 

02 

0.267E 

03 

0.236E  03 

0 

*  1I7E 

04 

0 

•  700E 

02 

0.303E 

03 

0.273E  03 

0 

•  1  18E 

04 

0 

•  100E 

03 

0.328E 

03 

0.298E  03 

0 

•  1  I9E 

04 

0 

•  1  50E 

03 

0 • 358E 

03 

0.331E  03 

0 

.  I  19E 

04 

0 

•  250E 

03 

0  «  40  1  £ 

03 

0.376E  03 

0 

•  120E 

04 

0 

•  350E 

03 

0 * 432E 

03 

0.410E  03 

0 

•  I  20E 

04 

0 

•  SOOE 

03 

0  o 467E 

03 

0.449E  03 

0 

.  120E 

04 

0 

•  70  OE 

03 

0 • 503E 

03 

0.491E  03 

0 

•  120E 

04 

0 

•  100E 

04 

0 • 545E 

03 

0.542E  03 

0 

.  1  19E 

04 

0 

«  1 50E 

04 

0.596E 

03 

0.612E  03 

0 

•  l  I7E 

04 

0 

•  200E 

04 

0.636E 

03 

0.672E  03 

0 

•  1  14E 

04 

0 

•  2S0E 

04 

0.668E 

03 

0.731 E  03 

0 

.  109E 

04 
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TABLE  IX  (Concluded) 


OUTPUT  (SELECTED) 

QRAT 

or 

GFLOW 

0.0 I l 1 4235 

0.00500000  0 

.35917528 

S 

Q/A 

TB 

NBOIL 

0. 150E- 

-01 

0.195E  04 

0.703E  02 

0 

H 

EXP9 

TF2  <  N I  ) 

N 1 

0. 1  15E 

03 

0. I 1 1 E  01 

0  a 860E  02 

1 

S 

Q/A 

TB 

NBOIL 

0.21BE 

00 

0*1 16E  06 

0.907E  02 

H 

EXP9 

TF2  <  N 1  ) 

Nl 

0. 124E 

03 

0.100E  01 

0.1 19E  03 

1 

S 

Q/A 

TB 

NBOIL 

0.41 5£ 

CO 

0.324E  06 

0.1 37E  03 

0 

H 

EXP9 

TF2 (N [ ) 

NI 

0. 179E 

03 

0  o 1 OOE  01 

0.1 93E  03 

1 

S 

Q/A 

TB 

NBOIL 

0.61 5E 

00 

0.474E  06 

0.1 67E  03 

0 

H 

EXP9 

TF2(N[ l 

NI 

0.206E 

03 

0.1 OOE  01 

0.248E  03 

2 

S 

Q/A 

TB 

NBOIL 

o.aisE 

00 

0.597E  06 

0.191E  03 

0 

H 

EXP9 

TF2  f  N I  > 

NI 

0.487£ 

03 

0.1  OOE  01 

0.281 E  03 

2 

S 

Q/A 

TB 

NBOIL 

0. 102E 

01 

0.671E  06 

0.206E  03 

0 

H 

EXP9 

TF2<NI ) 

NI 

0.712E 

03 

0.1 OOE  01 

0 .30! E  03 

3 

S 

Q/A 

TB 

NBOIL 

0.122E 

01 

0  «  739E  06 

0.209E  03 

0 

H 

EXP9 

TF2  (  N 1  ) 

NI 

0.703E 

03 

0 o 903E  00 

0.313E  03 

3 

S 

Q/A 

TB 

NBOIL 

O.  152E 

01 

0.760E  06 

0.21  IE  03 

1 

H 

EXP9 

TF2 (NI ) 

NI 

0.391 E 

03 

0  c  877E  00 

0.3I4E  03 

4 
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APPENDIX  G 

DEVELOPMENT  OF  DIFFERENCE  EQUATIONS 
FOR  GAS  FILM  COOLING  ANALYSIS 

The  difference  equation  for  the  momentum  (or  force)  balance , 
Gq.  (6),  page  14,  is  developed  by  referring  to  Pig.  2  and  evaluating 
the  terms  in  Eq.  (5). 

The  pressure  force  is  the  difference  in  pressure  multiplied 
by  the  annular  area  of  the  layer  projected  in  the  axial  direction, 
thus 


Pressure  forces  in 

d  direction  on  «  ird{  Sj  (pj  -  p2)  (86) 

coolant  layer  1 

The  inertia  forces  are  the  net  rate  of  momentum  flow  into 
the  volume  element  as  a  result  of  the  axial  mass  flow  rate  and  the 
velocity  of  the  fluid.  Thus, 


Inertia  forces  in 
A  direction  on 
coolant  layer  1 


w{ 

"  —  (V'U  “  V'2T) 

BC 


(87) 


The  shear  forces  on  the  top  and  bottom  of  the  coolant  layer 
are  calculated  from  the  shear  stress  equation 

dv 


pe 


dy 


(88) 
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Considering  t  as  the  current  of  momentum  and  dv  as  the  potential  for 
momentum  flow,  Eq.  (88)  may  be  .written  as 


dv 

T  "  (dy/pe) 


(39) 


By  analogy,  the  resistance  to  momentum  flow  is  given  by  the 
denominator  of  Eq.  (89).  Since  resistances  in  series  are  additive 
the  shear  stress  on  the  top  of  the  first  coolant  layer  is  given  In 
finite  difference  form  as 


vii  +  v2i  v'n  +  V‘*21 

- 2 - y - 

Ttop  df 

2(^61  + 


(90) 


Evaluation  of  the  shear  stress  on  the  bottom  of  coolant  layer  1  is 
determined  in  a  similar  manner.  The  shear  forces  are  found  by 
multiplying  the  stresses  by  appropriate  areas. 


Shear  forces  on  the 
top  and  bottom  of 
coolant  layer  1 


irdf  A  (vi i  +  V2i  -  v'2i  -  v'ij) 


Pl«l 


JL 

pici 


rdj  A  (v"i2  +  v'22  ”  v>ll  ~  v"2i) 

*“  jT  ~jT 

Pi^i  P2e2 


(91) 
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Equating  the  sum  of  Eqs.  (86) (87) (91)  to  zero  yields  the  momentum 
balance  equation  In  finite  difference  form. 

The  difference  equation  for  the  steady  state  energy  balance , 

Eq.  (22),  page  20,  is  developed  in  a  similar  manner.  Refer  to  Fig.  2, 
•page  16. 

The  enthalpy  change  in  the  axial  direction  of  a  layer  is 
determined  from  the  ideal  gaB  laws  as 

Enthalpy  change  in 

A  direction  of  ■  wfc^T7!^  -  T/2^)  (92) 

coolant  layer  1 

The  change  in  kinetic  energy  is  determined  from  the  square 
of  the  velocities  as 

Kinetic  energy  change 

W1  2  2 

in  A  direction  of  «  2Jg”  "  v'2l  )  (93) 

coolant  layer  1 

The  thermal  energy  conduction  term  is  calculated  using  the 
temperature  difference  between  layer.  The  resistance  to  the  flow  of 
heat  is 

A 

Thermal  resistance  *  - r 

pcEijr 
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Thus 


Thermal  energy  transfer 
from  gas  layer  1  to  ■  ■ 

coolant  layer  l  by  diffusion 


Tlx  +  T2X 

(-*.1 

1 


+ 


■l  +  1 

— 2 - — ) 

_!L  T..| 

2pfci Ei^y 


(94) 


and  a  similar  equation  Is  written  for  the  thermal  energy  transfer 
between  coolant  layer  1  and  coolant  layer  2. 

In  evaluating  the  work  at  the  boundaries  of  the  Increment 
only  shear  stress  effects  are  considered. 

Thus 


Viscous  work/unit  volume 


,dVv 

“ld?> 


2 


(95) 


Apply  Eq.  (95)  to  the  first  coolant  layer,  noting  that  the 
velocity  gradient  In  difference  terms  is  established  using  three 
layers 


work/ unit  volume 


2 

Vil  +  V2l  v'l2  +  v'2z 

- 2 — - 2 - 

6!  5f  S£ 

T  +  +  “T 


Thus  evaluating  y  and  multiplying  by  the  layer  volume  of  coolant 
layer  1  gives 


work  at  boundaries 
on  coolant  layer  1 


TfdfAfijfpie^vii  +  v2i  -  v'22) 

J8c(2$i  +  $2 +  fiTP 


(96) 


from  shear  stress  effects 
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